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Q_i ABSTRACT linear modification of the Schrodinger equa- 
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Quantum mechanics is an extremely sue- tion, which resolves these problems, while 

o3 

cessful theory that agrees with every ex- giving the same experimental results as quan- 
ta" 1 

periment. However, the principle of linear turn theory in the microscopic regime. Two 

superposition, a central tenet of the the- underlying theories for this phenomenology 

ory, apparently contradicts a commonplace are reviewed: Trace Dynamics, and grav- 
en 

observation: macroscopic objects are never ity induced collapse. As one approaches the 
found in a linear superposition of position macroscopic scale, predictions of this pro- 
states. Moreover, the theory does not ex- posal begin to differ appreciably from those of 
£> plain why during a quantum measurement, quantum theory, and are being confronted by 

deterministic evolution is replaced by prob- ongoing laboratory experiments that include 
abilistic evolution, whose random outcomes molecular interferometry and optomechanics. 
obey the Born probability rule. In this article These experiments, which test the validity of 
we review an experimentally falsifiable phe- linear superposition for large systems, are re- 
nomenological proposal, known as Continu- viewed here, and their technical challenges, 
ous Spontaneous Collapse: a stochastic non- current results, and future prospects summa- 
rized. We conclude that it is likely that over 
the next two decades or so, these experiments 
can verify or rule out the proposed stochastic 
modification of quantum theory. 
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I. INTRODUCTION 

Quantum theory has been extremely suc- 
cessful in explaining results of experiments, 
ranging from the spectrum of black-body ra- 
diation, atomic spectra, molecular chemistry, 
atomic interferometry, quantum electrody- 
namics and nuclear physics, to properties 
of lasers, superconductivity, semi-conductor 
physics, Bose-Einstein condensation, Joseph- 
son junctions, nano-technology, applications 
in cosmology, and much more. The theory 



is not contradicted by any experiment. Yet, 
there is one apparently innocuous observed 
phenomenon the theory seems unable to ex- 
plain, and in fact seems to contradict with. 
This is the observed absence of superposi- 
tion of different position states in a macro- 
scopic system. Quantum theory, by virtue of 
the principle of linear superposition, predicts 
that a microscopic object such as the electron 
can be in a superposition of different posi- 
tions at the same time, and this is of course 
observed, for example in the famous double- 
slit interference experiment. Moreover, the 
theory in principle makes no distinction be- 
tween microscopic objects and macroscopic 
ones, and predicts that large objects can also 
be in more than one place at the same time. 
But this is not what we observe. A table 
for example, unlike the electron, is never ob- 
served to be 'here' and 'there' simultaneously. 

Why should this be so? The present re- 
view article is devoted to discussing one pos- 
sible proposed resolution, known as Models of 
Spontaneous Wave Function Collapse, which 
is experimentally falsifiable. Namely that, al- 
though quantum theory is extremely success- 
ful in the microscopic domain, it is an ap- 
proximation to a more general theory. This 
general theory is capable of explaining the 
absence of macroscopic superpositions. It 
goes over to quantum mechanics in the mi- 
croscopic limit, and to classical mechanics in 



the macroscopic limit, but differs from both where S is the action of the system (Landau 



quantum and classical mechanics in the in- and Lifshitz, 1976). 

termediate [mesoscopic] regime which marks In con trast, the quantum dynamics of this 

the transition from the micro- to the macro- sys tem is described by first converting the 

world. A large number of experiments world- q . and p . to operators qi, Pi satisfying the 

wide are operating or are being planned, to commutation relations [q i; p;] = ih and then 

test the validity of linear superposition in the proposing that the operators evolve via the 

mesoscopic domain, and in this article we will Heisenberg equations of motion 

review the proposed modification to quantum % r , % r 

qi = -r[q*,H], Pl = --[ Pi ,H]. (4) 

mechanics, and the laboratory experiments 11 11 

i.i f i -r ,i- i Quantum dynamics is equivalently described 

which can falsity this proposal. ^ J ^ J 

by the time evolution of the system's wave- 
function if), which is a normalized element 
of a Hilbert space and obeys the norm- 
preserving Schrodinger equation 

ih^ = Hip, J dqiP*iP = 1 . (5) 

In the Heisenberg picture the relation be- 
tween quantum and classical mechanics is 
expressed by replacing operators by ordi- 
nary functions, and the commutators in the 



A. The relation between non- 
relativistic quantum mechanics and 
classical mechanics 

The classical dynamics of a system of par- 
ticles having a Hamiltonian H is described in 
phase space (qi,Pi) by Hamilton's equations 
of motion 



. _ OH . OH 
Qi — ~Q^/> Pi ~ ~~Qq~. (• ' equations of motion by Poisson brackets. A 

or via Poisson brackets more insightful comparison is obtained in the 

Schrodinger picture, and for the purpose of il- 

qt = {qi, H} , pi = {pi, Hj . (2) 

lustration it is adequate to consider the case 

The state of the system at an initial time tn c ■ i , . , { ■ 

J u of a single particle of mass m moving m one 

is a point in the phase space, and the equa- i- • ( u-^+ucu-j- 

v fry dimension, tor which the bchrodmger equa- 

tions of motion determine the location of the , • u , • ■,■ 

tion can be written m the position represen- 

system point at a later time. An equivalent tation> &fter defining ^ ^ js/h^ as 

description of the dynamics is through the qc I / dS\ 2 ih d 2 S 

" + V (l)-7r-^- (6) 



Hamilton- Jacobi equation dt 2m \dq J 2m dq 2 

qo / QQ\ m the approximation in which the last 



dt ~ V 9qJ K > terming can be neglected, this equation re- 
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duces to the classical Hamilton- Jacobi equa- classical solution, quantum mechanics pre- 
tion ([TJ) diets that the sum of these two wave-packets 



dS 1 f dS\ 2 is also a solution, and in principle such so- 

\ y / lutions should be observed in nature. How- 
provided the quantity S is assumed to be real ,. J , . , , , 

ever, according to classical mechanics, such a 

and identified with the action of the system. ., . . , , . . 

superposition is not a solution of the equa- 

This essentially corresponds to the limit S 3> . c • -x 1 j • ^ 

tions of motion, nor is it observed m the 

h. [We will not consider the more precise . n , , AT . , 

L macroscopic world around us. JNaively, we 

treatment where S is separated into its real , n . . n . , , , . , 

believe that classical mechanics, which ap- 

and imaginary parts, as it is not crucial to ,. ,. ... 

plies to macroscopic systems, is a limiting 

the present discussion! . . , , , 

1 case of quantum mechanics, and hence quan- 

There is thus a well-defined sense in which turn mechanics should apply to large systems 
the Schrodinger equation goes over to the as well. Why then we do not observe macro- 
Hamilton- Jacobi equation in the limit, and a scopic superpositions [such as a table being 
description of the dynamics in Hilbert space 'here' and 'there' at the same time]? 

gets replaced by a description in terms of evo- • , , , , , , , , , u 

& 1 J y One might argue that even though the 

lution of position and momentum coordinates Hamilton- Jacobi equation is non-linear, its 

in phase space. Yet, there is a profound as- non . linearity cannot be used to deduce the 

pect which gets lost in the limiting process. u a u c 

^ to to ^ observed absence of macroscopic superpo- 

The Schrddinger equation is linear: if ^ and sitionS) bec&use ^ theQry ig af _ 

^ 2 are two solutions of @ then the linear ^ aR an approximation . The last term in 

superposition c^i + c 2 ^ 2 is also a solution, £qn _ g ? howsoever smaU) is a i ways non . 

where Ci and c 2 are complex coefficients. On , , , , , , , 

1 z 1 zero and present, and can be used to trans- 

the other hand the Hamilton- Jacobi equation form back tQ ^ Hnear Schr5dinger equ ation. 

Q is non-linear : if Si is a solution corre- At & fundamental level; the descrip ti on of 

spending to one space-time trajectory, and S 2 ^ dynamicS; even for a macroscopic das . 

is a solution corresponding to another space- gical object; ig ^ termg of ^ wave . function 

time trajectory, then clearly a.S, + a 2 S 2 is of ^ quantum gtate; and not in terms of 

not a solution of this equation. ,. , . , . tt . 1+ 

H the action which appears m the Hamilton- 

In particular, if ip\ is a wave-packet which Jacobi equation. Hence superpositions must 

is peaked around one classical solution and be there. Nonetheless, one is left with the dis- 

ip2 is a wave-packet peaked around another comforting feeling that the prediction of the 



Hamilton- Jacob i equation regarding position observable O. Immediately after interaction, 
superpositions seems to be at variance with the combined state of the system and appa- 
quantum theory, and in accord with what ratus is 

is actually observed. Thus one needs to ex- , , A , , fQ \ 

V> = CylpyAy + c 2 ip 2 A 2 (8) 

plain the following^] : why is it that macro- 

where c\ and c 2 are complex coefficients pro- 

scopic objects which obey the rules of classi- 

portional to the relative amplitudes for the 

cal mechanics are not found in superposition 

system to be in the two states ipy and ijj 2 . 

of different position states, in spite of quan- 
tum theory suggesting otherwise? There is According to quantum mechanics, this 
no unique universally accepted answer to this state ip of Eqn. @ should evolve linearly 
question. In this sense this is an unsolved by way of Schrodinger evolution, and the lin- 
problem. ear superposition of the two parts should be 

The absence of macroscopic superposi- preserved. But that is not what is observed 
tions is of course at the heart of the so-called in a quantum measurement. The outcome of 
quantum measurement problem (lBassi andl the measurement is either pointer at position 



Ghirardil [2000} . Suppose a quantum system M (and hence system is driven to state ^i) 



which is in a superposition of two eigenstates or pointer at position A 2 (system is driven 

ipi and ip 2 of a physical observable O interacts to state ip 2 ). Repeated measurements on the 

with a classical measuring apparatus A. Let same initial quantum state yield outcome ipy 

us say that the state ipy of the quantum sys- or ip 2 with relative probability |ci| 2 : |c 2 | 2 . 

tern corresponds to a pointer position state This is the Born probability rule. The pro- 

Ai of the apparatus [meaning that if the sys- cess of measurement destroys linear superpo- 

tem had been in the state ipy and interacted sition of the initial states ipy and ip 2 . This 

with the apparatus, the pointer would result indeed has to do with the fact that the appa- 

in the position Ay and we would interpret ratus [which is a macroscopic object] is never 

that the observable had the value Oy). Simi- simultaneously observed in a linear superpo- 

larly, the pointer position A 2 corresponds to sition of pointer position states Ay and A 2 . 

the system state i/j 2 and a value 2 for the To the extent that we do not understand why 

~ . . . macroscopic objects are not found in super- 

One should keep in mind the difference between 

the conceptual issue raised here, and the purely posed states, we do not understand why the 

technical fact that performinff an experiment which , u i -j.- 

H & ^ measurement process breaks superposition. 

tests macroscopic superpositions (table "here" + 

table "there") is practically unfeasible. Perhaps even more remarkable is the 
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emergence of probabilities. The Schrodinger 
evolution is deterministic, and so is the clas- 
sical evolution according to the Hamilton- 
Jacobi equation. In our discussion above, on 
the transition from the Schrodinger equation 
to the Hamilton- Jacobi equation, nowhere 
did we encounter probabilities. And for good 
reason. Because the initial state is always ex- 
actly specified [including at the start of the 
measurement process, as in Eqn. (|8])], unlike 
in classical probability theory, where prob- 
abilities arise because of uncertainty in our 
knowledge of the initial state of the system. 
Thus the status of probabilities in quantum 
theory is absolutely unique, and besides ex- 
plaining the absence of macroscopic superpo- 
sitions one must also explain why during a 
measurement probabilities arise, in violation 
of deterministic linear superposition, and the 
quantum system is driven to one or the other 
outcome in accordance with the Born rule. 

Another important and related unsolved 
problem is the following: when do we call a 
physical system a quantum system and when 
we do call it a classical measuring appara- 
tus? In other words, where is the quantum- 
classical divide? How much mass, or degrees 
of freedom [say number of nucleons] should 
an object have, before it qualifies to be an 
apparatus? Of course, in order for it to be 
called an apparatus, different pointer posi- 
tions should never be simultaneously real- 

7 



ized, but one does not know at what mass 
scale this transition from micro to macro 
[and the concurrent breakdown of superposi- 
tion] takes place. Interferometry experiments 
have shown that quantum theory, and hence 
linear superposition, holds for molecules at 
least as large as those having about a thou- 
sand atoms [hence a molecular mass of 1CT 21 
grams]. Efforts are afoot to push this test 
limit up to objects of about a million atoms 
[1CT 18 grams]. On the other end, classical 
behavior [absence of superpositions of states 
corresponding to different positions] is known 
to hold out down to about a microgram [10 18 
atoms]. There is thus an enormous desert 
of some fifteen orders of magnitude, where 
linear quantum superposition yet remains to 
be tested experimentally. Does quantum me- 
chanics hold at all scales, including macro- 
scopic scales, and is there a way to under- 
stand the absence of macroscopic superposi- 
tions while staying within the framework of 
quantum theory? Or is it that somewhere in 
the grand desert modifications to quantum 
theory start becoming significant, so that lin- 
ear superposition becomes more and more of 
an approximate principle as the size of a sys- 
tem is increased, until for large objects the 
superposition of states corresponding to dif- 
ferent positions is no longer a valid principle? 
What exactly is the nature of the quantum- 
to-classical transition? A large number of 



ongoing and planned experiments worldwide 
are gearing up to address this question. 

The paradoxical issue of deterministic 
evolution followed by a peculiar probabilis- 
tic evolution during a measurement was of 
course well-appreciated by the founding fa- 
thers of quantum mechanics. Over the last 
eighty-five years or so since the discovery of 
the Schrodinger equation, extraordinary the- 
oretical effort has been invested in trying to 
find an answer to what is generally known as 
the measurement problem in quantum me- 



function collapse, such as the EPR paradox, 
whether this 'violates the spirit of relativity' 
or whether there is a need for a radical change 
in our ideas about space-time structure. 

B. Proposed resolutions for the quan- 
tum measurement problem and for the ob- 
served absence of macroscopic superposi- 
tions 

1. The Copenhagen Interpretation 



chanics QAlbert[ |1992| |Bell[ |1987b[ |Ghirardi[ The Copenhagen interpretation flBohr 



2005 




Leggett 


2002 




2005 




Maudlin 




2011 



1928) [reprinted in (Wheeler and Zurek 



Wheeler and Zurek, 1983). In the next sub 



section we give a brief overview of a few of the 
major categories of the explanations, keeping 
in mind that the modern outlook is to discuss 
this problem not in isolation, but in conjunc- 
tion with the question of lack of macroscopic 
superpositions, and SIS db part of the much 
broader investigation of the exact nature of 
the quantum-to-classical transition. 

Our review of the measurement problem 
will be almost exclusively confined to the con- 
text of non-relativistic quantum mechanics, 
as the relativistic version seems not within 
reach at the moment (though it does not 
seem symptomatic of a deep incompatibil- 
ity of modified quantum mechanics and rel- 
ativity). Thus we will not discuss issues 
raised by the instantaneous nature of wave- 



1983)] postulates an artificial divide between 
the micro-world and the macro-world, with- 
out quantitatively specifying at what mass 
scale the divide should be. Microscopic ob- 
jects obey the rules of quantum theory [su- 
perposition holds] and macroscopic objects 
obey the rules of classical mechanics [super- 
position does not hold]. During a measure- 
ment, when a micro-system interacts with 
a macro-system, the wave-function of the 
micro-system 'collapses' from being in a su- 
perposition of the eigenstates of the measured 
observable, to being in just one of the eigen- 
states. This collapse is postulated to happen 
in accordance with the Born probability rule, 
and no dynamical mechanism is specified to 
explain how the collapse takes place, in ap- 
parent contradiction with the linearity of the 
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Schrodinger equation. level one should prescribe a physical mech- 

von Neumann gave a more precise form to anism which causes the so-called collapse of 

this interpretation by explicitly stating that the wave-function. 

evolution in quantum theory takes place in The Copenhagen interpretation does not 

two ways : (i) deterministic evolution ac- solve the quantum measurement problem, 

cording to the Schrodinger equation before nor does it explain the absence of macro- 

a measurement, as well as after a measure- scopic superpositions. 

ment, and (ii) non-deterministic, probabilis- The "histories" approach is an observer in- 

tic evolution [the projection postulate] during dependent generalization of the Copenhagen 



a measurement (von Neumann, 1955). interpretation wherein the notions of appa- 

At a pragmatic level, this can be taken ratus and measurement are replaced by the 

to be a perfectly valid set of rules, in so far more precise concept of histories. In this 

as the goal is to apply quantum theory to approach, the reduction of the state vector 

sub-atomic, atomic and molecular systems, appears as a Bayesian statistical rule for re- 

and to compare with experiments the predic- la ting the density matrix after measurement 

tions based on theoretical calculations. How- to the density matrix before measurement 



(Griffiths 


2002 


Hartle 


1992 


Omnes , 


1992 



tions raised in the previous section, by sim- 1994, 1999 ). 

ply raising unresolved issues to the level of 

postulates. The interpretation creates an ill- „ n , 

c r 2. Veconerence 

defined micro-macro separation, which is be- 
ing challenged by modern experiments which The phenomenon of decoherence, which 
are verifying superposition for ever larger sys- is observed in laboratory experiments, high- 
tems. There is no precise definition as to lights the role played by the environment 
which systems classify to serve as a 'classical when a quantum system interacts with a 
measuring apparatus'. Even though there is measuring apparatus, during the process of 
a sense in which the Hamilton- Jacob i equa- measurement. By environment is meant the 
tion is a limit of the Schrodinger equation, system of particles which surround the appa- 
no attempt is made to explain the apparently ratus. More precisely one could define the en- 
different predictions of the two theories with vironment as the collection of particles which 
regard to absence of macroscopically different is present within a radius cT of the appara- 
position superpositions. At a fundamental tus, where T is the duration of the measure- 



ment: these are hence particles which can Decoherence destroys interference amongst 



causally interact with and influence the ap- alternatives, which is what Eqn. (12) signi 



paratus during a measurement. fies, but because it operates within the frame- 
To illustrate the effect of decoherence, we work of linear quantum mechanics, it cannot 
assume that the system on which measure- destroy superposition. Since loss of super- 
ment is to be made is a two state system ini- position is what is seen during a measure- 
tially in the state ment, decoherence does not explain the mea- 

j,(t = 0) = dV>i + c 2 ^ 2 . (9) 



surement process. What (12) implies is that 
decoherence forces quantum probability dis- 
tributions to appear like classical probabili- 
ties (weighted sums of alternatives); however 
this is neither necessary nor sufficient to ex- 
plain the outcome of an individual measure- 
= i>(t = 0) (fioA 4>oe- (10) ment. [The issue of non-observation of su- 

Over time, as a result of interaction, this Perposition for macro-systems becomes even 
state evolves into the state more acute in the case of isolated systems 



Denoting the initial state of the apparatus by 
4>oa and the initial state of the environment 
by 4>oe, we can write the net initial state as 
the direct product 



$(t) = c 1 tp 1 (p(t) EA1 + C 2 l\) 2 <\){t) EA2 . (11) 



as there are no environment degrees of free- 
dom to be traced out. Then the theory seems 
Here, <f,{t) EA1 and <j){t) EA2 denote macroscop- unable to explain the breakdown of supe rpo- 

lcally distinguishable entangled states of the sition for i so i a ted macroscopic systems, such 
apparatus and the environment. ag the universe as a whole i 



As demonstrated below, during measure- 
ment, the process of decoherence operates in 
such a way that very quickly, the inner prod- 
uct 

(<P(t) EAl \<p(t) EA2 ) ^0 (12) 
starting from the value unity at t — 0. The 
final state is reduced to a statistical mixture 
of states with relative weights |ci| 2 : |c 2 | 2 . 



The loss of interference can be understood 
as a consequence of the interaction of the 
very large number of particles of the envi- 
ronment with the apparatus. Assuming that 
the measurement starts at t = 0, the prod- 
uct {(pit) EA i\(p{t) EA2 ) , which is one at t — 0, 
rapidly goes to zero. To see this, one notes 

that in general a particle of the environment, 

This by itself does not explain why during ,. , .,, , , , 

J r j a sa y ^ e p ar ticle, will be scattered by the 

i measurement , , A f , , c ■, . , i- e 

state A\ of the apparatus to a final state dit- 

$(t) — > ipi4>(t) EA i or $(t) — >■ ip 2 4>(t) EA2 . ferent from the one to which it will be scat- 
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tered from the apparatus state A 2 . Thus, the 
product 

(E^EJit) = II, i(E(t = 0)\S A1 S A2 \E(t = 0)>, 

is made up of an ever increasing number of 
quantities, each of which is smaller than one, 
Sai and Sa2 being scattering matrices de- 
scribing the action of the apparatus on the 
environment. Hence this product can be 
written as exp(— At) and goes to zero for large 
t, A being the decoherence rate. Because the 
environment has a very large number of par- 
ticles, this cross-product between the two en- 
vironment states is very rapidly suppressed, 
and is responsible for the emergence of the 



property described by Eqn. (12). The deco 



herence time-scale A 1 is much smaller than 
the duration T of the measurement. 

The above discussion is partly based on 



the article by Adler (Adler, 2003) where a 



more detailed description of decoherence in 
the context of measurement can be found. 
There is a vast literature on decoherence, 
including the experiments and models by 



3. Many-Worlds interpretation 



The Many-worlds interpretation was in- 



vented by Everett (Everett III, 1957) to 



(Brune et al. 1996; Gerlich et al, 2007 Har- 



ris and Stodolsky, 1981), books by (Breuer 



and Petruccione |2000 Joos et al. 2003 



Schlosshauer, 2007) and the seminal papers 



(Caldeira and Leggett, 1981 Joos and Zeh 



1985 Zeh, 1970) and reviews (Schlosshauer 



2005 Vacchini and Hornberger] 2009 Zurek 



1991 2003), and (Bacciagaluppi, 2007) 



counter the Copenhagen interpretation. Ac- 
cording to Everett, evolution during a mea- 
surement is also Schrodinger evolution, and 
there is no such thing as a non-deterministic 
probabilistic evolution during the measure- 
ment. Thus in this interpretation the state 
^ evolves during a measurement according 
to the Schrodinger equation. Why then does 
it appear as if only one of the two outcomes 
has been realized? The answer is that the 
state continues to be of the form 

^ = cifaAiOx + c 2 ip 2 A 2 2 (13) 

where 0\ (0 2 ) is the state of the observer 
where the observer detects the system and 
apparatus in state one (two). The two parts 
of this state exist in two different branches of 
the Universe. 

Despite appearances, there is no logical in- 
consistency in the above interpretation, as 
has been argued by Everett : it is merely 
the assertion that Schrodinger evolution is 
universally valid at all scales, and the break- 
down of superposition during a measurement 
is only apparent, not real. The hard part 
is to explain the origin of probabilities and 
the Born probability rule. If the evolution 
is deterministic through and through, why 
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should there be a definite probability asso- 
ciated with an outcome? In our opinion, de- 
spite extensive investigation, this problem re- 
mains unsolved in the many-worlds interpre- 



tation (Barrett and Byrne, 2012 Deutsch 



1998 DeWitt and Graham , 1973; Hsu ( 2011 



Kent 1990 Putnam 2005 Saunders et al. 



2010 Tegmark 2007 Vaidman, 2002; Wal- 



lace, 2003) 



4- Decoherence + Many-Worlds 

Decoherence by itself does not solve the 
measurement problem, because it does not 
destroy superposition. However, one could 
propose that the decohered alternatives both 
continue to co-exist in different branches of 
the Universe, in the sense of the many-worlds 
interpretation, and these branches do not in- 
terfere with each other because decoherence 
is operating. While this merger helps both 
the decoherence and the many-worlds pic- 
ture of a measurement, the origin of the Born 
probability rule continues to be not under- 
stood, and as of now is essentially added as 
a postulate. 

This is perhaps today the 'establishment 
view', wherein one believes that one does not 
need to modify quantum theory in order to 
explain measurement. Its major weakness 
though is that it is not experimentally falsi- 
fiable. What experiment can one perform in 



order to find out whether the other branches 
of the many- worlds exist or not? In the ab- 
sence of such an experiment, we have at hand 
another interpretation of the same quantum 
theory, an interpretation which cannot be ex- 
perimentally distinguished from the Copen- 
hagen interpretation. 

For discussions on decoherence in the con- 



text of many- worlds see (Bacciagaluppi 



2001). 



5. Bohmian Mechanics 

Bohmian mechanics is a quantum the- 
ory of particles in motion. The positions 
of the particles of an iV-particle system are 
Qk k = 1,...,N moving in physical space. 
The role of the wave function, being gov- 
erned by Schrodinger equation, is to direct 
the motion of the particles. The theory is 
deterministic, randomness enters like in clas- 
sical mechanics via typicality. It is shown 
that the outcomes in measurement experi- 
ments are governed by Born's statistical rule. 
The equation of motion for the particles is 
given by Vk = dQk/dt where 

dQ k h 



dt rrik 



ImV Qk log ip(Q 1 , Q 2 , Qn, t). 



Bohmian mechanics is a quantum theory in 
which the collapse of the wave function is 
effective, in contrast to collapse models, so 
that macroscopic interference is in principle 
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possible. Predictions of Bohmian mechanics lapse is shown to obey the Born probabil- 

agree with those of orthodox quantum me- ity rule. The universal dynamics is stochas- 

chanics, whenever the latter are unambigu- tic: the outcome of a measurement is ran- 

ous. Bohmian mechanics would be falsified if dom and unpredictable, but the mathemati- 

collapse models were experimentally verified, cal structure of the dynamics is such that re- 

For literature on Bohmian mechanics see peated measurements on an ensemble of iden- 

(Bohm, 1952a|b Bohm and Hiley 1995 Bub, tically prepared quantum systems are shown 

19971 IDiirr and Goldsteinl 120121 IDiirr et al.\ to yield different outcomes, in relative fre- 



1992 Diirr and Teufel, 2009 Holland, 1993). quencies which obey the Born rule. 

The universal dynamics must be non- 

a r. . _ ,r • .. , linear, in order to allow for the breakdown 

6. Quantum theory is an approximation to a ' 

, of superposition during a measurement. Yet, 

more general theory 1 1 ° 

the non-linearity must be extremely negligi- 
It is proposed here that the measurement ble in the microscopic domain, so that the 
problem and the apparent inability of quan- experimentally observed linear superposition 
turn theory to explain the absence of macro- in microscopic quantum systems is repro- 
scopic superpositions are a consequence of duced. The new dynamics must be stochas- 
trying to apply the theory in a domain where tic; but once again, stochasticity must be 
it is not valid. It is proposed that there is a negligible for microscopic systems, so that 
universal dynamics, to which quantum the- the deterministic Schrodinger evolution pre- 
ory and classical mechanics both are approx- vails. Thirdly, as one progresses from mi- 
imations. In the domain of a quantum mea- croscopic to macroscopic systems, the univer- 
surement, the universal dynamics differs from sal dynamics must allow for non- unitary [but 
quantum dynamics, and operates in such a norm-preserving] evolution : this is essential 
way that interaction of the quantum system so that stochastic evolution can cause all but 
with the apparatus causes a collapse of the one outcome to decay exponentially, some- 
wave-function from a superposition to one of thing which would not be permitted during 
the eigenstates. The collapse is a physical, unitary evolution. Again, non-unitarity must 
dynamical process, and hence the universal be utterly negligible for microscopic systems, 
dynamics provides a physical explanation for Thus, the universal dynamics possesses a set 
the ad hoc collapse postulated by the Copen- of parameters, whose effective values are de- 
hagen interpretation. Furthermore, the col- termined for the system under study in such 
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a way that for microscopic systems these pa- 
rameters take values so that the dynamics is 
experimentally indistinguishable from quan- 
tum dynamics. Similarly, for macroscopic 
systems, there is an amplification mechanism 
built in the equations, such that the dynam- 
ics coincides with classical dynamics. For sys- 
tems that are mesoscopic [neither micro nor 
macro] the dynamics differs from both classi- 
cal and quantum, and is hence experimen- 
tally distinguishable from quantum theory. 
The properties of non-linearity, stochastic- 
ity and non-unitarity also ensure position lo- 
calization for macroscopic objects and hence 
dynamically explain the observed absence of 
macroscopic superpositions. It is clear that 
the universal dynamics is not tied to or in- 
vented for explaining just the measurement 
process or absence of macroscopic superpo- 
sitions - these two phenomena just happen 
to be special cases where the new dynamics 
plays a vital role in understanding the phe- 
nomenon. We say that the universal dynam- 
ics, which describes the behaviour of micro, 
meso and macro objects, is intrinsically non- 
linear, stochastic, and non-unitary. 

Over the last two decades or so there has 
been a significant progress in developing phe- 
nomenological models of such a universal dy- 
namics. At the same time, one would like to 
know if there are underlying theoretical rea- 
sons [new symmetry principles for instance] 
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which compel us to consider a generaliza- 
tion of quantum theory of the kind mentioned 
above, thereby lending an inevitability to the 
phenomenological models that have been pro- 
posed. There has been important progress 
on this front too. Thirdly, there have been 
important technological advances which are 
now permitting a host of experiments to be 
carried out to test these phenomenological 
models, and verify their predictions against 
those of quantum theory. Needless to add, 
all these three facets are best described as 
'work currently in progress'. The purpose of 
the present review is to present a state-of- 
the-art description of (i) the phenomenolog- 
ical models for the universal dynamics, (ii) 
the underlying theories, and (iii) ongoing ex- 
periments which aim to test these models and 
theories. It is our hope that a review of this 
nature will further stimulate the cross-talk 
between phenomenologists, theorists and ex- 
perimentalists in this field, thereby helping 
the community to sharply focus on those as- 
pects of phenomenology and experimentation 
which might be most directly accessible and 
feasible in the near future. 

Phenomenological models of modified 
quantum mechanics 

From early times, an aspect which has re- 
ceived considerable attention, is possible non- 
linear modifications of quantum theory, and 
this is not necessarily because of the mea- 



surement problem. Most fundamental differ- 
ential equations which describe physical phe- 
nomena are non-linear, with linearity being 
a convenient approximation in some appro- 
priate limiting cases. Why then should an 
equation as fundamental as the Schrodinger 
equation be a singular exception to this rule? 
[It is of course known that there are very 
strong bounds on non-linearity in the atomic 
domain, see for instance the experiment de- 



tic non-linearity to explain measurement, the 



pioneering paper is due to Pearle (Pearle 



scribed in (Bollinger et al, 1989)]. Non- 



linear quantum theories may be classified as 
deterministic non-linear, and stochastic non- 
linear. For discussions on deterministic non- 
linear quantum mechanics the reader is re- 



ferred to the works by Weinberg (Wein 



berg 1989afb ), by Goldin (Doebner and 



Goldin,1992 Goldin, 2000), and (Bialynicki 



Birula and Mycielski, 1976). It has often 



been suggested, and demonstrated, though 
perhaps not universally so, that deterministic 
non-linear modifications result in superlumi- 



nal propagation (Ghirardi and Grassi, 1991 



Gisin 


1990; 


Polchinski 


1991 



with the fact that stochasticity appears to 
be an essential ingredient for explaining the 
origin of probabilities, has meant that inves- 
tigations of a universal dynamics have tended 
to focus on stochastic non-linearities; see for 



Diosi , 


1988alb; 


Gisin , 


1981. 


1984 



1989 Gisin and Rigo, 1995 Weinberg, 2011) 



With regard to the application of stochas- 



1976 ) - the paper is aptly titled "Reduction of 
the state vector by a non-linear Schrodinger 
equation". Pearle proposed to replace the 
Schrodinger equation by a non-linear one, 
during measurement, and that certain vari- 
ables which take random values just after 
the quantum system interacts with the ap- 
paratus, drive the system to one or the other 
outcomes, thus breaking superposition. For 
the choice of these random variables he sug- 
gested the phases of the state vectors imme- 
diately after the measurement. An appro- 
priate assignment of the probability distri- 
bution of these phases over the allowed pa- 
rameter space leads to the Born rule. It is 
noteworthy that this assignment of the prob- 
ability distribution is something which has 
to be put in by hand, keeping in mind what 
probability rule one wants to emerge. This is 
one aspect where phenomenology and under- 
lying theories need to do better even today: 
there should be a fundamental reason for the 
probability distribution over the stochastic 
variables, which inevitably implies the Born 
rule. Two important missing pieces, in order 
to consider the proposed dynamics a univer- 
sal dynamics for all physical systems, were 
the preferred basis on which the wave func- 
tion should collapse, as well as the trigger 
mechanism. Both limitations are overcome 
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by the GRW model. Further investigations 



by Pearle were reported in (Pearle, 1979 



1982, 1984, 1989b) and reviewed in (Pearle 



1999a). 



The next major advance came from Ghi- 



rardi, Rimini and Weber (Ghirardi et al. 



1986) in a seminal paper titled "Unified dy- 



namics for microscopic and macroscopic sys- 
tems" and the model has come to be known 
as the GRW model. There were two guiding 
principles for this dynamical reduction model 
(also known as QMSL: Quantum Mechanics 
with Spontaneous Localization): 

1 . The preferred basis - the basis on which 
reductions take place - must be chosen in such 
a way as to guarantee a definite position in 
space to macroscopic objects. 

2. The modified dynamics must have lit- 
tle impact on microscopic objects, but at 
the same time must reduce the superposi- 
tion of different macroscopic states of macro- 
systems. There must then be an amplifica- 
tion mechanism when moving from the micro 
to the macro level. 

The reduction is achieved by making the 
following set of assumptions: 

1. Each particle of a system of n distin- 
guishable particles experiences, with a mean 
rate Ag RW , a sudden spontaneous localization 
process. 

2. In the time interval between two 
successive spontaneous processes the system 



evolves according to the usual Schrodinger 
equation. 

In their model, GRW introduced two new 
fundamental constants of nature, assumed to 
have definite numerical values, so as to re- 
produce observed features of the microscopic 
and macroscopic world. The first constant, 
-^grw ~ 10 16 seconds, alluded to above, de- 
termines the rate of spontaneous localization 
(collapse) for a single particle. For a compos- 
ite object of n particles, the collapse rate is 



(A 



GRW 



n 



seconds. The second fundamental 



constant is a length scale Tc ~ 10 -5 cm which 
is related to the concept that a widely spaced 
wave-function collapses to a length scale of 
about rc during the localization. 

A gravity based implementation of the 



GRW model was studied by Diosi (Diosi 



1989) and generalized by (Ghirardi et al. 



1990a) 



The GRW model has been upgraded into 
what is known as the CSL (Continuous Spon- 
taneous Localization) model by Ghirardi, 



Pearle and Rimini (Ghirardi et al., 1990c). 



In CSL a randomly fluctuating classical field 
couples with the particle number density op- 
erator of a quantum system to produce col- 
lapse towards its spatially localized eigen- 
states. The collapse process is continuous in 
time, and this allows to express the dynam- 
ics in terms of a single stochastic differential 
equation, containing both the Schrodinger 
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evolution and the collapse of the wave func- ment, the emergence of the Born rule, the ab- 

tion. The narrowing of the wave function sence of macroscopic superpositions, and the 

amounts to an increase in the energy of the excellent matching of the linear theory with 

particle, and actually amounts to a tiny vio- experiments for microscopic systems, 
lation of energy conservation. Various studies and arguments suggest 

An outstanding open question with regard that the structure of this equation is very 

to the dynamical reduction models is the ori- rigid and tightly controlled, once one as- 

gin of the random noise, or the randomly flue- sumes [as is true here] that the evolution 

tuating classical scalar field, which induces is norm-preserving, and secondly, superlumi- 



collapse. nal propagation is not possible (Adler, 2004 



The current status of the Spontaneous Gisin 1989). There is then a unique relation 



Collapse models is discussed in detail in Sec- between the coefficient a/X of the diffusion 

tion II. [stochastic] term and the coefficient —A/2 of 

A modern approach to stochastic reduc- the drift term : Drift Coefficient = -2 (Dif- 

tion is to describe it using a stochastic non- fusion Coefficient) 2 . This is the well-known 

linear Schrodinger equation, an elegant sim- martingale structure for a stochastic differen- 

plified example of which is the following one- t ial equation. 

particle case [known as QMUPL : Quantum In QMUPL, stochastic fluctuations take 

Mechanics with Universal Position Localiza- place only in the time direction and hence 



tion (Diosi, 1989)] [See Section II for details] : there is only one free parameter, i.e. A. In 

contrast, in the CSL model the stochastic 



#(t) 



^Hdt + V\(q {q) t )dW t fluctuations exist over space too, and hence 



\(q-{q)t?dt 



there is a second free parameter r c (analo- 



gous to GRW) which defines the scale of spa- 

q is the position operator, (q) t is its expecta- tial localization. Of course, in the QMUPL 

tion value, and A is a constant, characteris- and CSL models the stochastic process acts 

tic of the model, which sets the strength of continuously, unlike in GRW, wherein the 

the collapse. W t is a Wiener process which stochastic jumps are discontinuous and dis- 

describes the impact of stochasticity on the crete. In fact the QMUPL model can be 

dynamics. As for the GRW and CSL mod- understood as a scaling limit of the GRW 



els, this equation can be used to explain the process (Diirr et al, 2011) [the collapse fre 



collapse of wave-function during a measure- quency goes to infinity and the spread r c goes 
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to zero in such a way that their product re- variables and their conjugate momenta, and 

mains a constant]. one should first know the Hamiltonian or 

Part of the experimental effort on test- the action [for a path-integral formulation] of 

ing quantum mechanics, discussed in detail the classical system. This is unsatisfactory, 

in Section IV, is devoted to testing the va- In the canonical formulation, one then pro- 



lidity of equations such as (14) above, and poses canonical commutation relations such 

measuring / setting bounds on the rate con- as [q, p] = ih in an ad hoc manner. Why 

stant A and the length scale r c . should these be the relations, unless one al- 

Underlying Theories ready knows that they lead to results which 

Phenomenological models of dynamical match with experiments? It would be desir- 

wave-function collapse propose an ad hoc able to derive quantum theory from a starting 

modification of quantum mechanics, albeit point which is not classical mechanics, and 

retaining certain features such as norm- then obtain classical mechanics as an approx- 

preservation and no superluminal propaga- imation [and explain quantum measurement 

tion. In principle, there should be strong in the process]. The theory of Trace Dynam- 

underlying theoretical reasons which make a developed by Adler and collaborators does 



compelling case for a modified quantum the- well in progressing towards this goal ( [Adler 



ory, rendering the phenomenological models |1994| |2004[ |2006| |Adler and Millard[ [1996J). 



inevitable. Here we mention three different Trace Dynamics [TD] assumes that the 
theoretical developments in this connection, underlying theory is a classical dynamics of 
two of which arise from attempts to remove Grassmannian matrices, living on a given 
one or the other fundamental incompleteness space-time. However, this classicality does 
in the formulation of quantum theory, and no t mean that TD is a 'hidden variables' the- 
the third investigates how gravity might play ory . f or the eventual description is at an av- 
an effective role in wave-vector reduction. eraged level, where no reference is made to 
Trace Dynamics the matrices which have been coarse-grained 
Classical mechanics is supposed to be a over. The matrices satisfy the standard La- 
limiting case of quantum theory. And yet, grangian and Hamiltonian Dynamics, but 
in its canonical formulation, quantum the- as a consequence of global unitary invari- 
ory assumes a prior knowledge of classical ance, the theory possesses a remarkable addi- 
dynamics! In order to 'quantize' a system, tional conserved charge, not present in point- 
one should know the classical configuration particle mechanics. This is the Adler-Millard 
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charge ( Adler and Millard] |1996[ ) ory is shown to emerge as the thermodynamic 



q _ p.] _ ^^{g- p } (15) a PP rox i m ation to an underlying classical dy- 

i j namics of Grassmann matrices possessing a 

where the first sum is over commutators of global unitary i nvar iance. 



bosonic matrices, and the second sum is 
over anti-commutators of fermionic matrices. 



Perhaps the greatest asset of TD is to be 

able to go beyond this stage and address the 

[See Section III for details! . This conserved , ^ 

L 1 quantum measurement problem m a natu- 

charge, which has the dimensions of action, , A , . 

° ' ' rai manner. Quantum theory emerges m the 

plays a central role in the emergence of quan- , , , . , . c 

^ J thermodynamic approximation of the statis- 

tum theory at a coarse-grained level. , , c , , , n . 

J ° tical mechanics of the underlying matrix me- 

Assuming that these matrix degrees of chanics. Next, it is pertinent to consider the 

freedom are at a level sufficiently 'micro- impact of Brownian motion fluctuations - re- 

scopic' [e.g. at the Planck scale] that we do markably these modify the Schrodinger equa- 

not observe them in our routine laboratory tion and provide the necessary stochastic ele- 

experiments, a statistical thermodynamics of nient for the collapse process to operate, and 

this matrix dynamics is constructed. An for the origin of probabilities. Subject to cer- 

equipartition theorem for the thermodynam- tain crucial assumptions for which one would 

ically averaged quantities is derived, which eventually like to find a theoretical basis, 

results in the Adler- Millard charge being uni- the modified Schrodinger equation is a non- 

formly distributed across the averaged com- linear and non-unitary [but norm-preserving] 

mutators, each of which is assumed to equal stochastic equation of the type used in the 

Planck's constant. This is the origin of the CSL model. In this way, Trace Dynam- 

quantum commutation relations. As a con- ics, through its thermodynamic limit and the 

sequence of the assumed invariance of ther- associated statistical fluctuations provides a 

modynamic averages under constant shifts theoretical underpinning for the phenomeno- 

in phase space, a Ward identity is derived, logical collapse models. 

which under suitable assumptions shows that TD is perhaps the most well-developed un- 

the thermally averaged q's and p's satisfy derlying theory one has at present for col- 

Heisenberg equations of motion. A relativis- lapse phenomenology. Hence in Section III 

tic quantum field theory is arrived at, and a we give a detailed presentation of the physics 

non-relativistic Schrodinger equation holds in and mathematics of TD, leading to wave- 

the finite particle limit. Thus quantum the- vector reduction, and we also point out the 
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open problems of TD which remain to be ad- theory which does not refer to classical time, 

dressed. Such a reformulation can be shown to be 

Quantum theory without classical spacetime the limiting case of a non-linear theory, with 

Quantum theory requires an external clas- the non-linearity becoming important at the 

sical time for describing evolution. This is of Planck mass scale - The non-linearity is pos- 

course so obvious and essential that it is al- sibl y stochastic, and could have implications 

most never stated explicitly! However this for resolution of the quantum measurement 

dependence on an external classical time is Problem. Tentative heuristic discussions to- 

perhaps the greatest incompleteness of quan- wards th is investigation have been given in 



turn theory. Such a time is part of a classi- ( Singh 2006 2009 ) 



cal spacetime geometry which is produced by A detailed systematic program to develop 

classical matter fields according to the laws of a formulation of quantum theory without 

general relativity. But classical matter fields classical time and to study its impact on 

are a limiting case of quantum fields. If there quantum measurement has recently been be- 



were no classical fields in the Universe, but gun, and is qualitatively described in (Singh 



only fields subject to quantum fluctuations, 2011). The key symmetry principle here is 



there will be no definite metric available to that basic laws should be invariant under co- 
describe the spacetime geometry. An argu- ordinate transformations of non-commuting 
ment due to Einstein, known as the Einstein coordinates. The motivation being that 



hole argument (Christian, 1998) then implies if quantum fluctuations destroy a classi- 



that if the metric is subject to quantum flue- cal spacetime manifold, a possible replace- 
tuations, there is no longer available an un- ment for ordinary spacetime could be a non- 
derlying classical spacetime manifold. It is commutative spacetime. This approach pro- 
then not possible to describe quantum evolu- poses to generalize Trace Dynamics by raising 
tion. time, and space, to the level of matrices [op- 



We see once again that via its dependence erators]. This has been done in (Lochan and 



on external time, quantum theory depends Singh| 2011) and it has been shown that by 



on its classical limit [the required presence defining a non- commutative space-time met- 

of a Universe dominated by classical mat- ric as a Trace over the space-time operators, 

ter] . This is unsatisfactory from a fundamen- a Poincare invariant dynamics can be con- 

tal point of view, and hence there must ex- structed. We call this a generalized Trace Dy- 

ist an equivalent reformulation of quantum namics. Evolution is described with respect 
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to the scalar constructed by taking Trace over The program described here aims to ad- 
the non-commutative metric - this is the ana- dress a limitation of Trace Dynamics - a ma- 
log of the ordinary proper time. trix treatment for matter fields while leav- 
The next step is to construct, a la TD, in S spacetime as having point structure, thus 
a statistical mechanics for this generalized leavin § spacetime untouched. We regard 
matrix dynamics and obtain the equilibrium such a limitation as one which should be ad- 
thermodynamic approximation - this yields dressed " in the P rocess we see that remov - 
a generalized quantum theory which has an in S time from Quantum theory drives us to a 
energy-time commutation relation and a gen- starting P omt [generalized Trace Dynamics] 
eralized Schrddinger equation with an opera- whose eventual outcome is a possible resolu- 
tor time as one of the configuration variables. tion for the measurement problem. We thus 
This is the sought for reformulation of quan- wish to assert that there is a dee P connec - 
tum theory which does not refer to an ex- tion between the problem of time in quan- 



ternal classical time flLochan et al.\ 120121). tum theor y and the measurement problem in 



If the Universe is dominated by macroscopic quantum theory ( jSingh[ |20T2|) . Addressing 



objects, the consideration of Brownian mo- the former wiU Possibly compel us to con- 

tion fluctuations should yield position local- sider a modification of quantum theory and 

ization and the concurrent emergence of a that modification will have a bearing on the 

classical space-time. This is the classical Uni- measurement problem. 

verse, dominated by classical macroscopic ob- Since this work is at present in an early 

jects and in possession of a classical space- stage of development, we will not discuss it 

time. This Universe has a 'sprinkling' of an y furt her in the remainder of the review, 
quantum fields and non-relativistic quantum Gravity induced wave-function collapse 

systems. On the backdrop of this classical The fact that the fundamental mass scale, 

Universe one can postulate standard quan- Planck mass Mpi = (hc/G) 1 ^ 2 ~ 1CT 5 grams 

tum theory [now that an external time is is not far from the scale where the micro 

given] and then proceed to implement the to macro transition takes place has often in- 

program of Trace Dynamics to derive quan- trigued some physicists. Mass seems to have 

tum dynamics from matrix mechanics, for something to do with deciding which objects 

this 'sprinkling' of quantum matter fields on are quantum, and which are classical, and 

the classical spacetime background, and to mass also produces gravity. Could gravity 

resolve the attendant measurement problem, thus play some role in causing wave-function 
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collapse, and in localization of macro- ob- 
jects? The idea that gravity might somehow 
be responsible for wave-function collapse has 
been seriously pursued by Karolyhazy and 



unified mathematical description has taken 
shape only over the last three decades or so. 
And yet, none other than Einstein himself 



(Schilpp, 19490 saw it this way early on, and 



collaborators (Karolyhazy, 1966; Karolyhazy 



et al, 1986), by (Diosi, 1987), and by Pen 



rose (Penrose, 1996). Penrose's proposal is 



also the subject of an important ongoing ex- 



periment aimed at testing it (Marshall et al. 



2003). These issues will be discussed in Sec- 



tion III. 

Experimental Tests 

The Copenhagen interpretation was a 
need of the times when it was proposed: pio- 
neering experiments were being carried out 
for atomic systems. The measuring appa- 
ratus was a classical object, and the Born 
probability rule had to be invoked to ex- 
plain the random outcomes of measurements. 
For some, this dual aspect of quantum the- 
ory - unitary evolution followed by wave- 
packet reduction - was the 'truth' in quan- 
tum theory; this is how nature is. For others, 
this was completely unacceptable, and rein- 
terpretations and new mathematical formu- 
lations such as many-worlds, Bohmian me- 
chanics, and decoherent histories, were devel- 
oped. However, the idea that quantum the- 
ory may be an approximation to a holistic 
theory which better explains both the uni- 
tary and reductionist aspects as limits of a 



had this to say about quantum theory: 

. . it would, within the framework of fu- 
ture physics, take an approximately analo- 
gous position to that of statistical mechanics 
within the framework of classical mechanics". 

In the light of the theory of Trace Dynam- 
ics and models of spontaneous wave function 
collapse, these words are prophetic. These 
modern ideas suggest the emergence of prob- 
abilities as a consequence of thermodynamic 
averaging in a deterministic theory, and the 
related significance of stochastic fluctuations. 
Above all, their predictions for results of ex- 
periments differ from the predictions of quan- 
tum theory. The difference will be far too 
small to be detectable for an atomic system, 
but starts becoming significant as the size of 
the system is increased. The best example of 
an experiment which could detect such a dif- 
ference is double slit interference. If an object 
of mass m is directed at a suitably prepared 
double slit, with appropriate slit width and 
separation between the slits, quantum theory 
predicts that an interference pattern will be 
seen on the screen no matter what the value 
of m. Not so, say collapse models. Accord- 
ing to these models, the superposition state 
created after the object has passed the slits 
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lasts only for a finite time r, where r de- namics took over from classical dynamics in 

creases with increasing m, and its value can the atomic domain, and (iii) Newton's inverse 

be calculated precisely from a given theoret- square law of gravitation was replaced by the 

ical model. Thus, according to these mod- laws of general relativity for strong gravita- 

els, if the time of travel from the slits to the tional fields. 

screen is greater than r, superposition will interference experiments with matter have 
break down before the screen is reached, and & fascinating history5 with a quest develop- 
no interference pattern will be seen. This is ing over decadeg to tegt SU p erposit i on using 
perhaps the cleanest confrontation that spon- larger and larger objectg _ Qver dghty yearg 
taneous collapse and gravity collapse models nave passed s i nce tne c i ass i c experiment by 
make with experiment. A successful diffrac- Davisson and Qermer in 1927 where interfer- 
tion experiment in the right mass domain will ence wag demoris trated with electrons. Land- 
irrefutably confirm or rule out these models. markg Qn the way included connrmat ion of 

One should of course stay cautioned interference for Helium [1930] and neutrons 

against assuming that quantum theory will [1988]. A great modern breakthrough came 

be successful through and through, and that in 1999 when interference was demonstrated 

interference will be seen for all values of m. for Ceo in the famous fullerene experiment 



The fact that a theory is extremely success- (Arndt et al. 1999). [There seems to be 



ful in one part of the parameter space should prevalent a popular belief in some quarters 

not be taken as guarantee that it will con- that the discovery of quantum superposition 

tinue to be successful in a different part of in a molecule as 'large' as fullerene means the 

the parameter space - in the present instance end of all theories which predict breakdown 

the absence of macroscopic superpositions al- of superposition for large systems. This of 

ready provides reason for caution. And there course is not true - breakdown of superposi- 

are historical examples of long-standing sue- tion and the quantum-to-classical transition 

cessful theories eventually turning into ap- is expected around 10 6 amu to 10 9 amu.] This 

proximations to more general theories when opened the door for larger molecules, and 

their extrapolation into a new part of the pa- today, a decade later, interference has been 

rameter space failed to be confirmed by ex- demonstrated for molecules with 7,000 nucle- 



periment: (i) classical mechanics became an ons (Gerlich et al. 2011). Proposed future 



approximation to special relativity at speeds interferometry experiments plan to push the 
close to the speed of light, (ii) quantum dy- limit to macromolecules with a million nucle- 
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ons and beyond, going up to molecules with a Apart from direct laboratory experiments, 
100 million nucleons. Doing so involves over- collapse model parameters are also con- 



coming great technological challenges ( Horn- strained by their effect on known measure- 



berger et al. t 2011), and there are many or- ments. Section IV discusses the various 

ders of magnitudes in the mass scale yet to experimental tests of the phenomenological 

be covered. But we certainly live in excit- models, 
ing times where predictions of collapse mod- 
els and gravity based models are being tested 

C. Plan and outline of the article 

by these experiments, and constraints are be- 
ing put on model parameters. See Section IV Sections II, III and IV are the main parts 
for details. Q f ^he review. Sec. II reviews phenomeno- 
Also, this is perhaps a good place to clear logical models of spontaneous wave function 
another misconception regarding the domain collapse, which explain the absence of macro- 
over which quantum mechanics has been scopic superpositions, via a stochastic non- 
tested. Various remarkable macroscopic in- linear modification of the Schrodinger equa- 
ternal states have been achieved experimen- tion. Sec. Ill gives a review of Trace Dy- 
tally, in which an enormous collection of in- namics and gravity-induced-collapse as pos- 
ternal degrees of freedom behave as a collec- sible underlying theories for the phenomenol- 
tive one-particle coherent state. We have in ogy discussed in Sec. II. Sec. IV is a re- 
mind of course systems such as superconduc- view of the techniques and results of ongoing 
tors, superfluids and Bose-Einstein conden- and planned experiments which are testing 
sates. The existence of such states however the proposed phenomenological models. Sec. 
does not explain why macroscopic objects are V provides a critique of the current under- 
not found in superposition of position states, standing on the theoretical and experimental 
Quantum mechanics may yet have to be mod- front, and lists open problems, 
ified so that the modified theory can explain Sec. II begins by introducing spontaneous 
the absence of position superpositions, but collapse, and recalls the various collapse mod- 
the modified theory will certainly continue to els that have been proposed. The original 
successfully explain a collective phenomenon GRW model is then introduced. This is then 
such as superconductivity. In other words, followed by a detailed review of the QMUPL 
the discovery of superconductivity does not model, which is applied to show how stochas- 
solve / trivialize the Schrodinger cat paradox! ticity induces collapse, and how the Born 
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probability rule is derived. The possible ori- work of Penrose. 

gin of the noise field is discussed. The next Sec. IV on Experimental Tests starts 

sub-section discusses the most widely used by discussing the basics of the collapse the- 

(but physically equivalent to GRW) collapse ory necessary for performing and interpret- 

model, i.e. the CSL model, and its general- ing the diffraction experiments with macro- 

izations. Lastly, the current understanding of molecules. Matter- wave interferometry and 

the numerical values of the two parameters of optomechanics experiments with mechani- 

the collapse model is reviewed. C al cantilevers are reviewed in detail. Cav- 

Sec. Ill reviews Adler's Trace Dynamics ity optomechanics with micro-spheres and 

as a candidate fundamental theory for Spon- nanoparticles is discussed, followed by a re- 

taneous Localization. The fundamental ma- view of new developments which combine 

trix degrees of freedom of the theory are in- optical tweezing techniques with near-field 

traduced, and their dynamics described. The matter-wave interferometry. The challenges 

conserved charges of the theory, including the proposed to these experiments by various 

all-important Adler-Millard charge, are de- kinds of decoherence are considered. The 

rived. This is followed by the construction current bounds on collapse model parameters 

of the statistical mechanics and the canoni- coming from the diffraction experiments and 

cal ensemble for thermodynamic equilibrium from other measurement processes are sum- 

for the theory. Following this, an impor- marized. 
tant Ward identity, which is an analog of the 
equipartition theorem, is proved. It is shown 
how the commutation relations for quantum 
theory, and the Schrodinger equation, emerge 

at this coarse-grained level, from the micro- A T . . . „ , „ ,,„ 

to ' A. Introducing spontaneous collapses 

scopic theory. Finally, consideration of fluc- 
tuations described by Brownian motion leads Quantum Mechanics, in its standard text- 

to generalization from the Schrodinger equa- book formulation, refers only to the outcomes 

tion to the stochastic non-linear Schrodinger of measurements, but it has nothing to say 

equation, which makes contact with the CSL about the world as it is, independently of 

model. Subsequent sections describe the any measurement or act of observation. This 

gravity based models for collapse, based on is a source of serious difficulties, which have 



II. SPONTANEOUS COLLAPSE MOD- 
ELS 



the work of Karolyhazy et al., Diosi, and the been clearly elucidated e.g. by J. Bell (Bell 
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1990): It would seem that the theory is ex- to apply to anything but highly idealized In- 
clusively concerned about 'results of measure- oratory operations, are we not obliged to ad- 
ments', and has nothing to say about any- mit that more or less 'measurement-like' pro- 
thing else. What exactly qualifies some phys- cesses are going on more or less all the time, 
ical systems to play the role of 'measurer' ? more or less everywhere? Do we not have 
Was the wave function of the world waiting jumping then all the time? 
to jump for thousands of millions of years un- The basic idea behind the dynamical re- 
til a single- celled living creature appeared? Or duction program is precisely this: sponta- 
did it have to wait a little bit longer, for some neous and random collapses of the wave func- 
better qualified system ... with a Ph.D.? tion occur all the time, for all particles, 
Measuring devices, like photographic whether isolated or interacting, whether they 
plates and bubble chambers, are very sophis- form just an atom or a complex measuring de- 
feated and highly structured physical sys- vice. Of course, such collapses must be rare 
terns, which anyhow are made of atoms; we and mild for microscopic systems, in order 
then expect them to be ultimately described not to alter their quantum behavior as pre- 
in quantum mechanical terms by means of dieted by the Schrodinger equation. At the 
the Schrodinger equation. What else should same time, their effect must add up in such 
we expect, taking into account that physi- a way that, when thousands of millions of 
cists are trying to describe even the entire particles are glued together to form a macro- 
universe quantum mechanically? But if we scopic system, a single collapse occurring to 
describe measurements in this way, then the one of the particles affects the global sys- 
theory does not predict any definite outcome, tern. We then have thousands of millions of 
at the end of the process. The Schrodinger such collapses acting very frequently on the 
equation is linear, the superposition principle macro-system, which together force its wave 
holds, and it does so in such a way that all function to be very rapidly well-localized in 
possible outcomes are there simultaneously space. 

in the wave function, but none of them is se- On the mathematical level, the program is 

lected as the one which actually occurs. Yet, accomplished by modifying the Schrodinger 

if we perform a measurement, we always get evolution, introducing new terms having the 

a definite outcome. So we have a problem following properties: 

with Quantum Mechanics. • They must be non-linear: The new dy- 

Continuing quoting Bell: If the theory is namics must break the superposition 
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principle at the macroscopic level and 
guarantee the localization of the wave 
function of macro-objects. 

• They must be stochastic: When de- 
scribing measurement-like situations, 
the dynamics must explain why the 
outcomes occur randomly; more than 
this, it must explain why they are dis- 
tributed according to the Born proba- 
bility rule. On top of this, stochasticity 
is necessary because otherwise the non- 
linear terms would allow for faster than 
light communication. 

• There must be an amplification mecha- 
nism according to which the new terms 
have negligible effects on the dynamics 
of microscopic systems but, at the same 
time, their effect becomes very strong 
for large many-particle systems such as 
macroscopic objects, in order to recover 
their classical-like behavior. 

• They must not allow for superluminal 
signaling, as one would like to preserve 
the causal structure of spacetime. 

Looking carefully at these requirements, 
one soon realizes that they are very demand- 
ing: there is no reason beforehand, that they 
can be consistently fulfilled. One of the great- 
est merits of collapse models is to have shown 
that this program can be implemented in a 



consistent and satisfactory way. 

B. The zoo of collapse models 

In the literature, different collapse models 
have been proposed. A first characterization 
depends on the choice of the collapse oper- 
ators, i.e. on the basis on which the wave 
function is localized. Some models induce 



the collapse in the energy basis (Adler, 2002 



20041 lAdler et al\ [20011 lAdler and Brun 



2001 Adler and Horwitz 2000 ), (Brody and 



Hughston 




2002 


Hughston 




1996; Milburn 



1991), others in the momentum basis (Be 



natti et al. 1988), or the spin basis (Bassi 



and Ippoliti 2004; Pearle 2012). However, 



only models which collapse in the position 
basis make sure that different macroscopic 
superpositions rapidly collapse towards lo- 
calized states. To understand this, one can 
think of a superposition of two spatially sep- 
arated states of a macroscopic object, which 
have the same (or very similar) energy. In 
this case, an energy-based collapse model 
would not be able to collapse the superpo- 
sition fast enough, because the superposition 
in energy is null or negligible. Such a model 
would not be able to guarantee that macro- 
objects always occupy a definite position in 
space. Only space collapse models make sure 
that macroscopic objects always behave clas- 
sically, and therefore we will consider only 
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them in the following. 

Space collapse models can be conveniently 
grouped depending on the properties of the 
noise, which is responsible for the collapse. 
A first distinction is between white and non- 
white models. In white-noise models, the 
collapse-noise is assumed to be a Wiener pro- 
cess in time, and the resulting evolution is 
Markovian. All frequencies of the noise con- 
tribute to the collapse with the same weight. 
Examples of collapse models of this type are 



the GRW model (Ghirardi et al, 1986), the 



CSL model (Pearle, 1989a), (Ghirardi et al. 



1990c), the QMUPL model (Bassi, 2005 



Diosi 1989). In non- white noise models, the 
collapse-noise is taken to be a generic Gaus- 
sian noise, with mean equal to zero, and a 
generic correlation function. The correspond- 
ing dynamics turns out to be non-Mar kovian, 
and they are more difficult to analyze. A 
model of this kind is the non-Markovian 
QMUPL model ( [Bassi and Ferialdi[ |2009a|b| ) , 
while the non-Markovian CSL model is still 



ture. The wave function collapses, but at the 
same time the energy of the quantum sys- 
tem increases steadily; no dissipative effects 
are taken into account. This is a well-known 
feature of collapse models. Mathematically, 
these models are characterized by the fact 
that the wave function and the collapse-noise 
are coupled through the position operator 
only. The GRW model, the (Markovian and 
non-Markovian) CSL model, the (Markovian 
and non-Markovian) QMUPL model all be- 
long to this group. In the finite temperature 
models instead, the collapse-noise behaves 
like a reservoir at finite temperature. The 
wave function still collapses, but now dissi- 
pative terms are included (through a position 
and momentum coupling between the wave 
function and the noise), which thermalize any 
quantum systems to the temperature of the 
noise. The only such model so far available 



is the non-dissipative QMUPL model (Bassi 



under development (Adler and Bassi, 2007 



2008). General non-Markovian collapse mod- 



et al., 2005b), though also the other models 
can be generalized in this sense. Recently, the 
QMUPL model has been generalized in order 
to include both non-Markovian and dissipa- 



els have been discussed in (|Bassi and GuT] tive effects ( |Ferialdi and Bassi} \20U) . 



2002 


Diosi et al. 


1998 


Pearle , 


1993 



1996]. 

A second distinction is between infinite 
temperature and finite temperature models. 
In the first type of models, the collapse- 
noise acts like a reservoir at infinite tempera- 



A final distinction is between first quan- 
tized models and second quantized models. 
Models of the first type consider only system 
of distinguishable particles; the GRW model, 
the QMUPL model and its non-Markovian 
and/or dissipative generalization belong to 
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this group. Models of the second type are for- 
mulated in the language of quantum field the- 
ory and include systems of identical particles. 



The Tumulka-GRW model (Tumulka, 2006a) 



and the CSL model belong to this group. 

We also mention the earlier contributions 
of L. Diosi ( p6silp88b|c| ), N. Gisin ( pirn" 



1984, 1989) and J. Percival (Percival, 1999 



Schack et ai, 1995) to developing stochasti- 



cally modified Schrodinger equations for de- 
scribing the process of wave function collapse. 

Some comments are in order. The first 
one is that all space-collapse model are qual- 
itatively equivalent: they all induce the col- 
lapse of the wave function in space, and 
the collapse is faster, the larger the sys- 
tem. Of course, they can differ also in a 
significant way in the technical details, as 
we will see. The second comment refers to 
the nature of the stochastic character of the 
collapse process. One way to look at it— 
which corresponds to the original attitude to- 
wards these models — is that nature is intrin- 
sically stochastic, therefore stochastic differ- 
ential equations are the natural type of equa- 
tions for describing the dynamics of physi- 
cal systems. A new way to look at it is to 
assume that there is a random field, filling 
space, which couples to quantum matter in 
a non-standard way, and is responsible for 
the collapse of the wave function. The new 
terms in the modified Schrodinger equation 



are meant to describe such a coupling. Since 
this noise fills the whole space, most likely 
it has a cosmological origin. According to 
this scenario, the physically most reasonable 
collapse model is a model, where the col- 
lapsing field is "cosmologically reasonable", 
e.g. it has a typical cosmological correlation 
function and a typical cosmological tempera- 
ture. This could be the case for the colored- 
noise and dissipative CSL model, which how- 
ever has not been formulated yet. What one 
can do, is to extrapolate predictions from the 
other models already available. Therefore, in 
the following we will focus our attention on 
two of the above mentioned models: the CSL 
model, the one that more closely resembles 
the physically most reasonable model; and 
the QMUPL model, which is less physical, 
but has already been generalized in order to 
include dissipation, as well as colored noises, 
and is relatively easy to analyze mathemati- 
cally. 

A third comment is about the origin of the 
noise field. The important thing to bear in 
mind, is that this field cannot be a standard 
quantum field, otherwise we would fall back 
in the realm of standard quantum mechan- 
ics, with the superposition principle and the 
measurement problem. This field couples to 
quantum matter through an anti-hermitian 
and non-linear coupling. The most intrigu- 
ing guess is that this noise has a gravitational 
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origin. In fact, a gravitational background is 
part of the standard cosmological scenario; 
gravity is non-linear; gravity has not been 
successfully quantized yet, and we do not 
know today what shape a quantum theory of 
gravity will eventually take. Gravity-induced 
collapse models have been formulated in the 
literature, and we will discuss them in Sec. 
III.B. 

A fourth comment is about the relativistic 
extension of collapse models. All models pre- 
viously listed are non-relativistic. Their gen- 
eralization to relativistic quantum field the- 
ories has not been successful. The reason is 
very simple to understand: the collapse of the 
wave function is an instantaneous process, or 



at least faster than light (Maudlin, 2011) 



This is a necessary requirement, in order to 
reproduce non-local quantum correlations en- 



coded in Bell inequalities (Bell 1987b), which 



have been verified experimentally. An instan- 
taneous collapse process is not welcome in a 
relativistic framework, hence the difficulty in 
formulating relativistic collapse models. We 
will discuss this issue in Sec. II. I 

C. The GRW model 

In order to appreciate how collapse mod- 
els work, and what they are able to achieve, 
we briefly review the GRW model, the first 
consistent model proposed. Though it is not 



expressed in terms of a compact stochastic 
differential equation, it has the advantage of 
being physically very intuitive. In presenting 
the model, we follow the exposition of J. Bell 
pen] |1987a[ ) [reprinted in pen] |1987b[ )] in 
terms of discrete jumps of the wave function. 

Let us consider a system of iV particles 
which, only for the sake of simplicity, we 
take to be scalar and spinless; the GRW 
model is defined by the following postulates: 

States. The state of the system is repre- 
sented by a wave function ■0(xi,x 2 , . . • Xjv) 
belonging to the Hilbert space C 2 (R 3N ). 
Spin and other internal degrees of freedom 
are ignored for the sake of simplicity. 

Dynamics. At random times, the wave 
function experiences a sudden jump of the 
form: 



V>i(xi,x 2 , . . .Xjv) — ► 

L n (x)^(xi,x 2 , . . .Xjv) 



-, (16) 



||L„(x)^(xi,x 2 , . . .Xjv) 

where ipti^i, x 2, • • • x tv) is the state vector of 
the whole system at time t, immediately prior 
to the jump process. L n (x) is a linear opera- 
tor which is conventionally chosen equal to: 
1 



L n teq 



-(q„-x) 2 /2r2. 



where tq is a new parameter of the model 
which sets the width of the localization pro- 
cess, and q n is the position operator associ- 
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ated to the n-th particle of the system; the 
random variable x corresponds to the place 
where the jump occurs. Between two consec- 
utive jumps, the state vector evolves accord- 
ing to the standard Schrodinger equation. 

The probability density for a jump taking 
place at the position x for the n-th particle 
is given by: 

p n (x) = ||L„(x)^ t (xi,X2, • • .Xjv)!! 2 , 

(18) 

and the probability densities for the different 
particles are independent. 

Finally, it is assumed that the jumps are 
distributed in time like a Poissonian process 
with frequency A GRW ; this is the second new 
parameter of the model. 

The standard numerical values for tq and 
Agrw are: 



A 



GRW 



HT 16 sec" 1 , 



rc — 10 rn. 



(19) 



We will come back to the issue of numerical 
value of these parameters in Sec. II. J. 

Ontology. In order to connect the math- 
ematical formalism with the physical world, 
one needs to provide an ontology, which is 
rather straightforward for collapse models. 
Let m n be the mass associated to the n-th 
"particle" of the system (one should say: to 
what is called "a particle", according to the 



standard terminology); then the function: 

...d 3 x A r|^(x 1 ,x 2 ,...x 7V )| 2 (20) 



represents the density of mass (Ghirardi 



et al, 1995) of that "particle" in space, at 



time t. 

These are the axioms of the GRW model: 
as we see, words such as 'measurement', 'ob- 
servation', 'macroscopic', 'environment' do 
not appear. There is only a universal dynam- 
ics governing all physical processes, and an 
ontology which tells how the physical world 
is, according to the model, independently of 
any act of observation. 

The GRW model, as well as the other 
dynamical reduction models which have ap- 
peared in the literature, has been exten- 



sively studied (see (Bassi and Ghirardi, 2003) 



and QPearle[ |1999aQ for a review on this 
topic); in particular — with the numerical 
choice for A GHW and rc given in (19) — the fol- 



lowing three important properties have been 
proved, which we will state in more quanti- 
tative terms in the following section: 

• At the microscopic level, quantum sys- 
tems behave almost exactly as pre- 
dicted by standard Quantum Mechan- 
ics, the differences being so tiny that 
they can hardly be detected with 
present-day technology. 
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• At the macroscopic level, wave func- 
tions of macro-objects are almost al- 
ways very well localized in space, so 
well localized that their centers of 
mass behave, for all practical purposes, 
like point-particles moving according to 
Newton's laws. 

• In a measurement-like situation, e.g. 
of the von Neumann type, the GRW 
model reproduces — as a consequence of 
the modified dynamics — both the Born 
probability rule and the standard pos- 
tulate of wave-packet reduction. 

In this way, models of spontaneous wave func- 
tion collapse provide a unified description of 
all physical phenomena, at least at the non- 
relativistic level, and a consistent solution to 
the measurement problem of Quantum Me- 
chanics. 

It may be helpful to stress some points 
about the world-view provided by the GRW 
model, and collapse models in general. Ac- 
cording to the ontology given by the third ax- 
iom, there are no particles at all in the theory! 
There are only distributions of masses which, 
at the microscopic level, are in general quite 
spread out in space. An electron, for exam- 
ple, is not a point following a trajectory — as 
it would be in Bohmian Mechanics — but a 
wavy object diffusing in space. When, in a 
double-slit experiment, it is sent through the 
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two apertures, it literally goes through both 
of them, as a classical wave would do. The 
peculiarity of the electron, which qualifies it 
as a quantum system, is that when we try to 
locate it in space, by making it interact with 
a measuring device, e.g. a photographic film, 
then, according to the collapse dynamics, its 
wave function very rapidly shrinks in space 
till it gets localized to a spot, the spot where 
the film is exposed and which represents the 
outcome of the position measurement. Such 
a behavior, which is added ad hoc in the stan- 
dard formulation of Quantum Mechanics, is a 
direct consequence of the universal dynamics 
of the GRW model. 

Also macroscopic objects are waves; their 
centers of mass are not mathematical points, 
rather they are represented by some func- 
tion defined throughout space. But macro- 
objects have a nice property: according to the 
GRW dynamics, each of them is always al- 
most perfectly located in space, which means 
that the wave functions associated with their 
centers of mass are appreciably different from 
zero only within a very tiny region of space 
(whose linear extension is of order 10~ 14 m 
or smaller, as we shall see), so tiny that they 
can be considered point-like for all practical 
purposes. This is the reason why Newton's 
mechanics of point particles is such a satis- 
factory theory for macroscopic classical sys- 
tems. 



Even though the GRW model contains ator Ti[-], i = 1,2, ...N appears for each 

no particles at all, we will keep referring to particle. For ordinary matter — and with a 

micro-systems as 'particles', just as a matter good approximation — one can separate the 

of convenience. center-of-mass motion from the internal mo- 



Though the collapse dynamics is expressed tion ( Ghirardi et al. 1986). The reduced den- 



entirely in terms of the wave function, not of sity matrix for the internal motion obeys the 
the density matrix, in order to eliminate any standard Schrodinger equation, while that for 



possible ambiguity about the nature of the the center of mass is equivalent to Eqn. (21) 



collapse, it is nevertheless convenient to look where now the collapse rate entering the def- 

at the collapse dynamics for the density ma- inition of the operator T[-} is iVA GHW , with N 

trix, to analyze specific features of the model, the total number of particles making up the 

The 1-particle master equation of the GRW object. This is a manifestation of the am- 



model takes the form (Ghirardi et al] 1986): plification mechanics, perhaps the most im- 



d , . ir TT r , , portant feature of collapse models: the wave 

p{t) = --[H,p{t)\-T[p{t)\, (21) 



dt ry ' h 
where H is the standard quantum Hamilto- 
nian of the particle, and T[-] represents the 
effect of the spontaneous collapses on the par- 
ticle's wave function. In the position repre- 
sentation, this operator becomes: 



function of an object collapses with a rate 
which is proportional to the size of the sys- 
tem. This is the mathematical reason why 
collapse models can accommodate both the 
quantum dynamics of microscopic systems 
(negligible collapse rate) and the classical dy- 
x|T[p(£)]|y) — (22) naimcs °f macroscopic systems (fast collapse) 



A G Rw[l-e-( x -^ 2 / 4 '-e](x|p(t)|y). 



within one unified dynamical principle. 



D. The QMUPL model 



As expected, the effect of the spontaneous 
collapse is to suppress the off-diagonal el- 
ements of the density matrix, with a rate We now focus our attention on the 



proportional to A GRW , depending also on the QMUPL (Quantum Mechanics with Univer- 

distance between the off-diagonal elements: sal Position Localizations) model. As pre- 

distant superpositions are suppressed faster viously anticipated, the reason is that this 

than closer ones. model has the virtue of being both physi- 

The many particle master equation is the cally realistic, though very simplified corn- 



generalization of Eqn. (21), where an oper- pared to the more realistic GRW and CSL 
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models, and mathematically simple enough 
to be analyzed in great detail. The axioms 
defining this model are the same as those of 
the GRW model, with the only difference that 
the dynamics is described by a stochastic dif- 
ferential equation. The one-particle equation 
takes the form (for simplicity, we will work 
only in one dimension in space): 



dipt 



h 



Hdt + V\(q - (q)t)dW t 



A 



(q - (q) t fdt 



(23) 



where q is the position operator of the parti- 
cle, {q) t = (ipt\q\ipt) is the quantum expec- 
tation, and Wt is a standard Wiener pro- 
cess. For simplicity, we work in only one spa- 
tial dimension, the generalization to three di- 
mensions being straightforward. The collapse 
constant A sets the strength of the collapse 
mechanics, and it is chosen proportional to 
the mass m of the particle according to the 
formulcEJ 



A 



711 



m 



Ao) 



(24) 



where m is the nucleon's mass and Aq mea- 



sures the collapse strength (Bassi, 2005). If 



we set Aq — 10 2 m 2 s , then the strength 



2 One should keep in mind that the collapse strength 
depends on the type of model. For the GRW 
model, A G rw is a rate. For the QMUPL model, 
A has the dimensions of an inverse time, times an 
inverse square length. The two constants are re- 
lated, by requiring that the collapse strengths ac- 
cording to different models coincide in the appro- 



priate limit ( Bassi and Diirr 2009 1 



of the collapse mechanism according to the 
QMUPL model corresponds to that of the 
GRW and CSL models in the appropriate 



limit (Bassi and Diirr 2009). Note also 



that the QMUPL model is defined in terms 
on only one parameter (A), while the GRW 
model (and similarly the CSL model) is de- 
fined in terms of two parameters (A GRW and 
rc)- 

We will come back to the numerical values 
of the collapse parameter in Sec. II. J. The 
generalization for a many-particle system can 
be easily obtained by considering the position 
operator qi of every particle, each coupled to 
a different Wiener process W t (i) . The struc- 
ture remains the same, with a sum to include 
the contribution to the collapse coming from 
each particle. 



As we expect, Eqn. (23) contains both 



non-linear and stochastic terms, which are 
necessary to induce the collapse of the wave 
function. In order to see this, let us consider 
a free particle (H = p 2 /2m), and a Gaussian 
state: 

ipt{x) = exp [-a t (x -x t ) 2 + ik t x + j t ] ■ 

(25) 

It is not too difficult to show that ipt{x) 



is solution of Eqn. (23), provided that the 



time dependent functions in the exponent 
solve appropriate stochastic differential equa- 



tions (Bassi, 2005). In particular, the equa- 
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tions for a t which controls the spread both in will get by using the Schrodinger equation in- 
position and momentum, for the mean posi- stead of Eqn. ( 23 ) . 
tion x t and the mean momentum k t are^] 



da* 



A (a t ) 

m 



dt, 



dx t = — k t dt H dW t , 

m 2af 



dk t 



af 



(26) 

(27) 
(28) 



Eqn. (26) is deterministic and easy to solve. 



The spreads in position and momentum: 




(29) 



are given by the following analytical expres- 
sions: 



O,, 



H cosh(a;t + (pi) + cos(ut + ip 2 ) 



mu sinh(u;t + (fi) + sin(u;t + (f 2 ) 



(30) 



O n 



' hmuj cosh(wt + (pi) — cos(ut + (p 2 ) 



2 sinh(wt + ifi) + sin(wt + ip 2 ) 



with: 



u = 2\l— ~ lO^s" 1 . 
m 



(31) 



(32) 



The two parameters <p\ and i p> 2 are functions 
of the initial condition. Also, note that set- 
ting Aq = will give the same results as one 



3 The superscripts "R" and "I" denote, respectively, 
the real and imaginary parts of the corresponding 
quantities. 



Eqns. (30) and (31) tell that the spreads 



in position and momentum do not increase 
in time, but reach an asymptotic final value 
given by: 



0-9(00 ) = 

and: 
Op(oo) = 

such that: 



h 




10 



-19 




m, (33) 



Kg m 



sec 



o q {oo)o p {oo) 



h 

71 



(34) 



(35) 



which corresponds to almost the minimum al- 
lowed by Heisenberg's uncertainty relations. 
As we see, the spread in position does not in- 
crease indefinitely, but stabilizes to a finite 
value, which is a compromise between the 
Schrodinger's dynamics, which spreads the 
wave function out in space, and the collapse 
dynamics, which shrinks it in space. For mi- 
croscopic systems, this value is still relatively 
large (a q (oo) ~ lm, for an electron, and 
~ 1mm, for a buckyball containing some 1000 
nucleons), such as to guarantee that in all 
standard experiments — in particular, diffrac- 
tion experiments — one observes interference 
effects. For macroscopic objects instead, the 
spread is very small (cr 9 (oo) ~ 3 x 10~ 14 m, 
for a lg object), so small that for all practi- 
cal purposes the wave function behaves like a 
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point-like system. Once again, this is how 
collapse models are able to accommodate 
both the "wavy" nature of quantum systems 
and the "particle" nature of classical objects, 
within one single dynamical framework. One 
should also note that, byproduct of the 
collapse in position, one has an almost per- 
fect collapse in momentum, compatibly with 
Heisenberg's uncertainty relations. 



able result has been proven: with probabil- 
ity one, any initial state converges asymptot- 
ically to a Gaussian wave function, having 
a fixed spread both in position and in mo- 



Eqn. (28) says that the mean momentum 
undergoes a diffusion process. For micro- 
scopic systems, such a diffusion is appreciably 
large: the wave function is kicked back and 
forth by the collapse noise. For larger ob- 
jects instead, the diffusion becomes weaker 
and weaker, to the point that at the macro- 
scopic level it is almost entirely negligible. 
The same is true for the mean in position, 



mentum, given by Eqn. (33) and (34) respec- 
tively. The collapse process not only localizes 
wave functions, but also smoothes all their 
bumps and eventually shapes them as Gaus- 
sian functions. The second comment is that 
the above results refer only to a free particle. 
Also the harmonic oscillator can be treated 
in a fully analytical way, but more general 
potentials require perturbative approaches, 
which have not been explored so far. 

To conclude this section, let us consider 
the many-particle equation: 

N 



according to Eqn. (27). In this way, collapse dipt 



models can explain both the stochastic na- 
ture of quantum phenomena and the (appar- 
ently) deterministic nature of classical ones. 
Moreover, the average momentum E[(p) t ] is 
constant ((p) t = hk t ), while the average posi- 
tion is given by E[(g) t ] = E[(p) t ]/m: the par- 
ticle, on the average, moves along a straight 
line, depending on its initial momentum. 

Two comments are in order. The first one 
is that the above results refer only to the 
special case of Gaussian wave functions, like 



i 

K 



(0 



i=l 



1 N 



1=1 



1>t 



(36) 



that of Eqn. (25). However, in (Bassi and 



Duerr 


2008 


Bassi et al. 


2010b 



where H is the quantum Hamiltonian of 
the composite system, the operators (i = 
1, . . . N) are the position operators of the par- 
ticles of the system, and (i = 1, . . . N) 
are N independent standard Wiener pro- 
cesses. 

As often in these cases, it is convenient 
to switch from the particles' coordinates 
(xi, x 2 , . . . xn) to the center-of-mass (R) and 



36 



relative (xi, $2, . . . x^r) coordinates: 

( N N 



R 



-E 

M ^ 

i=i 



i=i 



(37) 



let Q be the position operator for the cen- 
ter of mass and (i = 1 . . . AT) the position 
operators associated to the relative coordi- 
nates. It is not difficult to show that — under 
the assumption H = H CM + H iel — the dy- 
namics for the center of mass and that for 
the relative motion decouple; in other words, 
ip t ({x}) = tpf M (R) ® ipf({x}) solves Eqn. 



(36) whenever ifjf M (R) and ^ o1 ({x}) satisfy 
the following equations: 



% 

K 



N 



H rcl dt +J2V^i(Q t ~ (q i )t)dW 1 



(0 



i=i 



1 N 

~o 5Z Ai( ^ ~ (Qi)t) 2 dt 



i=i 



(3t 



and: 



h 



H CM dt + ^/\Z(Q - (Q) t )dW t 



(Q - (Q)i 



) 2 dt 



CM 



(39) 



with: 



N 



Ar 



71=1 



M 
m 



A . 



(40) 



the system behaves quantum mechanically, 
modulo tiny modifications given by the col- 
lapse process. The second equation describes 
the center-of-mass evolution, and here we can 
see once again the most important feature 
of collapse models: the amplification mech- 
anism. The collapse strength of the cen- 
ter of mass is proportional to the size (i.e. 
the number of constituents) of the system. 
For microscopic systems, A CM is similar to 
Ao, i.e. very weak; in these cases, the col- 
lapse is almost negligible. For macroscopic 
objects, A CM (~ iVA , with iV ~ 10 24 ) can 
be very strong, implying a rapid and efficient 
collapse of the wave function. It is precisely 
because of the amplification mechanism that, 
with a single choice of Ao, one can describe 
both quantum and classical phenomena at 
the same time. 

1. The quantum formalism derived 

Collapse models contain a unique and uni- 
versal dynamics, which applies to all phys- 
ical situations. Measurements play no spe- 
cial role in collapse models. It then becomes 
interesting and important to show how the 
entire phenomenology of quantum measure- 



The first of the above equations describes the ments emerges from the universal dynamics 
internal motion of the composite system: it of collapse models. To do this, we use the 
basically says that the internal structure of QMUPL model, because of its relatively sim- 
ple mathematical structure. We will show 
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that measurements always have a definite suring apparatus A, which is general enough 

outcome, are randomly distributed according to describe all interesting physical situations: 

to the Born rule, and that at the end of the we assume that the apparatus consists of 

measurement process, the wave function of a fixed part plus a pointer moving along a 

the micro-system collapses according to the graduated scale, in such a way that differ- 

von Neumann projection postulate. All these ent positions of the pointer along the scale 



features are included in Eqn. (23), without correspond to different possible outcomes of 

any need for extra axioms. the measurement. To simplify the analy- 

m, , . , sis, we study the evolution of the center of 

Ihe measurement setup we consider con- 

. , £ . . . . ,. mass of the pointer only, and disregard all 

sists of a microscopic system o interacting 

.,, . , j i'i j ti other macroscopic and microscopic degrees 

with a macroscopic system A, which acts like 

n , of freedom; accordingly, the pointer will be 

a measuring apparatus; both systems are de- 

., , . , . , , „ T treated like a macroscopic quantum particle 

scribed m quantum mechanical terms. We as- 

, , . n , n ., of mass m moving in one dimension only, 

sume that the measurement includes a finite 

r Ai-i j.i , whose state space is described by the Hilbert 

set of outcomes. Accordingly, we assume that 

sp 

We assume that the wave function of the 



the microscopic system S can be described by s P ace ^ A ^ 



a finite-dimensional complex Hilbert space. 
For the sake of simplicity, and without loss of P ointer of A is sub J ect to a spontaneous col- 



generality, we can consider the simplest case: la P se P rocess according to Eqn. g3p, while 
U s = C 2 , because the generalization of what the wave Unction of the microscopic sys- 
follows to C n is quite straightforward. Since tem S evolves according to the standard 
the most general self-adjoint operator O act- Schrddinger equation, since, as typical of dy- 
ing on C 2 can be written as namical reduction models, the stochastic col- 
lapse terms have negligible effects on micro- 

= o + |+)(+|+o_|-)(-|, (41) 

scopic quantum systems, .tor defimteness, let 
where |+> and |-> are the eigenstates of O, ug consider a pointer of mass m = i g ( Le ., 

while o + and o_ are its two real eigenvalues, a pointer made of an Avogadro number of 

for definiteness and with no loss of generality, nuc i eons ) 

in what follows we will take o± = ±h/2 and ,„ , , , , TT ,,, . TT , , r 

We take the total Hamiltoman H to be of 

O to be the z-component of the spin, S z , of £ 

F F ' ' the form: H = H s + H A + i^iNT- The first 

a spin 1/2 particle. . TT r 

' term is the quantum Hamiltoman tor the mi- 

We take the following model for the mea- croscopic system: we assume that the time 
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scale of the free evolution of the microscopic initial state of the micro-system S: 
system is much larger than the characteristic 

time scale of the experiment ("instantaneous t c +l~'~) + c -l )] ® 00 

measurement" assumption); accordingly we c + |+) £g> + + c_ |— ) £g> 0_, (43) 

take H$ to be the null operator. The see- 
where <p± are final pointer states spatially 

ond term is the quantum Hamiltonian of the 

translated with respect to the initial state 4>o 

pointer, which we take equal to that of a non- 

by the quantity ±{h/2) kT. We will see how 

relativistic free quantum particle of mass m: 

collapse models modify this linear evolution. 

= p /{2m), where p is the momentum 

operator. Finally, we assume the interaction The strength of the coupling constant k 

term H mT between the two systems to be of has to be chosen in such a way that the dis- 

the von Neumann type, devised in such a way tance kT between the initial state O of the 

as to measure the spin S z : pointer and any of the two final states <$>± is 

,, N . T , . macroscopic: for definiteness, let us choose 

H mT (t) = kAJ S z ®p, (42) 1 

Tlk — 1 cm s -1 , so that HkT = 1 cm. 

where k is a coupling constant and A T : t h-> 

a t • rp r j|4l 4.- i We take the initial states of the micro- 

is a j -normaiizedjj non-negative, real 

, j f r , • .j n i scopic system S and of the macroscopic appa- 

vaiued, function of time, identically equal r J r rr 

, , . j . . , i f ratus A to be completely uncorrelated, as it 

to zero outside a given interval of the ^ J 

( /, , i rri\ i -j j_x, j_* - A is customary and appropriate for the descrip- 

form (to, to + 1 ), i.e., outside the time inter- J ^ r r 

i r, ,i m rp i j • i • i tion of a measurement process. Accordingly, 
val of length 1 , say 1 = Is, during which the ^ & J 

, , i , i ,• we assume the initial state of the total system 

experiment takes place; we choose the time J 

origin in such a way that the experiment be- ^ + ^° be. 

gins at to = 0. As is well known in standard [ c +l+) + c -\~ )} ® 0o, (44) 

Quantum Mechanics, Hint generates the fol- , , , , , « , „ , , f , , 

^ ' ° where <po describes the ready state of the 

lowing type of evolution, depending on the ■ , A 

& J r & macroscopic apparatus A. 

Some considerations are in order, regard- 

4 By a T-normalized function, we simply mean . . . . 

mg the initial state <po of the pointer. In the 

r+oo i-to+T 

/ Aj dt = Afdt = T. previous section we showed that, according 

J — oo J to 

to Eqn. (23), the wave function for the cen- 
Note that depends also on the initial time to; — 

we will omit indicating this explicitly, when no con- tre of mass of an isolated quantum system 
fusion arises. reaches asymptotically (and very rapidly, for 
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a macro-object) a Gaussian state of the form 

l-i 



lG, 



X) 



1™\ 



exp 



(x — x t ) + ik t x 
(45) 



whose parameters at G C, and x t , k t , % , 
t G 1 satisfy a complicated set of non- 



linear stochastic differential equations (Bassi 



(modulo a time-dependent global phase fac- 



and Salvetti, 2007), with given initial condi- 



tor) with a q defined as in Eqn. (33), taking 
the value: a q ~ 4.6 x 10~ 14 m for m = lg. 
The dispersion of the Gaussian function of 
Eqn. (45) in momentum space is a p ~ 1.6 x 
10~ 21 kg m s _1 , as described in Eqn. (34). 



tions. In particular: 

= In |c±| (49) 
(of course we now assume that c± ^ 0). 



it 



In order to extract the relevant physical 
information, let us consider the differences 



In our measurement model, we assume 
that the pointer is isolated for the time prior 
to the experiment; during this time its wave 
function converges rapidly towards a state 

close to (45), which we therefore assume to be can eagily prove that ^ and ^ satisfy a get 



X t := xf — xj and K t := kf — k^, which rep- 
resent the distance in position and (modulo 
h) momentum space between the centres of 
the two Gaussian functions <pf and <f)J . One 



the initial state of the pointer. To summarize, 
we take as the initial state of the composite 
system S + A the vector: 

= [c+|+> +c_|+)] ®0 G . (46) 

We choose the natural reference frame where 
the pointer is initially at rest, so that k = 
m _1 , with the origin set up in such a way that 
Xq = m. 

We are now ready to solve the collapse 
equation. It is not difficult to show that, for 
the given initial condition, the solution takes 
the form 

z/>t= + (47) 



of linear and deterministic equations (Bassi 



and Salvetti, 2007): 



d 
dt 







—uj h/m 




x t 


+ 


HkAJ 


K t 




-2A 




K t 








(50) 



where 4> t have the form: 



exp \-a t 



x x t ) 



+ikfx+ 1 t+i6t] 
(48) 



where both the non-linear and the stochas- 
tic terms cancel out. The solution of the 
above system depends of course on the spe- 
cific choice for the function A^; a simple rea- 
sonable choice is the following: 

l t e [0,21 

1 (51) 

else, 

which — in a standard quantum scenario — 
means that during the measurement each 
term of the superposition moves with a con- 
stant speed towards the left and towards the 
right, respectively. According to this choice, 
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X t , given the initial condition Xo 
evolves in time as follows: 

for < t < T 



m, of one of the two Gaussian wave functions 



X, 



e~ wt ' 2 sin -t 

2 



2hK 



-Lit/2 



sin -t - e^ /2 sm-(t-T) 
2 2 V ' 



for t > T. 

(52) 

Since u^ 1 ~ 2.0 x 10 4 s is a very long time 
compared to the measurement-time, we can 



meaningfully expand Eqn. (52) to first order 
in ut: 

{tint for < t < T = (Hk)- 1 = 1 s, 
lcm for T < t < uj- 1 ~ 2.0 x 10 4 s 

(53) 

As we see, the distance between the two 
peaks increases almost linearly in time, reach- 
ing its maximum (1 cm) at the end of the 
measurement process, as predicted by the 
standard Schrodinger equation; after this 
time, their distance remains practically unal- 
tered for extremely long times, and only for 
t ~ 2.0 x 10 4 s it starts slowly to decrease, 
eventually going to 0. Note that such a be- 
havior, being determined by uj, does not de- 
pend on the mass of the pointer, thus a larger 
pointer will not change the situation. The 
moral is that X t behaves as if the reduction 
mechanism were not present (as if Ao = 0) so 
we have to look for the collapse somewhere 
else. 

As we shall see now, the collapse occurs 
because, in a very short time, the measure 



(<% or 07) becomes much smaller than the 
measure of the other component; this implies 
that one of the two components practically 
disappears, and only the other one survives, 
the one which determines the outcome of the 
experiment. Of course, this process is ran- 
dom and, as we shall prove, it occurs with 
a probability almost equivalent to the Born 
probability rule. 

The relative damping between the two 



Gaussian components of Eqn. (47) is mea- 
sured by the stochastic process: 

r 4 = 7 t + -7r- (54) 



Note that, according to Eqn. (49), e r ° = 
|c + |/|c_|. If, at the end of the measurement 
process, it occurs that r t 3> 1, it means 
that 0^ is suppressed with respect to <fif, so 



that the initial state (46) practically evolves 



to |+)®0^; the opposite happens if T t — 1. 



In (Bassi and Salvetti, 2007) it is shown 



that T t satisfies the nonlinear stochastic dif- 
ferential equation: 

oT t = XX? tanhT t dt + VXX t dW u (55) 

to be solved with initial condition T = 
In |c + /c_|. To proceed further with the anal- 
ysis, it is convenient to perform the following 
time change: 

t — > s t := A f X 2 dt', (56) 
Jo 

which allows us to describe the collapse pro- 
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cess in terms of the dimensionless quan- Eqn. (55) reduces to: 
tity s that measures its effectiveness. Using 



tanh T s ds + dW s 



(59) 



Eqn. (52), one can solve exactly the above 



integral and compute s function of t. 
Such a function however cannot be inverted 
analytically in order to get t from s. There- 



this equation belongs to a general class of 
stochastic differential equations whose prop- 



erties are known in detail (Gikhman and Sko- 



fore, we use the simplified expression (53) 



rokhod, 1972). Here we report the main re- 



in place of the exact formula (52) to com 



pute the integral, an expression which, as we 
have seen, represents a very good approxima- 
tion to the time evolution of X t throughout 
the whole time during which the experiment 
takes place. Accordingly, we have: 

\h 2 K 2 



suits. 

Collapse time. The collapse time is the 
time when |T S | > b (where b is some fixed 
number much larger than 1), i.e. the time 
when one of the two terms of the superpo- 
sition becomes dominant with respect to the 
other: 



s t 



t = t s ~ 



t 6 ~ 2.0 x 10 17 (t/s) 

o 

for <t < T = 1 s, 



(57) 



'COL 



s ~ (1.7 x 10" 



s s 



\h 2 n 2 

for < s < \h 2 K 2 /3 = 2.0 x 10 17 . (58) 

Note that, according to the above equations, 
the physical time t depends on s through the 
inverse cubic-root of A, i.e. on the inverse 
cubic-root of the mass of the pointer; this 
time dependence of t on A is important since, 
as we shall see, it will affect the collapse time. 
We do not study the functional dependence 
between s and t for t > T since, as we shall 
soon see and as we expect it to be, the col- 
lapse occurs at times much smaller than T. 



inf{s : |r s | > b}. (60) 

This is a random variable (for each run of 
the experiment, the collapse time slightly 
changes), whose mean and variance can be 



exactly computed (Gikhman and Skorokhod 



1972). In particular, if we start with an 



equal- weight superposition (Tq = 0), then: 
EpfScoiJ - b and Vp^col] - b, where E P [-] 
and Vp[-] denote the mean and variance, re- 
spectively. If we transform back from s to 
the physical time t, we have the following 
estimate for the collapse time for the 1-g 



pointer (Bassi and Salvetti, 2007): 



1.5 x 10 _4 s. 



(61) 



Written in terms of the new variable s, 



(This value refers to b = 35.) The collapse oc- 
curs within a time interval smaller than the 
perception time of a human observer. More- 
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over, as proven in (Bassi and Salvetti 2007) 



T COL is proportional to the inverse cubic-root 
of the mass of the pointer: therefore, the 
bigger the pointer, the shorter the collapse 
time. With our choice for Ao, even for a 1-g 
pointer the reduction occurs practically in- 
stantaneously. 

It is important to note that, at time 
T CO l — 1-5 x 10~ 4 s, the distance between 
the two Gaussian components, according to 
Eqn. (53), is approximately X TcOL ~ 1.5 x 
10 -4 cm: this means that, with very high 
probability, the collapse occurs before the two 
components have enough time to spread out 
in space to form a macroscopic superposition. 
This means that, from the physical point of 
view, there is no collapse of the wave func- 
tion at all, since it always remains perfectly 
localized in space at any stage of the experi- 
ment. 

Collapse probability. Let us call P + the 
probability that T s hits the point +b before 
the point —6, i.e. the probability that <f)f 
survives during the collapse process so that 
the outcome of the measurement is "+/i/2". 
Such a probability turns out to be equal 



outcome of the experiment is u —h/2", is: 

1 tanh b — tanh T 
2~ 



(63) 



tanh b 

By taking into account that tanh b ~ 1, since 
we have assumed that > 1, and resorting 



to Eqns. (49) and (54), we can write, with ex- 



tremely good approximation (Bassi and Sal 



yettij [2007J): 

P \ ~ - 



[1 + tanhr 

e 27 + 
e 2 7 + _j_ e 2 7( 7 

- [1 — tannic 



:X(> 



\c+ 



e r ° + e~ r o 

|2 



(64) 



-r 



e L ° + e~ 

|2 



(65) 



e 27(T 4- e 2 To 

We see that the probability of getting one 
of the two possible outcomes is practically 
equivalent to the Born probability rule. On 
the one hand, this is an entirely expected 
result, since collapse models have been de- 
signed precisely in order to solve the measure- 
ment problem and in particular to reproduce 
quantum probabilities; on the other hand, it 
is striking that a very general equation like 



Eqn. (23), which is meant to describe both 



to (Bassi and Salvetti, 2007): 



1 tanh b + tanh T 
F+ ~ 2 ^nh6 ; (62) 

while the probability P that T s hits the 
point —b before the point +b, i.e. that the 



quantum systems as well as macroscopic clas- 
sical objects (i.e. all physical situations, at 
the non-relativistic level), when applied to 
a measurement situation, provides not only 
a consistent description of the measurement 
process, but also reproduces quantum prob- 
abilities with such a good precision. 
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State vector after the collapse. At time teraction is that the microscopic superposi- 

t > T COL the state of the composite system tion of the spin states of the micro-system is 

is: rapidly reduced right after the measurement. 

l\ _|_ e 2 ' v 1 Note finally that, after the collapse, the 
where e t = e^t-iT) and the normalized states of the micro-system and of the pointer 
Gaussian states 4>f are defined as follows: are de facto factorized: as such, after the 
measurement process one can, for all practi- 



2na 2 q 



cxp 



(x — X 



±\2 



, ( cal purposes, disregard the pointer and focus 

only on the micro-system for future experi- 

+ikfx + % Of] . (67) 

ments or interactions with other systems, as 



Let us assume that the collapse occurred in is the custom in laboratory experiments, 
favor of the positive eigenvalue, i.e. in such a 

lo conclude, we have seen how collapse 

way that T t > b for t > T COL ; it follows that: , , , „ 

models can describe quantum measurements 

ei<e fe ~0 if6^>l, (68) in a precise way, without ambiguities and 

and we can write, with excellent accuracy: paradoxes. We have seen that the stan- 
dard recipe for quantum measurements (def- 

% ~ \+)®4>t- (69) 

inite outcomes, the Born rule, the postu- 
We recover in this way the postulate of wave kte of wave function co n apse ) derives from 

packet reduction of standard quantum me- the dynamical equati on pfc and need not 



chanics: at the end of the measurement pro- be postulated in an ad hoc way . gut there 
cess, the state of the micro-system reduces ig somet hing more: it can be shown QBassi 
to the eigenstate corresponding to the eigen- 



et al. , 2007) that also the Hilbert space oper- 



value which has been obtained as the out- ator forma i ism _ according to whic h observ- 
come of the measurement, the outcome be- able quantities are re p resen t ed by self-adjoint 
ing defined by the surviving Gaussian com- operatorS) whose e igenstates and eigenval- 
ponent (0+ in this case). Note the impor- ueg haye the role ascribed to them by stan _ 
tant fact that, according to our model, the dard quantum mec hanics-can also be de- 
collapse acts directly only on the pointer of riyed from Eqn pK In other WQrdS) in cq1 _ 



the measuring apparatus, not on the micro- lapge modelg there ig Qnly the waye function 
system; however, the combined effect of the and a coUapge equation like Eqn . gg. every . 
collapse plus the von Neumann type of in- thing elge can be deriyed from it 
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E. On the origin of the noise field 

As we mentioned earlier, the QMUPL 
model can be generalized in order to include 



dissipative effects (Bassi et al. , 2005b) and 



non-white noises ( [Bassi and Feriddi[ |2009a]b 



Ferialdi and Bassi, 2011 ). According to these 



models, the noise field acquires a more physi- 
cal character: it can be assigned a finite tem- 
perature, and its spectrum is arbitrary; more- 
over, this noise is assumed to fill space. It 
becomes then natural to consider whether it 
can have a cosmological origin. At present it 
is too early to answer such a question, though 



some work has already been done (Adler 



and Bassi, 2008) and some people have al- 



ready suggested that it could have a gravi- 



tational (Diosi 1989 Feynman et al., 1995 



Karolyhazy et al, 1986 Penrose, 1996) or 



pre-quantum (Adler, 2004) nature. More 



over, it is still not clear why it has an anti- 
Hermitian coupling to matter, which is nec- 
essary to ensure the collapse^} However, one 
can meaningfully ask whether a noise with 
'typical' cosmological properties (in terms of 
temperature and correlation function) can in- 
duce an efficient collapse of the wave func- 
tion, where by 'efficient' we mean that the 
collapse is fast enough to avoid the occur- 



5 With an Hermitian coupling, one would have a 
standard quantum Hamiltonian with a random po- 
tential; the equation would be linear and no sup- 
pression of quantum superpositions would occur. 



rence of macroscopic superpositions. 

Let us consider a Gaussian state as in 



Eqn. (25), whose time evolution can be an- 



alytically unfolded in all models so far de- 
scribed. The quantity we are interested in, is 
its spread in space a t = {2^J~o^)~ 1 ^ 2 . This is 
plotted in Fig. 1. The top row shows the dif- 
ference in the evolution of the spread as given 
by the dissipative QMUPL model and by the 
original QMUPL model. Two temperatures 
have been considered: T = 2.73K (top row 
left), as for the CMBR and T = 2.73 x 10~ 3 K 
(top row right). In both cases the differ- 
ence is negligible. This means that even a 
rather cold thermal field, according to cos- 
mological standards, can induce an efficient 
collapse of the wave function, as efficient as 
with the standard QMUPL model. The bot- 
tom row shows the difference in the evolu- 
tion of the spread as given by the QMUPL 
model, and by the colored-noise model, with 
a noise having a frequency cutoff. Neither 
high-frequency cutoff affects the collapsing 
properties of the model in an appreciable way. 
These results can be compared with the be- 
havior of typical cosmological fields such as 
the CMBR, the relic neutrino background 
and the relic gravitational background. The 
spectrum of the first two have a cutoff (mea- 
sured or expected) at ~ 10 11 Hz, while the 
spectrum of the third one probably lies at 



10 10 Hz (jGrishchuk, 2010). All these cut 
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Differences in the spreads a after t = 10 2 s 



I^qmupl - ^d-qmuplI 
temperature T = 2.73 K 




|<7qmupl — o"d-qmupl| 

temperature 

T = 2.73 x 1(T 3 K 




I^QMUPL — <7CN-QMUPL| 

^^^^ frequency cut-off 

7 = io 10 s- 1 



2 3 4 5 6789 10 234 5 6789 10 

log(n), n =numbcr of particles 

FIG. 1 Difference between the spread a pre- 
dicted by the QMUPL model with that given by 
the dissipative QMUPL model (D-QMUPL), for 
T = 2.73K (top row left) and T = 2.73 xl(T 3 K 
(top row right); and by the colored-noise model 
(CN-QMUPL), with cutoff at 10 10 Hz (bottom 
row left) and 10 2 Hz (bottom row right). As the 
color bars on the right show, the whiter the re- 
gion, the greater the difference in the spreads. 
The initial value <7o = 5 x 10~ 7 m and the elapsed 
time t = 10 _2 s reproduce the typical geometry of 
the macromolecule diffraction experiments. At 
lower temperatures or lower cutoffs the wave 
function tends to collapse more slowly, which re- 
sults in a bigger difference with respect to the 
QMUPL model. Regarding the plots in the top 
row, the discrepancy manifests in the lower left 
corner of the plot for T = 2.73 x 1CT 3 K which 
disappears for T = 2.73K. Regarding the plots 
in the bottom row, it is manifest in the diag- 
onal strip for 7 = 10 2 Hz which decreases for 
7 = 10 10 Hz. This diagonal feature exactly lies 
on the ridge between the quantum and the clas- 
sical regime. The large discrepancy there is due 
to the missing high frequencies in the noise spec- 
trum of the colored-noise model. Without these 
high frequencies the colored-noise model is not 
able to reproduce the sharpness of the ridge pre- 
dicted by the QMUPL model. 



noise ensure a rapid collapse of the wave func- 
tion. While the collapse is robust over a large 
range of cutoffs, other effects, such as the 
emission of radiation from charged particles, 
highly depend on the spectrum of the noise 



correlator (Adler and Ramazanoglu, 2007). 



Therefore, the message which can be 
drawn is that a cosmological field with 'typ- 
ical' properties can induce an efficient col- 
lapse. A great challenge is to test the ex- 
istence of such a field. 

F. The CSL model 

The QMUPL model has the advantage of 
allowing for quite a rigorous mathematical 
analysis of the main features of collapse mod- 
els, as was shown in the previous sections. 
However it does not seem physically realis- 
tic, for two main reasons. The first one is 
that it is built for systems of distinguishable 
particles, and its generalization to identical 
particles does not seem straightforward. The 
second reason is that the noise field depends 
only on time, not on space, thus it cannot be 
immediately identified with a random field of 



Nature. The CSL model (Pearle, 1989a), 



(Ghirardi et at, 1990c) overcomes the above 



difficulties, and so far remains the most ad- 
vanced collapse model. In its mass propor- 



offs as well as that of ~ 10 Hz proposed tional version (Pearle and Squires, 1994), it 



in (Bassi and Ghirardi, 2003) for the collapse is defined by the following stochastic differ- 
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ential equation in the Fock space: ator, instead of the mass density operator]. 

dipt 



— —Hdt (70) As one can see fr° m Eqn. (70), in the CSL 

h 

^ . model the collapse operators are the density 

mo J ^ x (^( x ) (Mfaj/tjdWtfe) number operator r/>J (y) ipj (y), which means 



7 



2m\ 



dx(M(x) - (M(x)) t ) 2 dt 



that superpositions containing different num- 
bers of particles in different points of space 
as usual, H is the standard quantum Hamil- are suppressed. This is equivalent to collaps- 
tonian of the system and the other two terms ing the wave function in space, in a second- 
induce the collapse of the wave function in quantized language, 
space. The mass mo is a reference mass, 

The collapse occurs more or less as in the 

which as usual is taken equal to that of a nu- 

QUMPL model, though it is more difficult to 

cleon. The parameter 7 is a positive coupling 

unfold: an easier and more handy way to look 

constant which sets the strength of the col- 

at the collapse is through the density matrix, 

lapse process, while M(x) is a smeared mass 

in particular how its off-diagonal elements de- 

aensity operator: 

cay in time. Since ordinary matter is made 
M(x) = ^^mjiVj(x), (71) j us t °f electrons and nucleons and, according 



3 to Eqns. (70) and (71), the collapse effect on 

^Y?'( x ) = j dygiy — x)'0j(y)'0j(y), (72) electrons is negligible in comparison to the ef- 
fect on nucleons, we focus our attention only 
^j(y), V>j(y) bein g, respectively, the creation ^ nucleong 

and annihilation operators of a particle of 



type j in the space point y. The smearing According to Eqn. (|70j), the decay of the 



function g(x) is taken equal to off-diagonal elements of the density matrix 

1 -x 2 /24 / 7 on Pt = E [\A)(il>t\] (where E[-] denotes the 



(v^tt^c) stochastic average) of a many-nucleons sys- 

where r c is the second new phenomenological tern, in the position basis, is ( |Ghirardi etal. 
constant of the model. Wt (x) is an ensem 



1990c): 



ble of independent Wiener processes, one for d ... .... _„.,_,. .... 

— (x |p t |x ) = -r(x,x ) (x \pt\x ), (74) 

each point in space. [In the original, i.e. 'not 

^ OT j 1 .I 1 where x' = x',,x',,...x' lv (and similarly for 

mass proportional , LbL model, the integrals 1 ^ JV v 

1 j x"). In the above equation, we have neglected 

are proportional to the number density oper- ' 

the standard quantum evolution. The decay 
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function T is: 



r 2 ? X ^ + G ^ X * X ^ to Eqn. d76b their contribution to T is negli 



in a superposition are displaced by a dis- 



tance 



|x' — x"| <C r c , then according 



~ 2G(x^ - x" 



(75) 



where the indices i,j run over the A" nucleons 
of the system, and: 

1 



G(x) 



-x 2 /4rj 



(76) 



(47rr2,) 3 /2 

As a first observation, in the case of a sin- 
gle nucleon T reduces to: 



r(x',x"i 



7 



(47rr c ) 3 / 2 



1 _ e Hx'-x"| 2 /4r£ 



(77) 



gibly small. Thus, only superpositions with 
£ > rc contribute to V and trigger the col- 
lapse of the wave function. In such a case, 
the formula shows that for groups of particles 
separated (in each term of the superposition) 
by less than rc, the rate V increases quadrat- 
ically with the number of particles while for 
groups of particles separated by more than 
rc it increases linearly. Thus we have the 
following simplified formula for the collapse 



which is precisely the GRW 1-particle col- rate (lAdlerl 120071): 



lapse term (see Eqn. (22)). Accordingly, the 



A CSL n 2 Ar, 



two models give similar predictions regard- 
ing the collapse effects on systems containing 
just a few particles, while for many-particle 
systems important differences emerge, as we 
will see soon. The collapse rate is defined in 
terms of 7 as follows: 



(80) 



A 



7 



C'SL 



(4vrr c )3/2 • 



(78) 



In (Ghirardi et at, 1990c) the following 



choice for 7 was made: 7 ~ 10 30 cm 3 s 1 , 
corresponding to: 

A CSL ~ 2.2 x lO^V 1 . (79) 

Note the difference of about one order of mag- 
nitude between A GRW and A CSL . 

Several useful approximate formulas can 



where n is the number of particles within a 
distance rc, and A^ is the number of such 
clusters. We note that the quadratic depen- 
dence of T on the number of particles — which 
is absent in the original GRW model — is a 
direct effect of the identity of particles. This 
means that the identity of particles works in 
favor of the collapse. 

An estimate for small distances can be ob- 
tained by Taylor expanding G(x) as follows: 



G(x) 



(4vrr c ) 3 / 



.2 1 



which leads to: 

r(x',x") ~ 




4r c 



- x,- 



(81) 



(82) 



be obtained from Eqn. (75). The first one 



is for large distances. When the particles 



As we can see — and as expected — the col- 
lapse strength grows quadratically with the 
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superposition distance for small distances, 
like in the GRW and QMUPL models, the 
important difference here being that there is 
a quadratic dependence also on the number 
of particles. In both cases (large- and small- 
distance approximation), we see the amplifi- 
cation mechanism: the collapse rate increases 
with the size of the system. 

Another useful formula can be obtained 
for macroscopic rigid systems, for which the 
mass distribution can be expressed by a den- 
sity function D(x), averaging the contribu- 
tions of the single nucleons. In such a case 
the decay function T takes the simpler ex- 



Further properties of the CSL model 



pression (Ghirardi et ai, 1990c): 



r(X\X") =7 J dx[D 2 (x) 

- D(x) J D(x + X / -X")], (83) 

where X' and X" are the positions of the cen- 
ter of mass of the object, in the two terms of 
the superposition. The physical meaning of 



Eqn. (83) can be understood by making ref- 



erence to a homogeneous macroscopic body 
of constant density D. Then the decay rate 
becomes: 

T = -?Dn olJT , (84) 
where n OUT is the number of particles of the 
body when the center-of-mass position is X', 
which do not lie in the volume occupied by 
the body when the center of mass position is 
X". 



are discussed in (Ghirardi et ai, 1990c) 



and (Bassi and Ghirardi 2003). 



1 . On the nature of the noise field of the CSL 
model 

As anticipated, contrary to the QMUPL 
model, in the CSL model the noise field can 
be given a straightforward physical interpre- 
tation. In order to see this, it is convenient 



to rewrite the CSL dynamical equation (70) 
in the following equivalent form: 



dip t 



+ 



h 

7 
2ml 



Holt 



(85) 



dx(M(x) - (M(x)) t )dW t (x) 



rfxdy(M(x) - (M(x))< 



G(x-y)(M(y)-<M(y)) 4 )* 



where M(x) is still the mass density operator 



defined in Eqn. (71), while A^ (x) now is the 
standard number density operator: iVj(x) = 
ip\ (x) ipj (x), and G(x) is the same Gaussian 



function defined in Eqn. (76). The Wiener 



processes W t (x) are not independent any- 
more, instead they are Gaussianly correlated 
in space, the correlator being G(x). The 
"white" noise field w(t, x) = dW t(x) / dt cor- 
responds to a Gaussian field with zero mean 
and correlation function: 

E[w(t,x)w(s,y)] = 5(t- s)G(x-y). (86) 
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As we see, the noise field of the CSL model 
can be interpreted clS 8b (classical) random 
field filling space. It is white in time, for the 
very simple reason that white noises are eas- 
ier to analyze mathematically It is Gaus- 
sianly correlated in space, with correlation 
length equal to r c . As discussed in connec- 
tion with the QMUPL model, it is tempt- 
ing to suggest that such a field has a cosmo- 
logical nature, and preliminary calculations 
show that a noise with typical cosmological 
features yields a satisfactory collapse of the 



ter equation of the CSL model coincides with 
the 1-particle equation of the GRW model. 
This type of generalization should not present 
particular problems, however it has not been 
worked out so far. 

The second generalization consists in re- 
placing the white noise field with a more 
general Gaussian noise. As for the QMUPL 
model, the equations become non-Markovian, 
therefore difficult to analyze mathematically. 
The analysis has been carried out to the lead- 
ing perturbative order — with respect to the 



wave function ( |Bassi et al.\\2010a$ . However collapse parameters 7— in QAdler and Bassi 



at this stage this is only a (very tempting) 
speculation. Moreover, one has to justify the 
non-Hermitian coupling and the non-linear 
character of the collapse equations, which are 
necessary in order to obtain the effective col- 
lapse dynamics. 

2. Generalizations of the CSL model and rel- 
ativistic models 

Like the QMUPL model, the CSL model 
can also be generalized in several directions. 
The first generalization one would make, is to 
include dissipative terms, in order to cure the 
problem of energy non- conservation. In fact, 
in the CSL model also the energy increases at 
a steady rate, though such an increase is neg- 
ligible for all practical purposes. This can be 
easily seen by noting that the 1-particle mas- 



2007, 2008). The result of the analysis is the 



expected one: the collapse qualitatively oc- 
curs with the same modalities as in the white 
noise case, the rate depending on the cor- 
relation function of the noise. In particu- 
lar, the rate is robust against changes of the 
correlation functions, while other predictions 
are very sensitive to the form of the time 



correlator (Adler and Ramazanoglu, 2007). 



Much more work is needed in order to un- 
derstand the properties of the non-Markovian 
CSL model. 

The great challenge of the dynamical re- 
duction program is to formulate a consistent 
model of spontaneous wave function collapse 
for relativistic quantum field theories; many 
attempts have been proposed so far, none of 
which is as satisfactory as the non-relativistic 
GRW and CSL models. 
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The first attempt ( Ghirardi et al. 1990b 



Pearle, 1990) aimed at making the CSL 



model relativistically invariant by replacing 



to the noise which, when it is assumed to be 
white, is too violent, so to speak, and causes 
too many particles to come out of the vac- 



Eqn. (70) with a Tomonaga-Schwinger equa- uum. To better understand the situation, let 



tion of the type 
5a(x) 



h 



H(x) 



17) 



7 {£{x) — (C(x))) w(x) 



I (C(x) - (C(x))) 2 #x), 



where now the wave function is defined on an 
arbitrary space-like hypersurface a of space- 
time. The operator T-L(x) is the Hamiltonian 
density of the system (x now denotes a point 
in space-time), and C(x) is a local density of 
the fields, on whose eigenmanifolds one de- 
cides to localize the wave function. The c- 
number function w(x) is a random field on 
space-time with mean equal to zero, while 
the correlation function — in order for the the- 
ory to be Lorentz invariant in the appropri- 



ate stochastic sense (Ghirardi et al. , 1990b) 



must be a Lorentz scalar. And here the prob- 
lems arise. 

The simplest Lorentz invariant choice for 
the correlation function is: 

E[w(x)w(y)\ = 5 {4) (x - y), (88) 

which however is not physically acceptable 
as it causes an infinite production of energy 
per unit time and unit volume. The reason is 



us go back to the non-relativistic Eqn. (70): 



also there we basically have a white- noise pro- 
cess, which however is not coupled locally 
to the quantum field (rf(s, y)a(s, y), the cou- 
pling being mediated by the smearing Gaus- 
sian function appearing in the definition of 
iV(x). One can compute the energy increase 
due to the collapse mechanism, which turns 
out to be proportional to tq. Now, if we want 
to have a local coupling between the quantum 
field and the noise, we must set r c — > +0, in 
which case the energy automatically diverges 
also for finite times. 

The simplest way out one would think of, 



in order to cure this problem of Eqn. (87) 



is to replace the local coupling between the 
noise and the quantum field by a non-local 



one, as in the CSL equation (70); this pro- 



cedure would essentially amount to replacing 
the white noise field with a non- white one. 
In both cases we need to find a Lorentz in- 
variant function which either smears out the 
coupling or replaces the Dirac-<5 in the defi- 



nition of the correlation function (88). This 



however is not a straightforward task, for the 
following reason. 

One of the reasons why the third term 



that in Eqn. (87) the fields are locally coupled (7/2) (C(x) — (C(x))) appears in Eqn. (87) 
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is to guarantee that the collapse mechanism 
occurs with the correct quantum probabili- 
ties (for those who are experts in stochas- 
tic processes, the third term is such that 
the equation embodies an appropriate mar- 
tingale structure); if we change the noise, 
we then have to change also the third term, 
and it turns out that we have to replace it 



new way to tackle this problem has been re- 



cently proposed by D. Bedingham (Beding- 
harn||2011a|bfr . 



Dowker and Henson have proposed a spon- 
taneous collapse model for a quantum field 



theory defined on a 1+1 null lattice (Dowker 



and Henson, 2004 Dowker and Herbauts 



2004), studying issues like the non- locality of 



with a non-local function of the fields (Myr- 



vold et al. , 2009; Nicrosini and Rimini, 2003) 



But, having a non-local function of the 
fields jeopardizes the entire (somehow for- 
mal) construction of the theory based on the 
Tomanaga-Schwinger formalism, as the in- 
tegrability conditions are not automatically 
satisfied. More analysis is required, however 
it is very likely that the model will turn out 
to be inconsistent. 

What we have briefly described is the ma- 
jor obstacle to finding a relativistic dynamical 
reduction model. We want to briefly mention 
three research programs which try to over- 
come such an impasse. 

P. Pearle has spent many years in try- 
ing to avoid the infinite energy increase of 
relativistic spontaneous collapse models, e.g. 
by considering a tachyonic noise in place of 
a white noise as the agent of the collapse 



the model and the no-faster-than-light con- 
straint. More work needs to be done in try- 
ing to apply it to more realistic field theories; 
in particular, it would be important to un- 
derstand if, in the continuum limit, one can 
get rid of the divergences which plague the 
relativistic CSL model. 



In a recent paper (Tumulka 2006b), gen- 



eralizing a previous idea of Bell (Bell, 1987b), 



process (Pearle, 1999b), obtaining sugges- 



Tumulka has proposed a discrete, GRW-like, 
relativistic model, for a system of N non- 
interacting particles, based on the multi-time 
formalism with iV Dirac equations, one per 
particle; the model fulfills all the necessary 
requirements, thus it represents a promising 
step forward in the search for a relativistic 
theory of dynamical reduction. Now it is im- 
portant to understand whether it can be gen- 
eralized in order to include also interactions. 

Recently, a completely different perspec- 
tive towards relativistic collapse models 



tive results. Unfortunately, as he has re- emerged (Adler, 2004). If one assumes that 



cently admitted (Myrvold et al, 2009), this 



program was not successful. A promising 



the random field causing the collapse of the 
wave function is a physical field filling space 
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and possibly having a cosmological origin, 
then there is no need to make the equa- 
tions relativistically invariant. The noise 
would select a privileged reference frame — 
in pretty much the same way in which the 
CMBR identifies a preferred frame. Then 
the collapse equation, being a phenomeno- 
logical equation, need not be relativistic in- 
variant, being dependent on the noise. The 
underlying theory, out of which these equa- 
tions would emerge at an appropriate coarse- 
grained level, should respect the appropriate 
symmetries (Lorentz invariance or a possible 
generalization of it). The theory would ex- 
plain the origin of the noise field which, be- 
cause of the initial conditions, would break 
the relevant symmetry. This is only specula- 
tion at the moment. However it is a reason- 
able program, though difficult to carry out. 

G. Choice of the parameters 



Recently, a much stronger bound on the 



collapse rate has been proposed (Adler 



2007), namely: 



The choice (19) for A GHW (or (79) for 



A CSL ) and tq makes sure that collapse models 
agree with all observational evidence about 
the quantum nature of microscopic systems, 
while macro-objects are always localized in 
space and behave according to Newton's laws 
(within experimentally testable limits). It 



sets a quantum-classical threshold at ~ 10 



13 



nucleons (Bassi and Ghirardi 2003 Ghirardi 



A» 



2.2 x 10 



-8±2 S -1 



(89) 



corresponding to a threshold of ~ 10 5 nucle- 
ons. The underlying motivation is that the 
collapse should be effective in all measure- 
ment processes, also those involving only a 
small number of particles as happens in the 
process of latent image formation in photog- 
raphy, where only ~ 10 5 particles are dis- 
placed more than tq- In order for the collapse 
to be already effective at this scale, one has 
to increase the conventional CSL value A CSL 



by ~ 10 9±2 orders of magnitude. Eqn. (89) 



is also the value one has to take in order to 
make sure that a superposition of 6 photons^] 
reaching the human eye, collapses within the 



eye itself (Bassi et al, 2010a). However, the 
value A GRW makes sure that the collapse oc- 
curs before the signal reaches the brain and 



turns into a perception (jAicardi et oZ.|[1991J). 
Both values A 



A,, 



rMj ~ 2.2 x 10 -17 s -1 and 
2.2 x 10 _8±2 s _1 are compatible with 



known experimental data (Adler and Bassi 



2009). However, such a large discrepancy of 



et al, 1990c). 



~ 9 orders of magnitude shows that there 
is no general consensus on the strength of 
the collapse process and consequently on the 
scale at which violations of quantum linearity 

6 6 photons corresponds to the threshold of vision. 
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can be expected to manifest. Experiments 
will decide which value — if any — is the cor- 
rect one. 

III. UNDERLYING THEORIES 

It is beyond doubt that if Continuous 
Spontaneous Collapse models are the right 
way to resolve the problems of quantum the- 
ory, these models must derive from an under- 
lying physical theory based on new symmetry 
principles. Here, we review two such possibil- 
ities. Trace Dynamics is such a theory due 
to Stephen Adler and collaborators, which 
bears the same relation to quantum theory 
that thermodynamics bears to statistical me- 
chanics and the kinetic theory of gases. The 
other, which is not quite a theory yet, is the 
idea that gravity has a role to play in bringing 
about collapse of the wave-function. 

A. Trace Dynamics, Quantum Theory, 
and Spontaneous Collapse 

The commutation rules in canonical quan- 
tization are obtained starting from the Pois- 
son brackets in the classical theory. This may 
be somewhat unsatisfactory, since the classi- 
cal theory itself is supposed to be a limit- 
ing case of the fundamental quantum theory. 
In order to know the fundamental theory we 
should not have to know its limit. The quan- 
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turn commutations should be achieved in a 
more fundamental manner and this is what 
Adler's scheme of Trace Dynamics sets out 
to do. However, eventually the theory goes 
beyond deriving quantum theory from a more 
fundamental framework, and provides a plau- 
sible resolution of the quantum measurement 
problem which is experimentally testable. 

The physics of Trace Dynamics can be de- 
scribed in the following three well-laid out 
steps : 

(i) The classical theory, which is the New- 
tonian dynamics of Grassmann-valued non- 
commuting matrices; and which as a conse- 
quence of global unitary invariance possesses 
a unique non-trivial conserved charge of great 
significance [Sections 1.1.1 - 1.1.4 below]. 

(ii) The statistical thermodynamics of 
this theory, the construction of the canoni- 
cal ensemble and the derivation of equilib- 
rium. The derivation of an all important 
Ward identity, as a consequence of assumed 
invariance under constant shifts in phase 
space, from which there emerge, as thermo- 
dynamic averages, the canonical commuta- 
tion relations of quantum theory, the Heisen- 
berg equations of motion and the equiva- 
lent Schrodinger equation of quantum theory 
[Sections 1.1.5 - 1.1.7]. 

(iii) The consideration of Brownian mo- 
tion fluctuations around the above ther- 
modynamic approximation, the consequent 



non-linear stochastic modification of the 
Schrddinger equation, the resolution of the 
quantum measurement problem and deriva- 
tion of the Born probability rule [Section 
1.1.8]. 

The review below is based on the book by 



Adler (Adler, 2004). The interested reader 



should consult Adler's book and the refer- 
ences therein for further details. 



1. The fundamental degrees of freedom 

In this scheme the fundamental degrees of 
freedom are matrices living on a background 
spacetime with complex Grassmann numbers 
(or more precisely elements of a graded alge- 
bra Qc of complex Grassmann numbers) as 
elements. Grassmann numbers have the fol- 
lowing properties: 



8i6j-\-6j6i = 0, 
X =9 R + i"0j, 



01 = 0; 

{Xr,Xs} = 0. 



Thus we have a matrix field with the help of 
which to each spacetime point we can asso- 
ciate a matrix 

' Xu(x) Xi2(x) \ 

\ X2l(x) X22(X) 

M{x) = 

v / 

Some further important properties of these 
Grassmann (matrix-) elements are the follow- 



ing: 

(i) Product of an even number of Grassmann 
elements commutes with all the elements of 
Grassmann algebra. 

(ii) Product of an odd number of Grassmann 
elements anti-commutes with any other odd 
number product. Therefore we have two dis- 
joint sectors: 

Bosonic Sector- B: Consists of the identity 
and the Even Grade elements of the algebra: 

B = {I, XaXb, XaXbXcXd, ■■■■} 

Fermionic Sector-F: Consists of the Odd 
Grade elements of the algebra: 

F = {Xay XaXbXci }• 

Therefore, the fundamental degrees of free- 
dom of the trace dynamics theory are the ma- 
trices made out of elements from these sectors 

5/ G {M; M {j E B} 

X ie{M;M l3 eF} 

Bi,Xi e Gm 
where Qm is the graded algebra of complex 
Grassmann matrices. 

A general matrix can be decomposed as 

M = A 1 (e B I ) + A 2 (e xi) 

into Bosonic and Fermionic sectors. We note 
that being matrices the degrees of freedom 
of this scheme are non-commutative in na- 
ture. Dimensionality of the matrices can be 
arbitrary; however we work with finite dimen- 
sional matrices with the assumption that ev- 
erything done subsequently can be extended 
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to infinite dimensional matrices as well. 
Next, we define an operation Trace on this 
matrix field as follows 

Tr-.g M ^ g c , 

which is a map from the space of matrices 
[Gm] to the field of complex Grassmann num- 
bers [Qc] and is given by the sum of the di- 
agonal elements of a given matrix. There are 
some nice Trace properties satisfied by the de- 
grees of freedom of the theory 

TrBiB 2 = TrB 2 B 1 (90) 
TrxiX2 = ~Trx2Xi (91) 
TrBx = TrxB (92) 

One also has some interesting Trace tri-linear 
cyclic identities 

TrB 1 ((B 2 ,B 3 )) =TrB 2 ((B 3 ,B 1 )) 

= TrB 3 ((B 1 ,B 2 )) (93) 
TrB[xi,X2\ = Trxi[x2,B] 

= Tr X 2[xi,B] (94) 

Trxi{B,X2}=Tr{xi,B} X 2 

= Tr[ Xl ,X2]B (95) 
Tr X ((B 1 ,B2))=TrB2((x,B 1 )) 

= TrB 1 ((B 2 ,x)) (96) 
Trxi({x2,X3)) =Trx2((X3,Xi}} 

= 7V X 3«Xi,X2», (97) 

where (( )} can be commutator or anti- 
commutator. Next, we can verify that the 
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matrices have the following adjoint rule 

(Of ...Of y = i ■''■'■()■!; Of lf , 

where Qi is the grade (odd/even) of the ma- 
trix O. The anti-commutative feature of ma- 
trix elements induces non-triviality in the ad- 
jointness properties, as seen above. 

We now examine the dynamics of these de- 
grees of freedom and later will construct the 
statistical mechanics of a gas of such 'parti- 
cles' to find that the equations of quantum 
theory are identities, valid in the thermody- 
namic limit of this underlying theory. 

2. Classical dynamics 

We can construct a polynomial P from 
these non-commuting matrices (say O) and 
obtain the trace (indicated in bold) of the 
polynomial 

P = TrP. 

Trace derivative of P with respect to the vari- 
able O is defined as 

<5P 

SP = Tr-60, 

i.e. 5- variation in P should be written in a 
way that resulting 5- variation of O in each 
monomial sits on the right. Then terms com- 
ing on the left of SO are defined as the trace 
derivative. It should be mentioned that one 
always constructs P to be an even graded el- 
ement of Grassmann algebra. Moreover, 50 
and 5P are also taken to be of same type 



(Bosonic/Fermionic) as O and P respectively. Using the trace derivative we obtain the 

Thus |^ will be of the same type as O. equation of motion by extremizing the action 

For example, let with respect to variation in q r using the dif- 

p _ aq-qqq ferentiation technique described above 

5L d 5L 

be a polynomial, where A, B, C are operators §q r ~ dt Sq r ' 

which in general do not commute with each Above is a matrix equation; in component 

other or with the variable O. Then, form we can write down iV 2 Euler-Lagrange 

equations of motion 
5TrP = 5P = Tr[e AO BOCA(50) /(5L x dL 



+e OC CAOB(SO)}, \<Wij d (lr)ji' 

Further we define conjugate momenta as 



using trace properties (90), (91), (92) and 5~L 

i_r pr = 

(97). Hence, the trace derivative will be ^ r 



3. Lagrangian and Hamiltonian Dynamics 



rp Since the Lagrangian is Grassmann 

— = e AO BOCA + e cCAOB. 

OU even, momentum will be of same type 

Above, e X Y = +1 if either X or Y is bosonic, (Bosonic/Fermionic) as q r . In general the 
and it is equal to -1 if both X and Y are coordinates or the momenta do not commute 
fermionic. among each other, for these are all arbitrary 

matrices. The Trace Hamiltonian is obtained 

as 

H = ^^PrQr — L. 

r 

Therefore, the Hamiltonian equations of mo- 
Armed with these tools we write the La- tion are 

grangian of a theory as a Grassmann even _ _^ oH _ £ ^ ?gg\ 

8q r ' ' 8p r 

polynomial function of bosonic/fermionic op- 

where e r = ±1, depending upon whether r is 

erators {q r } and their time derivatives {q r } 

a bosonic or fermionic degree of freedom. 

with {q r ,q r € Gm}- We define the Trace La- 
grangian 

We define a generalized Poisson Bracket 

H{Qr},{(ir}}=TrL[{q r },{q r }}, 

over the phase space {q r ,p r f 
and subsequently the Trace Action . f 5B 5B 5A\ 



{A,B} = Tr^e r ( 

T' ^ 



5q r 5p r 8q r Sp r J 



S - / dtL. (10Q) 
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which satisfies the Jacobi identity the quantity Trace Fermion number 

{A, {B, C}}+{C, {A, B}}+{B, {C, A}} = 0. N = ^TrJ^h 



o — ' / j Lirj Vr\ 

For A[{q r },{p r },t] one can easily verify that, . g congervec i j e n = o 
dA 

A = -7^- + {A, Hj. (101) This conserved charge corresponds to U(l) 

We observe that the matrix dynamics ob- S au S e transformations of Fermionic degrees 
tained above is non-unitary in general. For 

operators that do not have explicit time de- q r — y exp{ia}g r : p r —y exp {ia}p r , 



pendence, (101 ) does not show a unitary evo- £ , , 

i ' for real and constant a and r 6 t . 

lution of the type ^ J £ . . 

Requirement of bilinear rermiomc sector of 

x r (t) = i[G,x(t)\. H and self adjoint kinetic part forces, 

However, for Weyl- ordered Hamiltonians, £? lrs = — B\ TS \ B 2rs = —B\ TS , 

trace-dynamics evolution and Heisenberg , + u . 

and p r = q<., resulting m 

unitary time evolution can be shown to be 
equivalent, on an initial time slice on which 



N = -iTr 



the phase space variables are canonical. which resembles the number operator for 

Fermionic degrees of freedom. 



4- Conserved parameters 



Consider the restriction to matrix mod- 
els in which the only non-commuting matrix 
quantities are the dynamical variables. Thus 
the trace Lagrangian and Hamiltonian are 



As an obvious result of (|101j), the trace constructed from the dynamical variables us- 



Hamiltonian itself is conserved, j n g only c-number complex coefficients, ex- 

pj = jjj, jjj = o (102) eluding the more general case in which fixed 

matrix coefficients are used. Then there is a 
Noether charge corresponding to a global uni- 
tary invariance possessed by the trace Hamil- 
tonian (or equivalently Lagrangian) i.e. 

H = Tr (PrVsB lrs + p r B 2rs q s ) H[{rfq r U}, {U ] q r U}} = H[q r , q r ], (104) 

+bosonic, (103) for some constant unitary matrix U. From 
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Moreover, for a Trace Hamiltonian restricted 
to a bilinear form in the Fermionic sector 
with a self-adjoint kinetic part 



(104) and (99), it can be shown that this For those canonical transformations, clearly, 



Ndether charge is {Qj Ga} = Q {m) 

^ ~ X^^ 1 "'^] ^^{QnPr}i (105) Alternatively we can interpret (108) by say- 



res reF 



which we call the Adler-Millard charge. This in S that G a ^ invariant under the action 

charge having the dimensions of action and °^ ^ Then, along these lines it can be 

being trivially zero in point particle mechan- shown that is Poincare invariant when the 

ics, makes all the difference between Trace trace Lagrangian is Poincare invariant, where 

Dynamics, and ordinary classical mechanics P° m care transformations are generated by 

of point particles where all position and mo- trace functional Poincare generators. These 

menta commute with each other. Note that generators are global unitary invariant when 

in (EM the individual {anti-}commutators the trace Lagrangian is Poincare invariant. 



take arbitrary values in time; yet the par- Hence > we can make use of this char S e C in 
ticular combination shown in this equation Poincare invariant theories, 
remains conserved ^ we consider a Lagrangian which has a 

fermionic Kinetic part given by 

Now, if the Fermionic degrees of freedom 
have the adjointness property of p r = q\ ^ kin = ^ r Ir^rsQs (109) 



and the Bosonic degrees of freedom are self- 



r,s£F 

where A rs is a constant matrix having the 

adjoint (or anti-self adjoint) then the con- , £ , , T ■ \ a 

v property (for real trace Lagrangian) A rs = 

served charge is anti-self adjoint and traceless A + ,-, ^ ■ ,■■>■> it. 

A [ sr , then C is still conserved but now it can 

C = —C\ TrC = 0. (106) have a self-adjoint part as well, 

Since C is the Noether charge correspond- C + & = — [q s qt,A rs \. (HO) 

ing to global unitary invariance of the matrix 



r,s£F 

Generically, for a continuous space- 
time based trace Lagrangian written 
in terms of trace Lagrangian density 
Ga = TrAC, with £({<?/(x)}, {d^qi(x)}) which is invariant un- 

der the following symmetry transformations 



model, it can be used to construct the gener- 
ator of the global unitary transformation 



{G A ,G s } = G [AiS] (107) 



as the algebra of the generators. Now we con- 



qi(x) -)■ qi(x) + a(x)Ai(x) 



sider those canonical transformations which d^qi{x) — > d^qi(x) + a{x)d ll Ai{x) 

have global unitary invariant generator G. +d fl a(x)A[(x), (111) 
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there is a local trace current under four space rotation 

5C 

Tr ^ Sd » Ql W ^ ^ ' x » ~* x 't* = x » + 

fOT Which «(^) = + (1/2) V E Xfrn^W, 

= 0, 

suggesting that there is a conserved charge for the antisymmetric infinitesimal rotation 

q_ f x j0^ parameter 9^ u . P° is the conserved trace 

Hamiltonian. P M and M* 1 " together form a 



For a global unitary and Poincare invariant 
theory the conserved charges are the follow- 
ing 

c = J d 3 xJ2(mpi-pm), (H2) 



complete set of Poincare generators. 



The matrix operator phase space is well 
P' 1 = J d^xT ^; (113) behaved. We can define on this phase space 



a measure 



M"" = J d 3 xM ^ (114) 

u ^ + , , t \ ' dpi = TT d(x r )^ n , (118) 

where the momentum conjugate to qnx) is 11 

r,m,n 

5 C 

Pl ^ X) = 5d oqi (xY (U5) where A = 0,1 for 

trace energy momentum density is , x _ f , v \i 

t-m" _ _ jy Q^ q .- fug) 

' ^ 5,9 ^ ^ v ' This measure is invariant under canonical 

Tf = diag(-l, 1, 1, 1) is the Minkowski met- transformations (Louiville's theorem), hence 
ric and under dynamic evolution of the system as 



time evolution is a canonical transformation 



M*"" = x»T Xv - x»T x » generated by dtU. 

+Tr J2sd^ l xZqm (U7) dp[{x r + Sx r }] = dn[{x r }\. 



lm 



Further, bosonic and fermionic measures can 
be separated 



is the trace angular momentum density, with 
Xim being the matrix characterizing the in- 
trinsic spin structure of the field q\ such that dfi = dfindfiF, 



which are separately invariant under the ad- 
jointness properties assumed above. 
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5. Canonical Ensemble rameters 

p = p(Tr XC; H, r; N, rj) 
where A, r and r\ are the Lagrange multipliers 



With the matrix equations of motion (J99J) conjugate to C, H and N respectively. One 

for time evolution in trace dynamics we study important aspect to note is that while H and 

the evolution of phase space distribution. We N belong to Q c , C E Gm- Hence, r,r) E Gc 

assume that a large enough system rapidly while A E Gm- The dependence of pon TrXC 

forgets its initial distribution and the time is motivated from global unitary invariance. 

averages of physical quantities are equal to If C has a self-adjoint part as well, one can 

the statistical averages over an equilibrium break it into its self-adjoint (sa) and anti-self 

ensemble which is determined by maximizing adjoint (asa) parts, 

combinatoric probability subject to conserva- 7 q „z,™ _, r„»„^,„= TT ~, N 

F J J p = p(Tr\ sa C sa ,Tr\ asa C asa ;H,T;N,r]). 

tion laws. If 

Next, we define ensemble average of an oper- 

dP = dp[{x r }]p[{x r }} 

ator U as 

is the probabiHty of finding the system in op- (0) = J ^O, 

erator phase space volume element dp[{x r }], which is a map from Q M to Gm- 

then, This ensemble average has nice properties 

J dP = 1. (i) When O is constructed only from phase 

For a system in statistical equilibrium, phase space variables {x r } then this ensemble aver- 

space density distribution is constant age depends only on the constant parameters 

p[{x r }} = 0. X = {X,T,f]} 

Hence, p depends only upon conserved op- (O) ~ -^b Op- 
erators, conserved trace functional and con- (ii) Since the integration measure is unitary 
stant parameters. By going to a frame where invariant and O is made up of only {x r } using 
the system is not translating, accelerating or c . number coefficients, under a global unitary 
rotating the charges associated with Poincare transformation 
symmetry can be put to zero. In that case 



F (X) = UFoiU^XU)^ 1 



p = p(C, H,N). 

In addition the distribution function of dy- 
namical variables can depend on constant pa- 
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(iii) As a consequence Fo(X) depends explic- Clearly, -D e // is determined by a single real 
itly on A only and hence commutes with A number with dimension of action 



[A, (O)} = 0. 
Taking a specific case when O = C 
[A, (C)} = 0. 



(C)=i eff h. (119) 

The constant h will eventually be identified 
with Planck's constant. 



Being {anti-} self-adjoint A can be diagonal- 

ized by a unitary transformation, hence so The traceless ( C ) implies 

will be (C). Again specializing to anti-self- Tri e ff = 0. 

adjoint A and (C) we get ™ , , . , ,, . , 

The above mentioned properties of z e ff along- 

(i) for a real, non-negative, diagonal maq- . t , ,. , r ,. . r 

with property (iv) forces the dimension of 

nitude operator D eff and unitary diagonal , , . , „ 

11 Hubert space to be even and uniquely faxes 

p/jase-operator i e ff, . 

l e ff to 

(C) = i eff D eff . = . [d . ag ^ _ x x _ x ^ x 

(ii) Since, C is traceless, AT , , . , , r , . , 

JNext, we obtain the functional form of p 

Tr(i e ffD e ff) = 0. through maximizing the entropy defined as 

(iii) anti-self adjointness of C is ensured with S = J d/ip log p, (120) 

i e ff = ~i\ff subject to the following constraints 



J dpp = 1, J dppC = (C) 



and 

(iv) As a consequence of this decomposition J dppH = (H), J dppN = (N), 
i 2 = -I. 

e ^ which gives 

For an ensemble symmetric in Hilbert space 

basis, i.e. the ensemble does not prefer Pi ~ ex p(—Tr\C — 

any state, the averaged operator should have rH — r]N — Trj, 

identical entries as eigenvalues. Therefore, 



D eff = hi. 



Zj = J dp,exp(—Tr\C — 
TK-riN-J2Trj r x r ), (121) 

r 

where we have introduced a book-keeping 
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matrix source term j r for each matrix vari- 
able x r of same type (B/F) and adjointness 
that can be varied and set to zero. This 
helps us in obtaining the ensemble properties 
of functions made explicitly of the dynamic 
variables {x r }. 
In that case, 

{0)j = J d m O, (122) 

with 

(O)avj = J dfipjO. 
Using this distribution and partition function 
we can evaluate ensemble averages 



<A> = 
(H> = 
(N> = 



5 log Zj 

sx 

dlogZj 

dr 
d log Zj 

dr/ 



(123) 
(124) 
(125) 



and the mean square fluctuations 

A J, = (Trpt) 2 log Z 3 (126) 



5X 



(dr) 



A2 _d 2 \ogZ J 

N = 



(<9??) 2 



(127) 
(128) 



with P being any arbitrary fixed anti-self- 
adjoint operator. 

We further study the structure of averages of 
dynamical variables in the canonical ensem- 
ble. This study is essential for subsequently 
making connection with the emergent quan- 
tum theory. According to the previous dis- 



cussion, 

(C) = ieffD eff , 

and A is related to (C) only using c-number 
coefficients. Since D e ff is a constant times 
identity, 

A = \i eff , (129) 
for a real c-number A. Now for a unitary 
matrix U e ff that commutes with i e ff 

[U eff , A] = =► U eff XUl ff = A. (130) 

Clearly, the presence of TrXC term in the 
partition function breaks the global unitary 
invariance since under, 

q r -> U*q r U, p r -> U^p r U, C ->■ U*CU, 

TrXC = XTri eff C ->■ XTrUi eff U ] C. (131) 
Thus, presence of this term breaks the global 
unitary invariance to {U e ff}. With this resid- 
ual invariance in canonical ensemble we de- 
fine an integration measure as 

dji = d[U e ff]dfi, 

with d[U e ff] as the Haar measure on the 
residual symmetry group and dfi is the mea- 
sure over operator phase space when over- 
all global unitary transformation U e ff is kept 
fixed. 

In this case, for a polynomial P e // which is a 
function of i e ff and dynamical variables {x r }, 

f d [UeffWpReff 



ReffAV ~ fd[u eff w P 

Now if we fix the U e ff rotation 



(132) 



q r = Ul ff q r U eff , p r = Ul ff p r U eff , 



63 



it results in, unitary fixing does not disturb z e //, for any 

D rrt d tt operator O made up of {x r \ using c- number 

K eff = U eff K eff U eff ■ 

coefficients, 

and hence, 

jd[U eff ]WR effAv U eff [A,O AV ]=0. 



R, 



effAV 



M^//] When O = C, 

In the above equation, 

_ / dppR eff 



[\,c Av ] = o. 



R 



effAV J dfip 

Using properties (1l35h , the most general C 



where R eff and p are obtained from R eff and commuting with A is 

p by replacing q, p therein by q, p as defined ~ 1 . . . 1 , . . 



above. 

Writing ie/y = ia^lx, where 1^ is a unit 
K xK matrix, a general NxN matrix M can 
be decomposed in the form (with N = 2K) 



(136) 

Thus, the ensemble average in the unitary 
fixed system is different from the canoni- 



cal ensemble average (119). However, since 



M = M e ff + M\2 (133) the unitary fixing has been done by restrict- 

w j|.j 1 ing only one canonical pair, in systems in- 

1 i volving large number of canonical pairs, the 

M eff = ^(^o + <r 3 )M + + -((7 - a 3 )M_ 

z z restricted average should be close to unre- 

anc ^ stricted average. Therefore, A± should be 

M 12 = a x M x + a 2 M 2 . smaU 

Here cr^z = 1,2,3 are 2x2 Pauli matrices 

and do = 12- The are four K x K 

matrices. These new matrices satisfy 

[i e ff,M e ff] 0, {i e ff,Mi 2 } (134) Thus far we have progressed from the clas- 
2i e ffM e ff = {i e ff,M e ff}. (135) s i ca l theory of Trace Dynamics, to develop- 
ing a statistical thermodynamics of this the- 
ory. We now have in hand the tools necessary 
to describe the emergence of quantum the- 
ory in this thermodynamic approximation. 
The first step is the derivation of a crucial 
general Ward identity. This identity should 



6. General Ward Identity and Emergence of 
Quantum Theory 



In irreducible systems unitary fixing can be 
done by fixing the global unitary rotation 
of one canonical pair of dynamical variables. 
With the restricted measure we need to know 
the restricted canonical average C A y. Since 
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be thought of as an analog of the equipar- Evaluating term by term in (139) we have^J 
tition theorem in statistical mechanics, and 



5x 



its implications in the present context are Tr\5 Xs C = Tr 

deeply connected with the existence of the 
Adler-Millard charge and its canonical aver- 
age (119). Averages now are defined with for s N = i u rs Trx r 5x 

r 

respect to the restricted measure. Under a 
constant shift of any matrix variable apart 
from the restricted pair 

x r ^x r + 5x r ,r^R,R + l, with, further, for a polynomial W 

j d(i5 Xr ( Pj O) = 0. (137) S Xs W = W s Ll 5x s W s Rl , 



S Xa ii = u rs Trx r Sx s 



5 x Tr{C, i efJ }W = Tr({i eff , W}6 Xs C + 

{C,i eff }5 Xs W), 



Using this and the definition of average given 

abo e o O C,i e ff W we get rnial. W^ l ^ Rl is the left (right) fraction of a 



J djX5 Xa [exp (—TrXC - rH - r]N 
~J2 Tr 3rXr)Tr{C,i eff }W}=0. (138) 



where I labels each monomial in the polyno- 
mial. Ws l ^ m is the left (right) fraction of a 
monomial. Collecting above terms and plug- 



Using the chain rule we have 



ging back in (139) with some manipulations 
leads to the generalized Ward identity 

( A «e//)i = ((-TXueff + iV^uXueff ~ 

j dfiexp {-TrXC — tH — 77N X s u us j seff )TrCi eff W eff + [i eff W eff , x ueff ] 

YV ■ WTltr a "H +E^( w s R i{^ie f f}W s Ll ) ) j = 

2_^Trj r x r )[{-TrX5 Xs C - t6 Xb H ~ t V z Jeff 

r 

-r}5 Xs N -Trj s 5x s ){C,i eff }W with £ u = 1(-1) for fermionic q(p), zero for 

+5 Xs {C,i ef f}W} = 0. (139) bosonic x u and 

^2u us u rs = S ur . (140) 



7 u> = diag(riB, ^b, Qf, ^f), with 
for bosonic and fermionic d.o.f. 
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The Ward identity can be written more com- ensemble average 

pactly as //^ \ • t 

{Ceff)A V = leffh. 

(px ueff )j - u us j se}f {TrCi e}f W eff )j = With thege assumptions the RHS of the iden _ 

s 

(141) tity ( fi~42| ) simplifies to 



where, . ^ / SW 

l efA W eff^ueff\ ~ ^us I ~F~7~ 

*Dx ue ff = (—Tx ue ff + ir]t; u x ue ff) anc i (141) implies {Vx ue ff)o = 0. If we con- 



xTrCi e ffW e ff + [i e ffW e ff, x U eff] + sider W 7 = if in the Ward identity, we obtain 

£^ (w^-{C,i eff }wA (142) p Xue// = i eff [H efh x ueff \ - hx eff 

s,l V 7 e // 

which gives the effective Heisenberg equations 



From (141) we see that for a polynomial S 

of motion for the dynamics when sandwiched 

made up of x re fj and c-number coefficients, 

between SL{Xteff) and SR{x te ff) and aver- 
(S L{xtef f){T}S(x re f f))S n(x te f f))o = 0, aged over zero-source ensemble. 

for left and right decompositions of the poly- For an arbitrary polynomial function P eff 
nomial S. made U P of x reff 

We now make the following realistic assump- 

(SL(xteff)PeffS R (x t eff)}o = 

tions 

(i) Support properties of x ueff and C eff are ( s i>( x teff^effti ^Heff, P e ff}S R (x teff )) , 
such that 

suggesting that within our assumptions, H e ff 
~ TX ueffTrCi e ffW e ff j s a constant of motion. Next, for W = a v x v 

in rtl42b can be neglected. for some c-number parameter a v we get 



(ii) Chemical potential rj is very small, such i eff Vx ueff = [x ueff , d v x veff ] - t eff hw uv a v . 
that the term (143) 
iri£ u x ueff TrCi e f f W eff Thus > when multiplied with S L (x teff ) and 

Snixteff) and averaged over zero source en- 



in ( 142 ) can be neglected. In fact, for bosonic 

semble EFFECTIVE CANONICAL COM- 

degrees of freedom this term vanishes and it 

MUTATORS EMERGE 

is taken to be small for fermionic degrees of 

({Queff,qveff)) = ((Pueff,Pveff)) = 0, 

(iii) When the number of degrees of freedom 

{\Queff,Pveff)) = ieffhSuv, (144) 

is large, C can be replaced by its zero-source 
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with (( )) being anti-commutator (commu- field theory 

tator) for u,v being fermionic (bosonic). It af( n ^ c/rv 1N 

; y 1 S{{x reff }) ^ S{{X reff }), 



is important to emphasize that these commu- 
tation relations emerge only upon statistical 
averaging, as a consequence of there being 
the conserved Adler-Millard charge. At the 
level of the underlying theory of Trace Dy- 
namics, the commutators/anti-commutators 
amongst the above operators are arbitrary. 



with X re ff being quantized operators in 
quantum field theory. Here, i e jf acts as two 
blocks of i and —i. With the assumptions 

(a) in the continuum limit the trace La- 
grangian is Poincare invariant, 

(b) the Hamiltonian H e ff is bounded from 
below by the magnitude of the correspond- 
Next, let W = G be a self-adjoint operator ing e ff ect j ve three-momentum operator P e// , 

such that G generates canonical transforma- and there ig a unique eigenvector ^ Q with low . 

^ lon est eigenvalue of H e ff and zero eigenvalue of 

h 1 Vx ue ff = i e ffh 1 [G e ff,x e ff]—5x e ff Peff, 

(145) we also have a correspondence between trace 

Thus, on sandwiching between S L (x teff ) and dynamics canonical averages and Wightman 

Sn(xteff) that do not contain x U) and averag- functions in emergent quantum field theory, 

ing, we see that infinitesimal canonical trans- 4(S{{x re}} })) Av % = (vac\S({X reff })\vac). 
formations at the ensemble level and within 
the above-mentioned assumptions are gener- 
ated by unitary transformations 



Using (143), (144) and the proposed cor- 
respondence we have at the quantum level 

77 ^ - pvn ( 7 - ,A [X ueff ,a v X veff } = i eff hu uv a v , (146) 

u can eff — yi e ffll ^jr e ff). 

„, r , , , , ,. , £ which gives the appropriate commutators at 

lherefore, we have at hand the essential fea- 

r , £11,1 mi r the quantum level, both for bosonic and 

tures of quantum field theory, lhe {anti- 

J . , , . i x- • fermionic degrees, 

jcommutator structure, time evolution m 

Heisenberg picture and unitary generation of Time evolution is given by 

canonical transformations emerge when we X ueff = i eff h~ 1 [H effl X ueff ], (147) 



carry out the statistical thermodynamics of 
the matrix variables. 



or equivalently for a polynomial S e ff of 

i X *ff} 

Now we make the following correspondences 

between operator polynomials in trace dy- S U eff — ieffh [H e ff,S ue ff]. (148) 

namics and operator polynomials in quantum For the fermionic anti-commutator to be an 
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operator equation, the following adjointness so that, 



assignment is required, 



Aejf - qreff *> Q reff - d^ f{t) = leffh -± ulff{t)Heff . (154) 

4eff = Preff <"> Preff = i e fM ef p(U9) 

From the Heisenberg picture time- 
such that independent state vector ip and time depen- 
{$ ue ff,^l eff } = H8 UV . (150) dent operator S eff (t), in the Schrodinger 

For the bosonic sector we obtain creation and 



picture we perform the construction 



annihilation operators A re tf, Al„ ff such that, , , . TT , . , 

U i>Schr{t) = U ef f(t)i> 

Queff = -^(A reff + A\ eff ) s eff schr = U eff (t)S eff (t)ul ff (t), (155) 
V 2 



P « e // = -ft. — ( A reff ~ A\ eff ) (151) which gives, 

V&eff 

and, i e ffft-^:ipschr{t) = H eff ijj Schr (t) 



eff 



d_ 

dt 



[Aueff, A veff ] = [Al ff , Al ff ] = 0, j t S eff s ch r = 0. (156) 



[Aue ff ,Al eff } = h6 uv . (152) 
Thus, we have the correct commutation/anti- 
commutation rules for bosonic/fermionic de- 
grees of freedom in both the i e ff = ±z sec- 
tors. 



For obtaining the Schrodinger equation, we 
make contact with space-time by taking the 
label r as x. In that case, the fermionic anti- 
commutator becomes 



{VffWMffW} = h5\x-y). (157) 
We have assumed that H e ff is bounded from 
Once we have the Heisenberg equations below havin S the vacuum state |V> as the 
of motion p7fc and p8l), we can make the lowest eigenvalue state and that * e// should 



transition to the Schrodinger picture as usual, annihilate it, 

without making any reference to the back- \P e y^|V) = 0. 

ground trace dynamics theory When the ef- _, „ 

1 heretore, 

fective Hamiltonian has no time dependence, 

we define (V\^ eff (x)^ eff (y)\V) = M 3 (f - y) (158) 

U eff (t) = exp (-iefffrHHeff), (153) 
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Similarly, for the bosonic operator A e ff 
(V\A eff (x)Al ff (y)\V) = h5 3 (x-y). (159) 



Thus, an analysis of (158) will analogously 



apply for (159) as well. Defining, 



m^ n {x) = (V\y eff {x)\n), 



we obtain 



and J d 3 y^* n (y)^ t 



(160) 



from ( 158 ). Then, using the Heisenberg equa- 



tion of motion, 



= (V\i e „h- 1 [H e „,y ell (x)}\n). (161) 
Again, defining 

we get 

[Heff^eff{x)} = -n ef f(x)^ eff (x), 



thus modifying (161) into a Schrodinger 
equation 



ieffh^nix) = U e f f{x)^> n {x) 



(162) 



However, in deriving (162) we have made 



certain approximations valid at equilibrium. 
More explicitly, we replaced C by its canon- 
ical average. If we also consider the fluctu- 
ations about the average quantities we have 
possibilities of obtaining a stochastic equa- 



tion of evolution by adding stochastic non- 
linear terms to the Schrodinger equation. 
Herein perhaps lies the greatest virtue of 
Trace Dynamics. By treating quantum the- 
ory as a thermodynamic approximation to a 
statistical mechanics, the theory opens the 
door for the ever-present statistical fluctua- 
tions to play the desired role of the non-linear 
stochasticity which impacts on the measure- 
ment problem. 

If we consider fluctuations in C to be de- 
scribed by 

AC ~ C - i eff h 

\{C,ieff} = -h+^{AC,i eff } 

l -{AC 1 i eff } = -h{K+N), (163) 

with fluctuating c-number /C and fluctuat- 
ing matrix A/", we obtain the modified lin- 
ear Schrodinger equation (restricting to the 
i e ff = i sector) 



$) = [ih- L {-l + (Ko(t)+i1Cx(t))}H t 



eff 



+ -i(M (t)+iM 1 (t))]\$), (164) 



where 0, 1 label the real and imaginary parts 
of /C and 

M(t) = J]fMf(^. 

r,l 

In the above equations m r is the rest mass of 

the r-th species such that the Hamiltonian is 
i 

H = ) J ) J -im r [^l u ipri] + const (165) 

r I 
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with I labeling a general complete basis set dom, much like in the case of Brownian mo- 
and we have used the correspondence pro- tion fluctuations. A Brownian motion is de- 
posed between trace-dynamics and quantum scribed by a stochastic process dW™ satisfy- 
mechanics, ing the following Ito table 

l*> = n»^/IV>- W ») 2 = 7„*, (166) 

As we will see next, with the added as- dW^dW™ = 0,m^n, (167) 

sumption of norm conservation through equa- dW™dt = dt 2 = (168) 



tion (164) Trace Dynamics connects with 

a CSL type non-linear stochastic equation. For our case we make the following identifi- 

The weak link however, in the chain, at the cations: 

present stage of our understanding, is the as- (i) For c- number fluctuations /C ,i 

sumption that norm is conserved. Rather 

, . ,. ... , fl ih^Kodt = ifinWi , (169) 

than being an assumption, this should fol- 
low from the underlying theory, and hope- ^ ^ l ~ P R( ^t > (1^0) 
fully with improved understanding this will 

with the following Ito table 

become possible in the future. On the other 

hand, the presence of anti-Hermitean mod- ( dW t R ) 2 = ( dw lf = dt,dW R dWl = 0. 
ifications in the Schrodinger equation is in- (ii) For the fluctuating matrix Mo,i having 
evitable from the Trace Dynamics viewpoint, spatially correlated noise structure 
since it is possible for C to have a self-adjoint 

part as well. ? M <& = % J d'xdW? (x)Ml(x), (171) 

~Mi = J d 3 xdW t R (x)M?(x), (172) 



7. Stochastic modification of the Schrodinger 
equation 



with, 



dWf{x)dWt(y) = jdtS 3 (x - y) (173) 
dW t R (x)dW t R (y) = 7dM 3 (f - y) (174) 
dW t I (Z)dW*(y)= 0(175) 
dWfdWKx) = dW*dW R {x) =0(176) 



We further motivate that the fluctuations 
can be described by linear superposition of 
white noise terms owing to the hierarchy be- 
tween the length scale associated with the C 
fluctuation and the length scale characteriz- dWfdW^x) = dW R dW R (x) = 0. (177) 

ing the emergent quantum degrees of free- 
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These identifications turn Eqn. (164) into a lapse models. Of course, at the present stage 
stochastic differential equation of understanding, it cannot pick out one col- 

lapse model out of the many discussed, nor 
— [ ih H e ffdt + i[3idW t ti e ff provide a theoretical origin for the values of 
+(3 R dW R H eff + i J d*xdWl{x)M[{x) the CSL parame ters A and r c . Nonetheless, 

+ / d 3 xdW R (x)Ai R (x)]\$>) ° ne canno ^ esca P e ^e profound and natural 

hypothesis that on one hand thermodynamic 
The above evolution is not norm preserv- equilibrium corresponds to quantum theory, 
ing. The idea is to define a physical state and on the other hand fluctuations around 
l^f) ^= -^j^yj with a conserved norm (as as- equilibrium correspond to stochastic modifi- 
sumption) which, along with the criterion of cations of quantum theory. Why the effect 
absence of superluminal signaling, after some of stochasticity must be larger for larger sys- 
calculation at length, gives terns remains to be understood. Nor is it un- 

derstood why norm should be preserved dur- 
d\q?) = [-ih- l H eff dt + ifaHeffdWt 1 ing evolution in Trace Dynamics: one should 

~2^^ff ~^ pRffleff — (H e ff)) 2 }dt no t have to put this in as an assumption into 
+f3 R (H e ff — (H e ff)) 2 dW R the theory, but rather have it come out of the 
+i J d z xM I {x)dWl{x) - ^dtx underlying theory as a consequence. 



d 3 x[M T (x) 2 + {M R {x) - (M R {x))) 



2l 



d A x{M R {x) - {M R {x))fdW R {x) 



i*>- 



For explicit demonstration of the collapse 
of the wave-function induced by stochasticity, 
we study a simplified version 
This equation is a stochastic non-linear 
Schrodinger equation which has the martin- 

gale structural of spontaneous collapse mod- + P]A 2 ])\^)dt + f3 R (A - (A))\V)dW t R 

els, and is capable of explaining state vec- + i(3jA\^)dW^ , (178) 

tor reduction. In this sense, Trace Dynamics 
is an underlying theory for spontaneous col- wr ^ n 

~ , . 1T , dp = ih~ l \p,H e ff\dt \$\ 2 \A, \A,p]dt 

8 For stochastic |*), H Lr ' //J 2 

+P R \p, [p, A}}dW R + ifii [A, p\dWl (179) 



dp=(d\*)){*\ + \*}(d{y\) + (d\*)){d{*\) 



Defining E[ ] as expectation with respect 
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to the stochastic process, (E[dWf-) = = Also, we obtain from ( |181 ) 

E\dWfX) and variance of A, cn T /mi ^ M^] 

E[V{t)] - l + 4/3|V[0]t' 

u ^ ^ " / ' anc l hence a time scale of reduction as T 

and using the Ito product rules gives 4/3 2 ? V r [0]. 

d£[V] = £[aV] = -A/3 2 R E[V 2 ]dt. (180) 

Therefore, 



We can also see that in such a reduction 
scheme the Born probability rule follows for 



E[V(t)] = E[V(0)] -4/3 2 f dsE[V(s) 2 }. the outcomes. To see that, let us take IT a 

Jo 

(181) as the projector into the a-th eigenstate of 

Using the inequality operator A, 

< E[(V - E[V}) 2 } = E[V 2 ] - E[V} 2 , n « = l a )( fl l- 

Now for any operator G commuting with 



this becomes 
E[V{t)]<E[Vm-^l I dsE[V{s)} 2 



H e ff and A 



o 



If the initial state of the system is 



Non-negativity of variance suggests that E W G)]=Tr(-ih^[G,H. lA E\p] 

E[V(oo)] = and again as V(t) is not sup- ~ ^W\G , A][A, E[p\})dt = 0.(182) 
posed to be negative anywhere this will en- 
force 

a 

As the variance in expectation of A goes to at t = when the stochastic evolution has 
zero asymptotically, the system in this way no ^. s ^ ar ^ ec i 
results in one of the eigenstateg^] of A. The 
demonstration of collapse using a system- 
apparatus interaction in the QMUPL model As we have argued, when the evolution is 
in Sec. II is a specific explicit application of driven b y A the s y stem results in a Particu- 
this general analysis. lar eigenstate |/) with some probability P f . 

Then, for the a-th eigenstate, 

E[(U a )] t=OQ = J2(fK\f)Pf = Pa- (184) 



E[(U a )} t=0 = (U a ) t=0 = \ Pia \ 2 . (183) 



9 An open question is to make sure that the preferred f 

basis for the collapse, as chosen by Trace Dynam- AT . A , , , . 

J J Mow, since A was taken to be commuting 

ics, corresponds to some sort of position basis, in 

order to guarantee that macroscopic objects are with H eff , we can choose their simultaneous 
always localized in space. eigenstates, which we can call \a). Therefore, 
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operators IT a constructed from these eigen- 
states will commute with H e ff and A, result- 
ing in the time- independence of £?[(n o )] as 



evident from (182). Therefore, 

E[{U a )] t=0 = (n a ) i=0 = E[(U a )] t=O0 (185) 

giving P a = \pia\ 2 - Thus, we have obtained 
the Born probability rule. 

We have seen that when treated fluc- 
tuation around the thermodynamic limit of 
trace dynamics theory, the emergent non- 
linear equation captures the essential features 
of CSL, and in a sense, can possibly be a the- 
oretical motivation for the phenomenological 
CSL equation of evolution. Of course at the 
present stage of understanding, trace dynam- 
ics makes no definite prediction for the actual 
numerical values of the CSL parameters, and 
this remains a challenge for the theory. 

B. Gravity induced collapse 

The general theory of relativity dictates 
that gravity is the curvature of space-time. 
This curvature is produced by classical ma- 
terial bodies. However, even the motion of 
classical bodies possesses intrinsic quantum 
fluctuations, and these fluctuations imprint 
a small uncertainty on spacetime structure. 
When one considers the motion of quan- 
tum mechanical objects on such a fluctuat- 
ing spacetime, the coherence of the quan- 



tum state can be lost, providing a possi- 
ble mechanism for wave-function collapse in 
the macroscopic domain, while leaving mi- 
crophysics untouched by gravity. Counterin- 
tuitive though it may seem, gravity possibly 
plays a profound role in bringing about wave- 
vector reduction, as the studies described be- 
low indicate. 



1. The model of Karolyhazy [K-model] 



The proposal of Karolyhazy: (Karoly 



hazy 




1966 





with a smearing of space-time which results 
from the fundamental uncertainty in quan- 
tum theory being forced upon space-time 
structure. It starts with a viewpoint that 
nature 'somehow' tries to reconcile classical 
general relativity with quantum mechanics 
as much as possible. Space-time has in gen- 
eral, a fairly definite metric structure mainly 
determined by classical massive objects with 
fairly definite positions. However, the met- 
ric should not be completely sharp, and must 
have an in-built haziness to avoid contradic- 
tion with the fundamental quantum aspect of 
massive objects (spread in position and mo- 
menta). Even a macroscopic massive body 
will have to satisfy 

h 



5x 5v > 



2m 
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where m is the mass of the body. The result- The uncertainties Ax and Axo are the un- 

ing haziness in the metric produced by the certainties in the top and bottom of the seg- 

body leads to a stochastic correction in the ments as well as in the length of the seg- 

evolution of state- vectors in quantum theory, ments. Minimum amount of uncertainty in 

The basic idea of the approach is that the len S th of segment will be introduced if 



when a wave packet of the center of mass of 
a body, sufficiently narrow in the beginning, 



we choose 

Ax = Ax . (187) 



spreads out in the Schrddinger evolution, into Clearl y the uncertainty in the length of seg- 
a space domain larger than a critical value ment decreases with increasing M and point- 
characteristic to the system), the coherence like description becomes progressively more 
between distinct parts of the wave function valid - Now ; the gravitational radius of the 
gets destroyed, owing to space-time haziness. mass M is bounded by the fact that it should 
This is interpreted as a signal for stochastic not be g reater than the spread Ax, 
reduction of the extended wave function to Ax ~ GM/c 2 . (188) 
one of its smaller, coherent parts. 



From Eqns. (186), (187) and (188), the 



Quantum imprecision of space-time uncertainty in the length of the segment is 

Let us consider a world-line segment s = ^' ^ 

cT, in a flat space-time. We wish to esti- , . x2 /A n2 ^ m 

F (As) 2 = (Ax) 2 = 7jjj- cT 

mate the precision with which we can realize ^ 



s. 



this segment. Thus, the segment of t-axis is 2Mc 2Asc 3 

to be realized by the narrowest possible tube , . > 2 f Gh\ 2 ^ 2/ o /<r ,„ N 

(As) 2 = — s 2/ \ (189) 



formed by a standing wave packet. Let, at 



V2c 3 J 



the start (i.e. at the bottom of the world- This relation, giving the minimum amount 
line segment) the width of the wave packet of uncertainty in space-time structure is of- 
be Ax . For mass M of the wave packet the ten known as the Karolyhazy uncertainty re- 
velocity spread is lation. Therefore, we should be careful in us- 

/\y = ^ ing classical space-time considerations once 

2MAx ' , , , r , 

the length of the segment starts approaching 



The corresponding spread at the end (i.e. at 
the top of the line segment) will be 



its uncertainty value; thus providing a critical 
length scale for the system. 



h 

Ax = AVT = — — cT. (186) AT i , . , 

2MAxoc Mext we consider a physical space-time 
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domain of nearly Minkowski metrics with a 
corresponding smear structure as argued in 



(Karolyhazy, 1966 Karolyhazy et al, 1986). 



We introduce a family {g^ u } of matter-free 
metrics very close to the Minkowski met- 
ric, where different /3 mark different mem- 
bers (hence different metrics) of the family. 
The proper length s = cT between two world 
points Xi and x<i will be defined as the mean 
value of the lengths sp corresponding to dif- 
ferent 

s = s p , (190) 
with the bar describing average over 0. The 
uncertainty in the line segment is defined as 

As = [( S -^)2]V2 ) ( i9i) 

In the family of metrics /3 = gives the 
Minkowski metric. In the present analysis at- 
tention will be confined to the case in which 
we do not have macroscopic bodies moving 
relatively to each other with a velocity near 
that of light. The co-ordinate system will 
therefore be assumed to be one relative to 
which all macroscopic bodies move slowly. 
This will enable us to confine our use of 
the set {sp} to non-relativistic many-particle 
wave equations in spite of the fact that by 
invoking curved manifolds we are employing 
the language of general relativity. 

Since we are considering only slowly mov- 
ing particles, t) C c, only (goo)p part of the 
metric will be required for the analysis. The 



general form of the metric in the family is of 
the form 

(goo)p(x) = -l+ 1 p(x); (/3^0). 

Since the space-time is matter-free apart 
from the test particle, we have 

□7/3 = 0. 

Now, the idea is to fix the set jp in such a 
way that the length of the world-line 

dx» dx vl 1/2 



sp 



dt 



9. 



^ dt dt 



(192) 



is averaged to (190), and the uncertainty ob- 



tained from (191) is the same as obtained in 



( 189 ) . We thus do not regard the functions 7^ 
as dynamical variables, rather we represent 
physical space-time by the whole set {g^ u } at 
once. In this spirit we construct 7s through 
their Fourier series, 



7/3 (x) 



-T 

3/2 



cp(k) exp [i(k ■ x — cut)] + c.c) , (193) 



where L is the length of an arbitrarily chosen 
large box (for normalization), 

k = —n and u = c\k\. 

LJ 

We now choose an integer iVg > 2 for each k 
and introduce a random variable ot(k), such 
that 

2tt 

a(k)e— [0,l,2,....,JVg-l]. 

k 

For a particular a{k) a particular Fourier co- 
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efficient Cp(k) is given as and a corresponding phase operator in the 

,t\ „ /lN r . /r \, /,^.n wave-function describing the quantum state. 

c^(fc) = /(fc)exp [za(fc)]. (194) 

The final results on wave- vector reduction are 

The unknown function f (k) is obtained from , , , , , . 

the same as those described below. 

the scheme proposed above and is found out 

, For considering wave propagation 
to be 

1 /3 

f(h\ — ( 7,-5/6 ac\K\ (Schrodinger type evolution) in this 'hazy' 

, ./ ,. r € ^ space-time, we introduce a family {i/jr} of 



using (189). The contribution to fk for large 

, £ , r . ,i wave-functions corresponding to the metric 
values of k comes from the requirement that 

t-i /FTnr& i ill tj . £ family {q^„\. For a single scalar elementary 

Fqn. (189) should be valid even if s is very ujjwj 

ii ^i i FTTtrX ■ c i • „ particle, via the relativistic Klein-Gordon 
small. Ulearly, (189) is not meaningful m the 

limit s — > and a cut-off is assumed: f(k) = e 1 ua ^ lon 

for k > 10 13 cm" 1 , s < l^cm. It is as- ~J=^W^^ -fad ~ ("jf J $ = °> 

serted that details of the cut-off are not im- we obtgdn the non . relativistic generalization 

portant, and only long-wave components are f) ( h 2 \ 

lh dt^ = [~^ 2 + V ^^- (196) 



relevant. This has been contested in (Diosi 



and Lukacs[ [l993| ) where it has been claimed The sma11 perturbation Vp is given as 

that this cut-off is at a very high physically Vp(x,t) = - . 

2 



unacceptable value of k and leads to absurd 
situations such as neutron star scale densities 
all over space. However, it seems that this 
objection can possibly be avoided by work- 

miC 2 -fp(xi,t) 



More interesting is the case for many parti- 
cles: the multi-particle equation, where Vp is 
replaced by 



ing entirely in real space, without going to Up({X},t) = 

Fourier space (Frenkel, 2002). The analysis 



2 

{X} = {x l }. (197) 



of Karolyhazy has been repeated by Frenkel, 

according to whom some of the Fourier sums rp , u v , , . a ± t u 

to lo see the qualitative effect of such a 

diverge in some intermediate expressions, but • , , • ... , £ 

to ^ ' smearing we start with an initial compos- 

"in the formulae for physical quantities these ., t , lT( rfv\ n\r iwu + 

v J M ite wave-function W ({A }, 0) for all the met- 

sums are convergent". In this work, the im- . r m/-. Af , , ,. A . a . 

to ' rics {gp }. After the evolution different 

pact of the Karolyhazy uncertainty relation is y p ({X},t) will become different. We can 

realized, not by introducing a family of met- , j 

' J write to a good approximation, 

rics, but by introducing a local time operator, 

*p({X},t)n* ({X},t)e^W\ (198) 
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with {V^} no longer behave as a single coher- 

, rrv-i .\ 1 [* J./TT i r +\ ent wave function. When the original coher- 

<MW^) = -W dtUf,({X},t). (199) 6 



ent set develops incoherent parts of compara- 
ble weights, it is taken as signal for stochas- 
tic reduction of {ipp} to a single cell. There- 

these two points in configuration space for r , . . , . , 

fore, this stochastic reduction scheme is gov- 

different 8. The answer will depend on 8 , , , , , . , . n 

erned by bchrodme;er evolution and stochas- 



Let us choose and fix an X\ and an X 2 , 
and calculate the difference in phase between 



and on time. The root mean square spread 



tic part comes through smearing of space- 
in the phase (average is over 8) can be es- , . , . m ,. 

v ' time metric, this process provides us with 



timated clS db function of {Xi,X 2 } and time 
t. The uncertainty in the relative phase de- 



a description of a physical phenomenon tak- 
ing place regardless of the presence of any 

pends only on the separation between the two , c j- 

observer, btill this formalism indicates to- 

points in configuration space, and for a suffi- , , . , c , 

wards but does not provide any formal em- 

ciently large separation can reach the value , , r . n c , , ,. . 

bedding of idea of stochastic jumps into evo- 



7T. 



lution in a consistent mathematical frame- 



Microscopic and Macroscopic behavior: work - lt can be heuristically argued that 

microscopic quantum particles will take as- 



For a single quantum particle of mass M 
for small values of a = \xi — x 2 \ the spread in 



the phase 

A(a) < vr, 
and only for a large critical value a c 

A(a c ) ~ 7r For an elementary particle of mass m, it 



tronomically large time before their wave- 
functions "spill over" a single coherence shell, 
making the possibility of stochastic reduction 
very remote. 



will be achieved. The spread in the phase can be shown that i MB ' 



ft 2 1 / L \ 2 ft 
and the separation for which the critical value flc ~ ~ J ) L; L « , (200) 

G m 3 \LpJ vac 



and the critical time of reduction can be 
shown to be 



is reached can be calculated, as described 
above. We next discretize the space in terms 

' coherence cells' of dimension a r . If initially 

r c ^— ±. (201) 
the particle is confined to a single cell, then 

For a proton one finds 

Schrodinger evolution will try to spread the 

1 x Tx- C Q>r ~ 10 ' 10 M s 

wave packet, resulting m the wave function 

extending over to different cells and the set thus showing that one can never observe 
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ma 2 . 



wave-packet reduction for a proton. The The reduction time is again given by (201). 
origin of the expression for the reduction For a ball of R — 1 cm and for terres- 
time lies in the fact that according to the trial densities this gives a c pa 1CT 16 cm and 
Schrodinger equation, a wave-packet initially r c pa 1CT 4 cm. The wave-function under- 
spread over a c will spread to a size 2a c over goes 10 4 expansion-reduction cycles per sec- 
time r c . When this happens, we could take ond, and at the end of each cycle the mo 
that as an indicator of loss of coherence, mentum performs a jump Ap c of the order 
and hence stochastic reduction. The dynam- h/a c which corresponds to a velocity shift of 
ics thus consists of cycles of deterministic the order a c /r c ~ 1CT 12 cm/sec. These re- 
Schrodinger evolution followed by stochastic peated kicks amount to an anomalous Brown- 
jumps - something fully reminiscent of the ian motion and a tiny associated energy non- 
GRW model; a comparison to which we will conservation of the order h 2 /Ma 2 , another 
return shortly. In fact, r c is analogous to A -1 feature in common with spontaneous collapse 
in GRW, and a c is analogous to r c . models. 

One can try to understand the transition 



For more complex systems, such as a 
macroscopic body, one works with the cen- 
ter of mass co-ordinate. However, care is 
needed because the gravitational perturba- 
tion described by the multi-particle potential 



region from micro- to macro- behavior. We 
have seen that for R pa lcm we have a c <C R. 



(197) depends on extended region of space. 



Furthermore, the expressions (202) and (200) 
become the same when a c = R. So one can 
now classify 

(i) a c > R (i.e. h 2 /G > M 3 R) micro- 
behavior regime 

(ii) a c pa R (i.e. h 2 /G pa M 3 R) transition 
region 

(iii) a c < R (i.e. h 2 /G < M 3 R) macro- 
behavior regime. 



Still it can be shown that only the phase of 
the wave-function of the center of mass is af- 
fected, and the already introduced concepts 
of coherence cells and coherence length a c can 
be applied to the center of mass coordinate. 
In such cases not only mass M of the system 
but the size R also enters in the expression 
for a c (as arbitrariness in metric will be expe- 
rienced throughout the size). It can be shown 
that 

-^2\ !/ 3 j^2/3 / jj> \ 2 / 3 ^ Setting a c = R in (202) and assuming den 

- 1 L ', L 



If a c ^> R it can be shown that (Frenkel 



2002) the expression (200) continues to hold, 



for a micro-object having an extended linear 
size R. 



GJ M \LpJ mc sity to be about 1 gram/cc for terrestrial bod- 
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ies gives for the transition region 



a tr « 10- 5 cm, r tr « 10 3 s, M* r ss 10 



eToO [1982b. 



-14, 



(202) 

It is significant that a tr coincides with the 
favored value for rc in the GRW and CSL 
model. The transition mass corresponds to 
about 10 10 amu. Note that because Planck 
length and the size of the body also enter the 
picture, the transition occurs at a mass much 
lower than the simplistic but much higher 
Planck mass (10~ 5 grams). 

Interestingly, the jump velocity Av c = 
a c /r c in a reduction cycle takes its maximal 
value in the transition region a c ~ R and 
decreases on either side away from this tran- 



sition region (Karolyhazy et al, 1986) 



Thus we obtain a transition point which, 
in principle can be tested upon. The mea- 
surement process can be argued as inter- 
actions resulting in significant change in 
mass distribution of the whole set up (sys- 
tem+surrounding) making the definite and 
different outcome states incoherent thereby 
reducing the state of set-up to a particu- 
lar outcome. Still the formal mathematical 
framework of this idea is missing making it 
a challenge to precisely calculate the charac- 
teristics of quantum state reduction. 

The K-model has also been discussed 



Comparison with GRW model 

There is a fascinating similarity between 
the K-model and GRW model, despite signif- 
icant differences in detail. The overall pic- 
ture of Schrodinger evolution interrupted by 
stochastic reduction is the same. In the K- 
model, the origin of stochasticity lies in the 
intrinsic uncertainty of space-time structure, 
whereas in GRW the origin is left unspeci- 
fied. However, both models have length/time 
scales (a c /r c in K-model, r c /\~ l in GRW). 
There are no free parameters in the K-model, 
whereas GRW introduce new parameters A 
and rc- Thus it is entirely possible that grav- 
ity might provide the fundamental underpin- 
ning for models of spontaneous collapse. Of 
course a mathematically rigorous treatment 
of gravity in the K-model remains to be de- 
veloped, but the physical principles and semi- 
rigorous results already obtained are highly 
suggestive by themselves. 

An important early study comparing the 



K-model and GRW was made in (Frenkel 



1990). It should be noted that while in both 



in (Frenkel 1977 1990, 1995, 1997, 2002 



Karolyhazy 1974 1990, 1995 Karolyhazy 



cases the reduction time decreases with in- 
creasing mass, the quantitative dependence 
is different. In the K-model, for a c > R the 
reduction time falls as 1/m 5 , and if a c <C R 
it falls as l/m 5//9 assuming a fixed density. 
In GRW, the reduction time simply falls as 
1/m, whereas we have seen that in the CSL 
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model the dependence is more complex. Sim- 
ilarly, a c falls with increasing mass. While r c 
in GRW does not depend on mass, the linear 
size to which the stochastic reduction con- 
fines an expanding wave-packet does depend 
on mass - this linear size is the analog of the 



coherence cell a c (Frenkel, 1990). In the light 



of modern experiments there is perhaps need 
for a more careful comparison between the 
quantitative predictions of the K-model and 
GRW/CSL. Also, a careful quantitative de- 
scription of the quantum measurement pro- 
cess, showing the emergence of the Born rule, 
seems to not yet have been developed in the 
K-model. A time-evolution equation for the 
density operator in the K-model, analogous 
to the corresponding equation in GRW, has 



been discussed in (Frenkel, 1990). See also 



turn probe over a time T, then its average 
g(r, t) over a volume V exhibits an uncer- 
tainty which is universally bounded by 

(5g) 2 > hG/VT (203) 

This is Diosi's analog of the Karolyhazy un- 
certainty relation, and the idea now is to 
see how how this intrinsic quantum impre- 
cision in the space-time metric affects the 
Schrodinger evolution of a quantum state in 
quantum mechanics. 

To this effect, Diosi introduces the con- 
cept of a universal gravitational white noise, 
by proposing that the gravitational field pos- 
sesses universal fluctuations [in other words 
the potential 0(r, t) is a stochastic variable] 
whose stochastic average equals, up to nu- 
merical factors of order unity, the intrinsic 



related discussions in (Unturbe and Sanchez- 



Gomez, 1992). 



2. The model of Diosi 



Diosi's approach (Diosi, 1987), while be 



ing similar to Karolyhazy's, is inspired by 
the famous work of Bohr and Rosenfeld 



(Wheeler and Zurek, 1983) which investi- 



gated the principles of measuring the electro- 
magnetic field by apparatuses obeying quan- 



tum mechanics. It was argued in (Diosi and 



Lukacs, 1987) that if a Newtonian gravita- 
tional field g = — V0 is measured by a quan- 



uncertainty given by (203) 

([V0(r, t)} 2 ) - [(V0(r, t))} 2 = const x hG/VT 

(204) 

From here, it can be shown that, assum- 
ing ((f)(r,t)) = 0, the correlation function of 
(f)(r, t) is given by 

(0(r,i)0(r',t')> = frGlr-iT 1 S(t-t') (205) 

The probability distribution of the stochastic 
variable 0(r, t) is completely specified by this 
correlation function if the distribution is as- 
sumed to be gaussian [gaussian white noise]. 

Next, one asks for the effect of the stochas- 
tic fluctuations in on the propagation of 
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the quantum state tp of a system whose evo- the master equation can be written as 
lution is assumed to be described by the ^ 
Schrodinger equation (X\p\X ) — ^("^1 Ho, p(t) 



ihij = (^H + J (f>f(r) d 3 r^j ^(t) (206) 



X') 

-[T d {X,X>)\-\X\p{t)\X>) (209) 



where /(r) stands for the operator of the local Just like in decoherence and in models of 

mass density of the system. spontaneous collapse, the second term on the 

, . , . , , right hand side destroys interference between 

ip is now a stochastic variable, and the cor- 

.. . , , , + the states \X) and \X') over the characteris- 

respondmg density operator p = (ip(t)ip J (£)) 

tic time r^, and this effect can become signif- 
icant if the difference between the mass dis- 
tributions f(r\X) and f(r\X') is significant. 

To estimate the scale of the gravitationally 
induced violation Diosi considers a dynami- 
cal system consisting of a rigid spherical ball 
of homogeneously distributed mass m and ra- 
dius R, so that the configuration X is repre- 



obeys the following deterministic master 
equation for the assumed gaussian white 
noise 



p= -~h 

G f f d 3 r d 3 r' 
2h 




H ,p(t) 



/(r), f(v>),p(t) 



(207) 



The second term on the right hand side is the sented b ^ the center of mass coordinate x. 

damping term which represents the universal Tne characteristic damping time r d is shown 

violation of quantum mechanics. ^° be 

To compute the nature of violation, de- ^(x,x') = h[U(\x - x'|) - tZ(O)]- 1 (210) 

note the configuration coordinates of a dy- where jj is the grav it a tional potential be- 

namical system by X, and denote the corre- tween two sp h e res, each of mass m and radius 

sponding mass density at a point r by f(r\X). R The maste r equation can now be written 

Given a pair of configurations, a characteris- as 
tic damping time Td(X,X') is defined by 

i( x | p | x ') = ^( V 2 -V /2 )(x|p|x'> 
G r f « . dV m 1 2m v A m 1 



[f(r 


X) - f(r 


X')][f(r'\X) - f(r'\X')} 




r — r'| 





^[[/(|x-x'|)-[/(0)]<x|p|x'> (211) 



(208) 



We define the coherent width I of a given 
state as the characteristic distance / = |x— x'| 
Introducing the coordinate eigenstates \X) above which the off-diagonal terms (x|p|x') 
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become negligibly small. The time-scale tun but the exact mathematical treatment re- 
over which kinetic changes are introduced mains to be found. Importantly, the tran- 
due to ordinary quantum evolution given by sition l crit = R happens at the same value 
the first term on the right hand side is of the l cr u = ti 2 /Gm 3 in both the models. Notice 
order ml 2 /h. A crititical length l CT u is defined though that for small masses l cr u is not in- 
by equating tki n (lcru) an d the damping time dependent of R, unlike in the K-model. For 
T~d(lcrit) a proton, taking R to be the classical radius 



mllJh = h[U(l crit ) - UiO)}- 1 (212) 

If the coherent width I of the quantum state 
is much smaller than the critical value Z cr # 
then the standard quantum kinetics domi- 
nates and damping is not effective. On the 
other hand if I 3> l cr it then the coherence 
of the state will be destroyed by the grav- 
itational damping term in the master equa- 
tion. l crit is the analog of the phase coherence 
length a c of the K-model, and the length pa- 
rameter rc of GRW. Also, there clearly are 
analogs of r^lcrit) in the other two models. 

One can now show that in two limiting 
cases l crit takes the following form: 



10 -13 cm, Diosi estimates l cr u to be 10 6 cm, 
which is curiously much smaller than the pre- 
diction 10 25 cm for the K-model. Also, the 
reduction time is 10 15 sec, much smaller than 
in the K-model. However, the models are in 
better agreement in the macro- region, and in 
Diosi's model too, the transition parameters 



are the same as that given by Eqn. (202). 

Subsequently, Diosi took the inevitable 
step of casting the master equation in 
the equivalent language of a stochastic 



Schrodinger equation (Diosi, 1989). He 
called this model QMUDL (Quantum me- 
chanics with universal density localization). 
It is similar to his QMUPL model, which we 
reviewed earlier in this article, except that 



hrit ~ (h 2 /Gm 3 ) 1/4 R 3/4 , if Rrn 3 > h 2 /G, the localization is not in the position operator 

~ (h 2 /Gm 3 ) 1/2 R 1/2 , if Rm 3 < h 2 IG g ' but in the mass densit y operator /(r) intro- 
duced above. The universal free parameter A 

(213) 

of QMUPL is now replaced by the gravita- 
These expressions are similar to, though not tional constant, so that the theory becomes 
identical with, those in the K-model. The parameter free. 



fact that they are similar but not identical As was discussed in (Ghirardi et al. 



suggests that the involvement of gravity in 1990a) the QMUDL model has certain lim 



wave-vector reduction is strongly indicated, itations - it cannot deal with point particles 
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(for which case it leads to divergent densi- tion \ip) = a\a) + b\0) of two well-defined 

ties) and restricts itself to extended objects, stationay states |a) and \(3), having same en- 

The model parameters are such that it leads ergy E If gravitation is ignored, as is done in 

to an unacceptably high rate of energy in- standard quantum theory, the superposition 

crease during reduction. Furthermore, for = a\a) + b\/3) is also stationary, with the 

microscopic dynamics the reduction and lo- same energy E, 

calization process can lead to unacceptable jfj ®^ — E \ip) (214) 

processes such as excitation or dissociation ^ 

of nuclei. To avoid these problems. Ghirardi However, the inclusion of gravitation 

et al. proposed a CSL type modification of raiseg a new ques ti on: w h a t is the mean- 

QMUDL, and the introduction of a new uni- ing of the Schr5dinger time-evolution oper- 

versal length parameter. It is suggested that ator d/m7 There wiU be a nearly classical 

it does not seem possible to have a parame- spacetime associa ted with the state \a), and 

ter free theory for reduction, such as gravity a KiUing yector associated with it which rep . 

induced collapse. An alternate way out of resents tbe time displacement of stationarity. 

the difficulties of the otherwise very attrac- And there wiU be a different nearly c l assi . 

tive model of Diosi has been suggested by cal spacetime assoc i at ed with the state \/3), 

Penrose - we will recall this proposal next, and a dlfferent Kming vector associat ed with 

but find that here one is faced with a pos- it which represents the associa ted time dis- 

sibly new set of difficulties. Thus it would placement of sta tionarity. The two Killing 

seem that at present CSL might be the best vectorg can be identified with each other only 

model at hand, even though its fundamental if the two space -times can be identified with 

origin remains to be understood, and it yet each other point by point However; the prin . 

may have a strong connection with gravity dple of general covariance in general relativ . 

whose proper implementation remains to be ity forbids ^ gince the matter distribu . 

achieved. tions associated with the two states are dif- 
ferent, in the presence of a background grav- 

o Tn, j i t u itational field. On the other hand, unitary 

j. 1 he model oj Penrose ' J 

evolution in quantum theory requires and as- 



Penrose (Penrose, 1996, 1998) addressed sumes the existence of a Schrodinger oper- 
the question of the stationarity of a quantum ator which applies to the superposition in 
system which consists of a linear superposi- the same way that it appies to the individual 
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states, and its action on the superposition is 
the superposition of its action on individual 
states. There is thus a conflict between the 
demands of quantum theory and of general 
relativity. 

A tentative resolution is to make an ap- 
proximate point-wise identification between 
the two spacetimes, which in turn corre- 
sponds to a slight error in the identification 
of the Schrodinger operator for one spacetime 
with that for the other. This corresponds to 
a slight uncertainty in the energy of the su- 
perposition, for which it is possible to make 
an estimate in the case when the superposi- 
tion amplitudes are nearly equal in magni- 
tude. In the Newtonian approximation, this 
energy uncertainty Eq is of the order of the 
gravitational self-energy of the mass distribu- 
tion in the two superposed states. In accor- 
dance with the Heisenberg uncertainty prin- 
ciple, the superposition lifetime can be taken 
to be h/Eo, beyond which time the super- 
position will decay. In concept and in de- 
tail, this is quite like the damping time in 
Diosi's model. It is not clear here though, as 
to how the Born rule will be recovered dy- 
namically. 

Penrose notes the commonality with 
Diosi's ideas, the difficulties encountered by 
Diosi, and the resolution proposed by Ghi- 
rardi et al. by way of introducing a funda- 
mental length scale. Penrose observes that 



essentially the same difficulty would arise 
in his own approach too, because if one 
were dealing with point particles, the gravi- 
tational self-energy difference can become in- 
finitely high, implying instantaneous reduc- 
tion, which is clearly unreasonable. While 
Ghirardi et al. avoid this problem by intro- 
ducing a new length scale, Penrose proposes a 
different way out. The way out is in particu- 
lar based on noting that one needs to specify 
which states are the basic [stable] states, to 
which superpositions of basic states decay. 

It is proposed that the basic stationary 
states to which a general superposition will 
decay by state reduction are stationary so- 
lutions of the so-called Schrodinger-Newton 
equation (SN-equation). This equation is ac- 
tually a pair of coupled differential equations 
which are set up as follows, for a quantum 
mechanical particle of mass m moving in its 



own gravitational field (Diosi, 1984) 

dib h 2 

ih^- = V 2 ^ + 

ot 2m 

V 2 $ = AnGmWl 2 



(215) 



This system of equations has been ana- 



lyzed in (Bernstein et a/. 1998 Giulini and 



Grofiardt 2011 Harrison et al, 1998 , 2003 



Moroz and Tod 1999 Ruffini and Bonazzola 



1969). These equations are closely related 



to the Schrodinger-Poisson equations which 



have been studied for much longer (Lange 
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etal, 1995). surprisingly low value is at variance with 

At this stage an important difference with the much higher estimates coming from sim- 

the models of Karolyhazy and Diosi seems to pie analytical estimates [and also from the 

be that that unlike in the latter two mod- work of Karolyhazy and Diosi] and prompted 



els, where an intrinsic uncertainty in space- flGiulini and Gro6ardt[ |2011| ) to look at the 
time structure is assumed, here the impact problem closely. 

on the evolution of the quantum state is due , , . . , A ,. „ 

Various numerical studies, as well as 

to the particle's own gravitational field. Also, , , , n , 

heuristic estimates, show that the ground 

the system seems to be set up deterministi- , . r , , , 

state energy is of the order 

cally and the presence of a stochastic element i r 12 ™ 5 

E~-l±^- (216) 
is not evident, at least a priori. Thus one 8 h 2 

could ask as to the origin of the stochastic Th e width a of the mass distribution in the 
feature which actually drives the system to ground state is 

ofe2 

one of the stationary states, and the accom- aQ ~ (217) 

panying Born rule. Also, if the evolution is 

which we immediately notice coincides with 

deterministic and non-linear, there appears 

the phase coherence cell length in microscopic 

to be present the possibility of superluminal 

limit of the K-model. 

propagation. 

These issues apart, the SN system of B ^ introducing a length scale / the SN 
equations yields some very interesting re- ec i uation can be written in terms of a dimen- 
sults. Spherically symmetric stationary solu- slonless coupling constant 

Cm 3 / ( I \ f m \ 3 
tions have been found and their stability has j( = 2 = 2 f ) f | (218) 

, . , h2 \ L pJ \ m P J 



been investigated. A comprehensive recent 



analysis is given in (Giulini and Grofiardt 



2011). Their study was motivated in re- 



sponse to (Carlip 2008 Salzman and Car- 



One considers the time-dependent SN equa- 
tion for initial values given by a spherically 
symmetric gaussian wave-packet of width a 



Vi(r,f = 0)= (7ra 2 )- 3/4 exp(-^) (219) 



lip, 2006) - the SN equation induces a grav- ^ v ' ' v / V 2a 



itational suppression of expanding Gaussian There are thus two free parameters a and m, 

wave-packets, and it was suggested by Car- and one asks for the regions in this parame- 

lip and Salzman that the suppression [and ter space where significant inhibitions of the 

hence wave-vector reduction] becomes sig- usual free quantum dispersion occur. In an 



nificant already at m ~ 1600 amu. This important analysis (Giulini and GroBardt 
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2011 ) give four different analytical arguments induced collapse. Other considerations of 



to show that inhibition of the dispersion gravity induced collapse have been made in 



(Ellis et al, 


1984 


) and ( 


Percival 


1995) 



coupling constant K of (218) becomes of or- brief but elegant summary of gravity models 



der unity. This conclusion coincides with and trace dynamics is given by {Diosi, 2005) 



that of Karolyhazy and Diosi, and we believe In our view, gravity induced collapse is 

it leads us to an important inference: a promising physical mechanism for a physi- 

The models of Karolyhazy, Diosi and Pen- cal realization of spontaneous collapse. Fur- 
rose all agree that if the width of the quantum thermore, Trace Dynamics and its extension 
state associated with an object of mass m be- to space-time structure [treating space and 
comes greater than of the order h 2 /Gm 3 , the time as operators] provide a plausible mathe- 
quantum-to- classical transition sets in. matical avenue for rigorously developing the 

For the experimentally interesting a = stochastic theory of gravity induced sponta- 

0.5 fim this gives m to be about 10 9 amu. neous collapse. 

These results are further supported by nu- 
merical investigations of the SN equation in 

IV. EXPERIMENTAL TESTS OF THE 
THEORETICAL PREDICTIONS 



(Giulini and GroBardt, 2011). The authors 



A. Introduction 



also note that the coherence time, the time 
beyond which collapse takes place, can be 
brought down by reducing the grating period 

in a molecule interferometry experiment. For We have considered two classes of under- 

instance, for a mass of 10 11 amu and grating lying theories for dynamical collapse: Trace 

period of 0.5 /mi they report a coherence loss Dynamics, and gravity induced collapse. The 

time of 300 ms. phenomenology of Trace Dynamics manifests 

It is significant that while the Penrose ap- itself through models of spontaneous col- 

proach does not directly address the emer- lapse. If spontaneous collapse or gravity in- 

gence of the Born rule, it correctly predicts duced collapse is a possible explanation for 

the regime where the quantum to classical the measurement problem, then the experi- 

transition takes place, in agreement with the mental predictions of these models differ from 

other gravity models. those of standard quantum theory. Bounds 

This essentially completes our brief review can be set on the parameters of these models 

of the three well-known models of gravity- by requiring that their predictions should not 
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disagree with those observations which are 
well-explained by the standard quantum the- 
ory. On the other hand one can perform new 
experiments, such as diffraction experiments 
with large molecules, for which the predic- 
tions of these experiments differ appreciably 
from those of quantum theory. The results of 
such experiments could vindicate the modi- 
fied quantum dynamics [and specific values 
of the associated parameters] or rule it out. 
This section reviews the bounds on model 
parameters which come from known physical 
and astrophysical processes, and from diffrac- 
tion experiments that have been carried out 
in the laboratory or are planned for the near 
future. Experiments are discussed through 
Sections IV.B to IV. G. Bounds on sponta- 
neous collapse models will be discussed in 
Sec. IV. H and IV. I and those on gravity 
based models in Sec. IV. J. 

As we saw in Sec. II, a large vari- 
ety of collapse models has been proposed: 
QMUPL, GRW, CSL, dissipative and non- 
dissipative, white, colored, Markovian and 
non-Markovian. The overall task of con- 
straining these models is very extensive, 
given their large variety, and considering that 
a large variety of observations [laboratory 
and astrophysical] as well as table-top exper- 
iments, has to be considered. The subject is 
in a state of rapid flux, and still in a develop- 
mental stage. Here we try to do as complete 



a job as possible, relying on the analysis in a 
host of important papers that have appeared 



in the last few years (Adler, 2007 Adler and 



Bassi 




2009 


Adler and Ramazanoglu 


2007; 








Feldmann and Tumulka, 


2012 


Nimmrichter 



et al, 2011b Romero-Isart 2011). It should 



be noted though that the CSL model has 
received the maximum amount of attention, 
and we will focus mainly on CSL. 

The CSL model introduces two new pa- 
rameters, the rate constant A and the corre- 
lation length r c . If spontaneous collapse is a 
correct theory of nature, the values of these 
parameters must follow from some underly- 
ing fundamental principles and/or be deter- 
mined by experiments. As mentioned earlier, 
GRW chose A ~ 10~ 16 sec -1 and r c ~ 10~ 5 
cm, in order to be consistent with observa- 
tions, while Adler chose A ~ 10 _8±2 sec -1 
and rc — 10~ 5 cm. However, there is room 
for more general considerations and for es- 
tablishing the allowed part of the parameter 
space in the A — r c plane. In Sec. IV. G we 
will consider bounds coming from cosmology 
and in Sec. IV. H we will summarize other 
physical processes which constrain rc and A. 

In Sec. IV. C we will summarize exper- 
iments which directly test quantum super- 
position by interferometers. These experi- 
ments also test collapse models and gravity 
models. The rationale is that if we can ob- 
serve superposition at mesoscopic and pos- 
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sibly even macroscopic scale, the quantum 
dynamics does not need alteration. If in- 
stead we experimentally observe a quantum 
to classical transition such as the collapse of 
the wave function while convincingly reduc- 
ing all potential sources of noise, this would 
strongly hint that an alteration of the fun- 
damental equations of quantum mechanics is 
needed. The CSL model with Adler's value 
for A predicts a quantum to classical tran- 
sition at only two or three orders of mag- 
nitude away from present molecule matter- 
wave experiments. This and new proposals 
for optomechanics experiments with trapped 
bead particles [see Sec. IV.D and IV.E] bring 
experimental tests of CSL (amongst other 
proposals) within reach. The main focus of 
the following sections will be on those table- 
top matter-wave and optomechanics experi- 
ments. 

While we do not aim to give a complete 
review of all possible experimental tests of 
collapse models, we present the strongest cur- 
rent bounds in Table 1. 

B. Possible experimental tests of CSL 
basing on quantum superposition 

Matter-wave interference experiments 
such as molecule interference are approach- 
ing the mass limit for the quantum to 
classical transition in a 'bottom-up' fashion, 
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starting with particles, where quantum 
superposition is existing and pushing the 
limit upwards step by step. Nano- and 
micro-mechanical devices cooled to the 
quantum-mechanical ground state within 
optomechanics approach the problem from 
the top, starting at very massive objects 
namely mechanical cantilevers of hundreds 
of nanometer and even some micrometer 
size, which are sometimes even visible with 
the naked eye. The range in which both 
types of experiments will probably meet 
[to combine techniques, knowhow and ideas 
to overcome experimental hurdles and to 
switch-off known decoherence mechanisms] 
is from 10 nm to 100 nm in size [mass - 10 6 
- 10 9 amu]. The experimental aim is to show 
quantum superposition by negativity in the 
Wigner function of the motional states or 
by proving the wave-nature of such particles 
by single-particle interference. Interestingly, 
this size range is vital to test non-standard 
quantum theories such as CSL and gravity 
induced collapse. 

1. Collapse theory for diffraction experi- 
ments 

In order to understand how collapse mod- 
els differ from standard quantum mechanics, 
when applied to interferometric experiments, 
let us consider once again the QMUPL model 



of Sec. II. E, due to its simplicity. The multi- which we are primarily interested in. In this 



particle dynamics is given by Eqn. (36). By case, the damping behavior is less trivial than 

using Ito calculus, it is easy to show that the that of the QMUPL model. As we have seen 

master equation for the statistical operator in Sec. II. I, for small distances there is a 

p t = E[|^t)(^ t |] is: quadratic dependence of the decay function 

d i 1 on the superposition distance \x — y\, while 

-rpt = -j-[H,Pt\ - r 2^MQi, VluPtW- 

1 i=i for large distances such a dependence disap- 

Suppose - for simplicity sake - that all par- pears . T h e intermediate behavior is not easy 

tides are identical, and that during the time to unfold, but it can be conveniently modeled 

of measurement, the free evolution (given by by the following ansatz. Recall first that for 

H) can be neglected. Then, according to the a sing i e constituent the master equation 

above equation, the density matrix in the po- d i 

—pt{x,y) = --[H,pt{x,y)]-TasL(x,y)pt{x,y) 

sition representation evolves as follows: 

implies that the decay function is 

Pt (x,y)=p (x,y)e- XN ^ 2t / 2 . (220) 

r cSL (x) = A[l-e-^ 2 ] (221) 

This equation contains everything there is 

, i , , , , a . r n j i for one single constituent. See Fig. [2l for a 

to know about the effect of collapse models & & i_i 

. , r , • • , , 1 , r plot of r cSL vs x. Here we see how the two 

on mterterometric experiments, at least from ° b 

,i i c ■ t-,-£c i fundamental parameters of the CSL model 

the conceptual point of view. Different mod- 

i j-rr i • i-i i l • l j i. *i rri enter into play. For a many-particle system, 
els differ only m the technical details, lne r J J ^ J ' 

i ,. , ni ,i , ■ j one makes an ansatz and assumes that the 

above equation tells that, m order to mea- 

1, rr , ,. , above expression for the decay function holds, 

sure a collapse effect, corresponding to a sig- r J ' 

■ c , , f , c n ., • except that one has to multiply A by the 

mfacant damping factor, the following criteria ^ r j j 

,i .4.1, 4. u- appropriate numeric factorial, as described 

must be met: the system should be as big as ^ ^ ' 

•ii/i Ar x ., i 1, i i i . in Sec. ILL The numerical factor is n 2 N 

possible (large JM ); it should be created m a 

til „ ... /i i i \ v. ■ "u where n is the number of nuclei [called clus- 

large superposition (large \x — y\), which L 

■ .jr. i -ii terl within a volume of linear size rn, and 

is monitored after a time as large as possible 1 ' 

/i ,\ mr- • ,i iii i ii' j. £ is the number of clusters in the many-particle 
(large t). lhis is the goal that all interfere- J ^ 

metric experiments aim at reaching, in order s y s ^ em - 

to test the validity of the superposition prin- In a11 mterferometric experiments so far 
ciple, and thus also of collapse models. realized, the period of the grating is compa- 

rable to rc ~ 100 nm. Therefore, to extract 
We now come back to the CSL model, the significant order of magnitude, it is sufli- 
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Wave-packet size x (in meters) 

FIG. 2 This plot aims to compare experiment with model to map the significance of different 
experiments to test collapse models. The CSL decay function T is shown against spatial dimen- 
sion/separation in case of four different mass systems (represented here as four different y— axes), 
using A = A A dicr- The x-axis is not continuous, so as to permit combination of all experiments of 
very different mass and size scales in one plot. We aim for experimental parameters to fit to the 
part of the T curve where we observe a significant chance to test collapse models. Essentially we 
need a good mixture of mass, spatial separation and duration of superposition. For example, for 
molecule (10 4 amu, experiment done) and OTIMA metal cluster (10 6 amu, experiment proposed; 
bead experiments such as MERID are at the same range in the plot) interferometer which both 
give spatial separation on the same scale as x = rc, decay rates are 2.2 x 10~ 2 s , and 2.2 x 10 2 s _1 
respectively. In the case of the cantilever (experiment proposed) the experimental bound to V 
would be quite strong, but the spatial separation is small and quite far away from x = tq- Atom 
interferometers (experiment done) have an exceptionally large spatial separation, but the mass is 
small and therefore the bound to V is very weak. 



cient to work in the regime x 3> Tq- Taking have: 



into account Eqn. (80) (Adler, 2007) and 
that in the case of macromolecules N — 1 
(typical molecule size being about 1 nm), we 
have: 

r cSL ^ An 2 , (222) 
Since no interferometry-based experiments 
have so far detected any spontaneous collapse 
effect, this implies that the damping factor 
exp[— r cSL t] must be insignificant. We then 



A<l/n 2 t, (223) 
where n measures the number of nucleons in 
the system, and t the duration of the exper- 
iment. This is the type of bound that in- 
terferometric experiments place on the col- 
lapse rate A. The experiments do not pro- 
vide a bound on the second parameter of 
the CSL model rc, for the reasons explained 
above. More general situations could be con- 
sidered, but they have not been analyzed so 
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far. It will be desirable to carry out 
ful analysis of the allowed part of the A — tq 
plane, based on the data available from ex- 
periments, and to understand what role the 
grating period and size of the macromolecule 
will eventually play in bringing experiment 
and theory closer. 

The latest situation on the results from 
diffraction experiments has been discussed 



in (Feldmann and Tumulka, 2012; Nimm 



richter et ai, 2011b; Romero-Isart, 2011) 



The strongest current bound on A seems to 



be from the experiment of (Gerlich et al. 



2011) which sets A < 10 5 sec 1 for 



n 



7, 000. Adler estimates that an experiment 
with n = 500, 000 will confront the enhanced 
CSL value proposed by him, based on reduc- 



tion in latent image formation (Adler, 2007) 



Interestingly, in the same paper he also pro- 
poses to test whether 'latent image forma- 
tion' constitutes a measurement, by using a 
photographic emulsion 'which path' de- 
tector in one arm of a quantum interferome- 
ter. 

Spatial or centre of mass motion superpo- 
sition is needed to be demonstrated in exper- 
iments to test the quantum to classical tran- 
sition. As described previously atoms are so 
light that even the very large areas in to- 
day's atom interferometers do not increase 
the chance to test CSL (see Fig. 2). On the 
other hand the very massive cantilevers do 



not possess a large enough spatial separation 
(spatial size of superposition) to become good 
test embodiments for the quantum to classi- 
cal transition. It seems that the size range 
of particles of 10 nm to 100 nm, which corre- 
sponds to a mass range of 10 6 amu to 10 9 amu 
are ideal for such tests in matter-wave inter- 
ference experiments. 

In this section we focus on possible exper- 
imental tests of the CSL model, while similar 
or quite different experiments are possible to 
test different collapse models. 

C. Matter- wave interferometry: 
Molecule interferometry 

Experiments with matter waves exist 
since 1927 when Davissson and Germer 
diffracted a beam of electrons. It was the 
first proof of de Broglie's hypothesis on 
particle-wave duality. Since then matter- 



wave interferometry of electrons (Hassel 



bach, 2010), neutrons (Rauch and Werner 



2000), atoms (Cronin et ai, 2009) and 



molecules (Hornberger et ai, 2011) has a 



long and successful history to investigate fun- 
damental physics, and has been applied for 



metrology and sensing (Arndt et ai, 2011). 



Interestingly, a recent interpretation of atom 
interferometry experiments resulted in a de- 
bate on the possible detection of gravita- 
tional red shift by such tabletop experi- 



91 



ments flMiiller et oZ.l|2010b. 



We are here interested in centre-of-mass- 
motion interferometry (or de Broglie inter- 
ference) of very massive particles, as these 
experiments are promising to test modifica- 
tions of Schrodinger dynamics such as col- 
lapse models predicting a quantum to classi- 
cal transition at mesoscopic length and mass 
scales. The appearance of a single-particle in- 
terference pattern demonstrates wave-like be- 
haviour of the particles and can be seen as an 
indication for superposition. The full beauty 
of this particle position superposition can be 
seen from reconstruction of the Wigner func- 
tion of the motional quantum state by tomog- 



as those would stick to the surfaces. Typ- 
ically beam splitters for molecules are real- 
ized by gratings. Gratings are nanofabri- 
cated highly ordered periodic structures of 
freestanding nanowires made from metal or 
semiconductor materials or realized by stand- 
ing light fields using the Kapitza-Dirac ef- 



fect (Kapitza and Dirac 1933). Today the 



raphy (Kurtsiefer et al. 1997) 



tightest bound for the quantum to classical 
transition comes from molecule interferome- 
try. We shall give a brief history of molecule 
interferometry before we describe more de- 
tails of the work horse of molecule interfer- 
ometry - the Talbot-Lau interferometer. 

Beams of small molecules were first scat- 
tered at surfaces in the experiments by Es- 



Technically, to perform de Broglie inter- 
ference experiments, one has to overcome 
challenges of preparation of intense gas-phase 
beams, of preparation of spatial and tem- 
poral coherence of the matter wave, and of 
the efficient detection of the particles. Cen- 
tral to all experimental demonstrations of 
matter wave interference are optical elements 
which serve to coherently manipulate wave 
phases, and in particular to divide the wave 
fronts, thus creating different possible inter- 
ference paths. While bulk and surface crys- 
tals are well-adapted to diffract electrons and 
neutrons with de Broglie wavelengths in the 
range of 1 pm to 10 pm, it is impossible to use 
the same structures for atoms or molecules 



termann and Stern in 1930 (Estermann and 



Stern, 1930) followed by interferometry ex- 



periments with di-atomic molecules in the 
1990s. In 1999, matter-wave interferome- 
try with large neutral molecules was first 
demonstrated with the C 60 fullerene in Vi- 



enna (Arndt et al. 1999). Fraunhofer far- 



field interference was shown by using molec- 
ular diffraction at a single nano-fabricated sil- 
icon nitride grating with a grating constant 
of 100 nm. The beam was collimated by a se- 
ries of 5 fivcv slits to a beam divergence smaller 
than the expected beam diffraction angle of 
about 10 /irad. Only very few molecules orig- 
inally in the beam reached the diffraction 
grating and the detector and typical count 
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rates were of only very few molecules per 
second with a detection efficiency of around 
10 %. The resulting long integration time to 
resolve the interference pattern makes such 
experiments susceptible to noise. Prospects 
for large particle far-field interferometry and 
the related Poisson spot experiments can be 



tance between the gratings. 



found elsewhere (Juffmann et al, 2012b) as 



well as new developments of promising tech- 



niques for far-field experiments (Juffmann 



et al, 2012a). 



Talbot-Lau interferometer: Later 
molecule interferometry experiments were 
done with a so called Talbot-Lau interfer- 
ometer (TLI) to increase the beam intensity 
of the diffracted beam. A TLI is operat- 
ing in the near-field diffraction regime de- 
scribed by Fresnel integrals, where the spa- 
tial period of the diffraction grating and the 
interference pattern are on the same size 
scale. The scheme has been introduced by 
Clauser to cope with beams of low intensity 
and low collimation in interferometry exper- 



iments (Clauser and Reinsch, 1992). An ad- 




vantage of a TLI with respect to a Fraunhofer 
single-grating far-field interferometer is that 
the scaling of the distance between the grat- 
ings (Talbot length: Lt) is inversely propor- 
tional to the de Broglie wavelength XdB but 
quadratic with the grating period d, Lt = 
d 2 /\dB- This helps to compensate for small 
de Broglie wavelength by increasing the dis- 



FIG. 3 Different configurations of the Talbot- 
Lau interferometer are shown. (A) Three 
material grating as experimentally realized in 
(Brezger et al, 2002). (B) Kapitza-Dirac- 



Talbot-Lau interferometer realized in (Gerlich 



et al., 2007| ). (C) Optical Time-Domain Ionizing 
Matter Interferometer [OTIMA] as proposed in 



(Nimmrichter et al., 2011a) 



In more detail, the three grating TLI oper- 
ates with weakly collimated molecular beams 
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with divergence of about 1 mrad and accepts 
a large number of molecules in the initial 
beam contributing to the final interference 
pattern. The first grating prepares the beam 
coherence, while imprinting a spatial struc- 
ture on the molecular beam (see illustration 
of the Talbot-Lau interferometer in Fig [3]) 
acting as an absorptive mask. The second 
grating - the diffraction grating - is then si- 
multaneously illuminated by 10 4 individual 
coherent molecular beams. The second grat- 
ing generates a self-image at the Talbot dis- 
tance Lt- Therefore each of the 10 4 initial 
sources will be coherently mapped to the Tal- 
bot distance after the second grating. This 
Talbot effect results in a self-image of the sec- 
ond grating at half-integer multiples of Lt- 
The Lau effect makes an incoherent summa- 
tion of the individual coherent beam sources 
located at the first grating to an integrated 
signal. Basically, the number of molecules 
contributing to the final interference pattern 
is multiplied by the number of illuminated 
slits of the first grating. The third grating is 
then placed close to this Talbot position af- 
ter grating two and scans over the diffraction 
pattern perpendicular to the molecular beam 
to enable integrated signal detection. This 
is enabled as the period of the scan grating 
exactly matches the period of the Talbot self- 
images of the diffraction grating. Talbot, Lau 
and Talbot-Lau effects have been nicely illus- 



trated by recent optical experiments (Case 



et ai, 2009). The successful implementa- 



tion of different Talbot-Lau interferometers 
for molecules has been summarized in a re- 



cent review article (Hornberger et al, 2011), 



where more and detailed information about 
techniques and requirements can be found. 

A recent version of the TLI is the so- 
called Kapitza-Dirac- Talbot-Lau interferom- 
eter (KDTLI), which has been used to 
demonstrate interference of a 3 nm long di- 



azobenzene molecule (Gerlich et ai, 2007). 



Here the second grating was realized by an 
optical phase grating, where molecules are 
diffracted at periodic optical potentials due 



to the Kapitza-Dirac effect (Kapitza and 



Dirac 1933). The use of light gratings 



avoids the dispersive van der Waals (vdW) 
or Casimir-Polder (CP) attraction between 



molecules and gratings (Hornberger et al. 



2009), which is known to phase shift the in- 



terference pattern but also to reduce the vis- 



ibility (Hackermiiller et a/. 2003) due to dis- 



persive effects for molecular beams with finite 
velocity spread. The interaction effect scales 
with the particle velocity and particle polar- 
isability as well as the dielectric properties of 
the grating material. The details of the in- 
teraction potentials and related vdW-CP ef- 
fects and how those can be investigated by 
molecule interferometry experiments are still 



under intense investigation (Buhmann et al. 
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2012 Canaguier-Durand et al, 2012a[b ). Es- ing constant. However any particle traveling 



timates show that even with improved veloc- over that distance, even if it starts at zero 

ity selection schemes, where the width of the velocity, will be accelerated to higher speed 

selected velocity is below one percent of the than 1 m/s (namely 7 m/s over 2.5m) by 

mean velocity at full width half maximum, it Earth's gravity g acting over that distance 

is expected to disable interference with parti- on the particle. Slowing or compensation of 

cles with masses beyond 10 5 amu. Presently, acceleration by additional carrier fields or in 

the largest de Broglie interfered particle is an space experiments would be needed to over- 

about 7000 amu massive perffuoro-alkylated come this limitation while presenting a signif- 



6*60 molecule (Gerlich et al., 2011), which icant experimental challenge. Therefore, TLI 



gives currently the strongest bound on col- experiments (where the speed of the particle 

lapse models. For comparison all present ex- or equivalently it's wavelength has to have 

perimental bounds are listed in Table 1. a certain value between the gratings) with- 

The specifications of a TLI can be esti- out compensation of Earth's gravity are lim- 

mated easily. For instance the specifications ited to a P article mass of around 1Q7 amu - 

for interference of 10 6 amu massive particle That limit exists for a11 Possible orientations 

in a Talbot-Lau type interference with grat- of the interferometer to g. While that is true 

ings of period rf=100 nm: at a Talbot dis- the alternative single grating far-field inter- 

tance of L T =2.5 cm a particle velocity of ferometry is not limited by g. Diffraction of 

v=l m/s would be needed to be constant the matter-wave at the location of the grat- 

over this distance L T . [The size of the grat- in S' the separation of maxima and minima 

ing opening is limited to about 50nm by the of the interference pattern does not depend 

size of the nanoparticles which have to trans- on the s P eed of the Particles, but only on the 

mit as well as by available technologies for distance between the grating and the particle 



grating realization by light and from mate- 



detector. 



rial nanostructures] . For a higher mass the Technical challenges for mass scal- 

particle would need to be slower at the same ing: The quest is for new technologies which 

Talbot distance or alternatively the Talbot can efficiently control and manipulate the 

distance would need to be extended for the centre of mass motion of heavy particles. The 

same particle speed. Simple estimates show mature techniques of ion manipulation and of 

that for a particle of 10 8 amu we find al- optical tweezing are of particular interest to 

ready L T =2.5 m at the same speed and grat- scale the mass up to particles of 10 nm to 
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1 urn (mass of 10 6 amu to 10 10 amu) in diam- cooling technique has been realized for op- 
eter. All experiments have to be performed tically trapped beads of 1/im in diameter in 



under ultra-high vacuum conditions to avoid the field of optomechanics (Li et al, 2011). 



decoherence by collision. We will come back We will come back to this in Sec. IIV.DI 



to this in Sec. IV. F In particular the chal- Furthermore particles need to be structurally 



lenges are: stable to survive launch and detection pro- 
1) Generation of intense particle beams: cedures. This includes techniques to gener- 
Particles need to be slow if massive to keep ate gas-phase particles such as by thermal or 
the de Broglie wavelength within the range laser Voiced sublimation, laser desorption or 
for experimental possibilities (not much ablation, but also sprays of particles from so- 
smaller than pm). The ideal particle beam ^tions, etc. and the subsequent manipula- 
has a high phase space density, which means tion of such Particles to meet the coherence 
that many (ideally all) particles propagate at requirements of matter-wave experiments, 
the same speed. The beam needs to be highly 2) Beam path separation: This is the need 
collimated, which means that the transverse for coherent beam splitters and other matter- 
velocity needs to be as small as possible, wave optical elements. While different reali- 
ideally zero. [Very high beam collimation sations of beam splitters are known for cold 



^< 10/irad) would enable the conceptually atoms (Cronin et al., 2009), material and 



simpler far-field single grating interferome- optical gratings are only existing option for 

try.] In the wave picture this means that the large particles. The challenging part is the 

transverse or spatial coherence needs to be realization of gratings with a high enough 

high. All this could be achieved by cooling precision in periodicity. The demand on the 

techniques, that effect the centre of mass mo- periodicity is very high for the TLI scheme, 

tion of particles, which yet have to be devel- where the average grating pitch has to be ac- 

oped for complex particles. Very interesting curate within sub-nanometer scales between 

and promising approaches have been followed all gratings. This can only be realized so far 

in the last few years. This especially includes by sophisticated optical interference lithogra- 



the collisional buffer gas cooling (Maxwell phy techniques. For far- field gratings the de- 



et al. , 2005 ) as well as optical cooling tech- mand is lower and electron-beam lithography 



niques (Shuman et al., 2010). Both tech- with alignment pattern to avoid stitching er- 
niques have so far been demonstrated for di- rors is possible for fabrication. The ability to 
atomic molecules. Interestingly, a feedback form laser light gratings from retro-reflection 
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or other superposition of laser beams de- cles of increased size and mass. To resolve the 

pends on the intensity and frequency sta- interference pattern a spatial resolution on 

bility of the laser. The power of the laser the order of the grating period is needed for 

needs to be sufficient to form an optical po- near-field interferometry experiments, which 

tential strong enough to act as a phase grat- is elegantly realized in the case of the TLI 

ing. This is on the order of some 1W con- by the third grating. Also a high spatial 

tinuous power for fullerenes. The limited resolution of the detection is needed if the 

availability of stable and medium power UV particle beam is not velocity selected before 

and XUV (wavelength < 200nm) lasers lim- entering the interferometer gratings. This 

its the fabrication of grating periods by op- is important to select the temporal coher- 

tical lithography as well as the optical grat- ence which is given by the distribution of de 

ing periodicity to about lOOnm (opening of Broglie wavelengths of matter waves emitted 

about 50nm). About an order of magnitude by the source. Fluorescent molecules can be 

smaller grating periods can be possibly fabri- detected with single particle resolution and 



cated by electron-beam lithography or direct sufficient spatial resolution ( Juffmann et al. 



focus ion beam (including the novel He-ion 2012a). 



direct write) milling. Another limitation is In the following, we will discuss different 

that the grating area has to fit the size of the alternative approaches on how to possibly im- 

particle beam diameter which is on the order plement experiments to probe the quantum 

of 1mm. Not many fabrication techniques are superposition of particles. We will summa- 

capable of manufacturing precise gratings on rize different proposals for such experiments, 
that size scale. However, in combination with 
an efficient detector this dimension can be de- 

1. Neutral particles vs charged particles 

creased. 

3) Efficient detection of large particles: Quantum superposition experiments per 
Ideally, we want single particle detection res- se need to avoid any decoherence effect which 
olution. For example in the most recent is able to read out which-way information 
molecule interferometry experiments detec- and to localize the particle. A neutral par- 
tion is realized by ion counting after electron ticle is a natural choice for superposition ex- 
impact ionization, which is known to have a periments as the number of possible inter- 
very low ionisation efficiency (10~ 4 ). This has actions, which would enable a readout of 
to be at least kept at the same level for parti- which-way information, is reduced in com- 
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parison to charged particles. That is espe- 
cially true for superposition of slow particles. 
Therefore all interference experiments with 
molecules have been performed with neutrals. 
On the other hand the centre of mass mo- 
tion of charged particles can be manipulated 
and controlled to a higher degree by exter- 
nal electric and magnetic fields. This would 
be handy to prepare coherent particle beams. 
Here we will discuss the benefits and possi- 
bilities for charged particle interferometry as 
neutrals have been covered in the previous 
section. From the matter wave point of view 
we have to achieve the same parameter val- 
ues: de Broglie wavelength, periodicity of the 
diffractive element, etc. to observe an inter- 
ference pattern. This especially means that 
for a given (high) mass the speed needs to be 
rather low: for m=10 6 amu requires around 
v=l m/s. 

Trapping charged particles such as elec- 
trons or ions in Paul and Penning traps has 





2005 




Mikosch 



et aL, 2010; Willitsch et ai, 2008) such as 



with buffer gas cooled polyatomic ions in 



multi-pole traps (Gerlich 1995; Gerlich and 



a long and successful history (Paul, 1990). 
It has been used for studies of fundamen- 
tal physics such as the precise evaluation 



of physical constants (Brown and Gabrielse 



1986), for quantum information processing 



with one or many ions (Duan and Monroe 



2010 Leibfried et ai, 2003 Singer et ai 



2010), in chemical physics to investigate the 
kinetics and dynamics of chemical reactions 
on the few molecule level under controlled 



Borodi, 2009). The obvious benefit of us- 
ing charged particles for matter-wave experi- 
ments is the higher control over the motion of 
the particles. Guiding, trapping, and cooling 
is possible even for massive ions. For instance 
200 bio-molecules of 410 amu have been co- 
trapped with laser cooled atomic ions (Ba + ) 
and cooled to 150mK ( Ostendorf et al. , |2006 ) . 
This sympathetic cooling via Coulomb inter- 
action of laser-cooled atomic ions with molec- 
ular ions has been demonstrated to be effi- 
cient, however a difficulty which remains is 
to realize an ion trap which is stable for both 
species. The mass over charge ratio m/q 
must not be too different for both particles, 
which demands also a high control on the ion- 
isation technique for atom and molecule. 

Most of these techniques aim to spatially 
fix the ion in the trap to increase interaction 
times for spectroscopic and collision stud- 
ies or to cool the ions, while we are inter- 
ested in well controlled centre of mass mo- 
tion for interference. It might be difficult 
to achieve a coherent centre of mass mo- 
tion manipulation, but seems not impossi- 
ble also with respect to exciting new guid- 
ing techniques such as the recently demon- 



strated microwave-manipulation (Hoffrogge 
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et al., 2011), the manipulation of ions by effect (Batelaan, 2007). Electron interferon! 



light (Schneider et al. 2010), or multi-pole etry has been for instance used to investigate 



trap techniques (Gerlich, 1995). 



the Aharonov-Bohm effect (Tonomura and 



To scale up the mass of ions for exper- 
iments in order to test collapse models the 
very mature techniques of gas-phase cluster 
sources are available, such as sputter mag- 



Batelaan, 2009). The general understanding, 



netron sources (Haberland, 1994) or other no 



which is supported by experiments on for in- 
stance image charge decoherence effects, is 
that ions have to be very fast to not deco- 
here via one of the multiple interaction chan- 



ble gas aggregation sources with pick up for nels with the environment ( |Sonnentag and 



large molecules (IGoyal et al.\ 119921 iToenniesI |Hasselbach[ |2007[ ). The challenge will be to 



and Vilesov, 1998). Beams of such sources 



are intense since they are cooled by the su- 
personic expansion and the mass of a sin- 



gle cluster can be 10 amu and beyond (von 



Issendorff and Palmer, 1999). In combina- 



tion with quadrupole mass filters, which work 
very similar to ion Paul traps, metal clusters 
of very narrow mass distributions can be real- 



ized with m/Am=25 (Pratontep et al. , 2005). 
Additional techniques will need to be realized 
for deceleration of such big clusters, but as 
long as the particle is charged a high degree 
of control is guaranteed. 

While this is true massive ion interference 
has yet to be shown to work. A recent review 
article on ion interferometry lists that so far 
only electrons and the He + ion showed quan- 



avoid and shield all possible interactions of 
ions with surrounding matter and fields, such 
as for instance the coupling of the ion to its 
own image charge in a metal surface. 

As for the neutral particles in the case of 
a TLI the acceleration by Earth's gravita- 
tional field has to be compensated. Guiding 
potentials have to be extremely flat to not 
influence the superposition state, to not lo- 
calize the particle. External electro-magnetic 
fields have to be shielded by a Faraday cage of 
the right dimensions and materials, where re- 
cent technological progress has been made for 
the stabilization of magnetic fields in atom 



experiments (Gross et a/. 2010). Very sta- 



tum interference (Hasselbach, 2010). Typ 



ically electrons are diffracted at bi-prisms 
or solid surfaces as if applied for hologra- 



phy (Tonomura, 1987), but also light grat 



ings are possible utilizing the Kapitz-Dirac 



ble electrical power supplies will be needed 
for the cold Paul trap for ion beam gener- 
ation and an ion guide field. Electric stray 
fields from patch effects of adsorbed atoms 
and molecules at the shielding and elsewhere 
may be avoided as well as time varying elec- 
tronic inhomogeneities in the shielding ma- 
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terial. Edge fields of the guiding electrodes ter aggregation source as mentioned before, 
and other parts inside the shielding have to A further cooling/ deceleration device will re- 
be carefully considered. duce the velocity of the big clusters which is 

However, a simple estimate shows that for an existing technology for charged particles, 

instance all applied voltages would have to A chopper modulated particle beam can be 

be stabilized to the level of below 10 _ 10 V used for mass as well as velocity selection of 

for the time of interference which seems to the clusters in combination with a time of 

be impossible to be achieved at the moment, flight mass spectrometer (TOF-MS) detector. 

At present only the neutrals show success for The main inve ntion is a neutralisa- 

large particle centre of mass motion inter- tion/ionization scheme implemented as the 

ference. On the other hand interference at- interferometer. The neutralization of the 

tempts with larger particles suffer from the c i uste rs to enable a coherent propagation 

non-existence of guiding, slowing and cooling of the superpos ition state is planned to be 

techniques for neutrals. Therefore a clever ac hi eV ed by light-matter effects directly at 

solution for now is to try to take the best the light gra ti ngs . The scheme makes use 

of both the worlds: manipulation of charged of a sequen ce of three vacuum ultra-violet 

particles and interfering after neutralization, (VUV, A=157nm) ns long light pulses to re- 

which we describe in the following section. alize the interferometer gratings. The en- 
ergy of a single photon of about 8 eV is 

2. The compromise - combination of sufficient to ionize or neutralize metal clus- 



techniques for charged and neutral particles: ^ ers by photo detachment (Haberland, 1994). 



OTIMA These processes are also applicable to large 



bio- molecular complexes (Marksteiner et al. 



A novel three light grating Talbot-Lau 2009). The light intensity pattern realized 



scheme in the time domain aims towards the by three retro-reflected laser pulses hitting 
interference of particles of up to 10 9 amu as the propagating particles transversely at pre- 
proposed by the Vienna molecule interferom- cisely timed locations with respect to each 



etry group and described in (Nimmrichter other realize the TLI gratings with grating 



et al, 2011a). This interferometer is called period of d = A/2. The standing wave nor- 

optical time-domain matter-wave (OTIMA) mal mode pattern forms on one hand the 

interferometer. The charged particles will gratings but is also a spatial resolved ioni- 

be provided by a mass filtered metal clus- sation/neutralization device: the intensity in 

100 



the antinodes is sufficient to ionize or neu- 
tralize the particles while it is not in the 
nodes. Therefore clusters which pass through 
the antinode will get ionized while others in 
nodes will not. This in combination with 
electrodes to divide the beams of neutrals 
and ions is the realization of an absorptive 
grating. The first and third gratings need 
to be absorptive gratings, which means they 
need to spatially mask out parts of the cluster 
beam and are realized by intensity dependent 
ionisation of the molecules. The second grat- 
ing needs to be a phase grating and will be 
realized by the optical dipole force acting on 
the particle by making use of the Kapitza- 



D. Optomechanics: Cantilever 

Here we describe an experimental ap- 
proach which is alternative to matter wave 
interferometry. While the aim to understand 
the limitations of quantum mechanics is as 
old as quantum mechanics itself, the first pro- 
posals for a table-top experimental test by 
using the superposition or other non-classical 
states of massive mesoscopic or even macro- 
scopic mirrors have been published in the 



Bose et al. 


1997. 


1999 


Marshall 



Dirac effect (Batelaan ; 2007). 



Charged clusters will be neutralized by the 
first grating, diffracted at the second phase 
grating and ionized again for detection at 
the third grating. This is a promising at- 
tempt to realize matter wave experiments 
with very massive particles to test collapse 
models. More details about the OTIMA ap- 



proach to test CSL can be found in (Nimm 



richter et al. , 2011b). The spatial size of the 
superposition is estimated by the grating con- 
stant and is on the order of the CSL param- 
eter, namely rc=100 nm. Fig. 2 illustrates 
the bound on Ycsl while choosing Adlerss 
value for A and a cluster mass of 10 6 amu. 



et a/.[|2003[ ). 

The mechanical motion - the vibration - 
of the mirror, which was later realized by 
a nano-mechanical or micro-mechanical can- 
tilever, has to be prepared in the quantum 
mechanical ground state which is modeled 
by a simple harmonic oscillator (fc^T < hu>). 
The readout of the vibrating mirror is done 
by coupling it to a sensitive optical inter- 
ferometer to compare the light phase with 
a stabilized cavity. The conceptual idea is 
to first prepare the mechanical oscillator in 
the vibrational or phononic ground state |0) 
by cooling and then generate a coherent su- 
perposition of or with the first excited vibra- 
tional state |1) by single photon excitations. 

A high mechanical as well high opti- 
cal quality (Q) factor is needed to reach the 
regime of low dissipation to strongly couple 
optics to mechanics and to cool the device 
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ultimately to the ground state. While opti- 
cal control of cantilevers was under investiga- 
tion for quite some time it was only in 2006 
that two groups reported the successful opti- 



cal cooling of mechanical cantilevers (Gigan 



et al, 2006 Schliesser et al, 2006). Interest 



ingly, the cooling mechanism is very similar 
to the optical cooling of atoms: The optical 
resonance - in most cases an optical cavity 
resonance - is slightly detuned to the cooling 
sideband of the mechanical resonance, which 
can be in the range of 1MHz to 10GHz. This 
opens a cooling channel for the mechanics 
of the cantilever though the optical leakage 
of the cavity. Achieved temperatures corre- 
sponded still to a high phononic occupation 
- to many vibrational states occupied, but it 
boosted the rapid development of an exciting 
new field of research, namely opto- mechanics. 



This is summarized elsewhere (Aspelmeyer 



eTal] [20T0| |Kippenberg and Vahala] [2008 



Marquardt and Girvin, 2009). Also the first 



schemes on how to generate and probe the 
superposition state of a cantilever appeared 



(Kleckner et al, 2008). 



Only a few years later ground state 
cooling of micrometer sized structures has 



been achieved by Cleland's group (O'Connell 



et al. , 2010 ) and by the groups of Aspelmeyer 



and Painter (Chan et al., 2011). Advanced 



nanofabrication technology enabled realiza- 
tion of structures with both high mechanical 



a) 
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FIG. 4 Optomechanics. Upper Panel: Proto- 
type of optomechanically cooled cantilever as re- 



alized in ( Gigan et al. 2006 ) . Quantum op 



tical detection techniques enable the sensitive 
read out of vibrations as they couple to light 
fields. Lower Panel: Mechanical Resonator In- 
terference in a Double Slit (MERID) as proposed 
in ( Romero-Isart 2011). The centre of mass 



motion of a single optically trapped nanoparti- 
cle is first cooled and then superimposed by an 
optical double potential. The interference pat- 
tern evolves in free fall after switching off the 
trapping field. 



and high optical Q factors (10 5 ) in addition 
to clever optical or electronic readout tech- 
niques. This opens the door to many exciting 
quantum information processing and sensing 
experiments in the near future, but let us go 
back to our initial question if those structures 
can test our beloved collapse models. 

These structures are very massive, 10 6 
amu to 10 15 amu depending on their size, but 
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the vibration amplitudes when compared to 
the parameter tq of the CSL model are very 
small. This limits their ability to test col- 
lapse models and the parameter range to test 
CSL by such systems is indicated in Fig. [2} 
To investigate this a little further we esti- 
mate the spatial size Xq of this position su- 
perposition state by using the size of the zero 
point motion of a simple harmonic oscillator 
xq = \Jh/2mu, where m is the mass of the 
cantilever and u> its frequency in a harmonic 
potential. This spatial size of the ground 
state at 25 fiK is l-10~ 15 m for a typical micro- 
mechanical oscillator of a mass of 50 ng res- 



onating at about 1 MHz (Gigan et al, 2006). 

However, spatial superpositions to tests 
CSL have to be on the order of 10 nm or 
larger (rc = 100 nm), which is roughly 
seven orders of magnitude away from what 
micromechanical oscillators can achieve at 
the moment. Mass or frequency or combi- 
nations of both have to be improved by that 
amount, which is very difficult as for most 
materials mass and resonance frequency are 
coupled and depend on the spatial dimen- 
sions of the cantilever. To see vibrational 
state superposition larger than the quantum 
mechanical ground state the optomechani- 
cal device has to be driven in an extreme 
regime: eight orders of magnitude in mass 
or frequency at optical finesse of 10 6 . But 
there are very interesting systems providing a 



larger zero point motion such as carbon mate- 
rials with exceptional mechanical properties 



(Iijima, 1991 Novoselov et al, 2004): e.g. in 



dividual single-wall carbon nanotube oscilla- 
tors generate xq =1 pm, with m=8 x 10~ 18 g 
and at u=100 MHz at a ground state temper- 



ature of T=2.5 mK (Sazonova et al., 2004) 



Such systems have been used for mass sens- 



ing with hydrogen mass resolution (Chaste 



et al., 2012). One big challenge remains for 
such carbon materials, which is their very 
small absorption and reflection cross sections. 
This means it is not clear how to realize the 
needed high optical quality factor for opto- 
mechanics. But hope is to cool via other in- 
teraction channels possibly in the electronic 



Brown et al. 


2007 


Chen et al. 


2011 



Eichler et al., 2011). 



Another difficulty for the test of collapse 
models by cooling mechanical cantilevers to 
the ground state |0), is that the light field has 
to be switched on all the time. Otherwise the 
substrate to which the cantilever is coupled 
will rapidly heat back which is probably much 
faster than the collapse time. There is not 
much time for 'free propagation' of the super- 
position. New ideas on pulsed optomechanics 
may help to prepare and reconstruct quan- 



tum states of the mechanical motion (JVanner 



et al. , 2011 ) faster. So ideally we would prefer 
to use a massive harmonic oscillator which is 
realized without a link to any substrate. This 
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is what we discuss in the next section. 

Optomechanical superposition using sin- 
gle photon postselection and their detection 
with nested interferometers has been dis- 



cussed in (Pepper et al. 2011, 2012). 



Interestingly, mirror stabilization ideas are 
linked to the much larger interferometers for 



the detection of gravitational waves (Bragin- 



sky et al. , 2002), while in a different parame- 



ter range due to the much higher mass of the 
mirrors in use. 

E. Micro-spheres and nanoparticles in 
optical potentials 

Here we describe a new and very promis- 
ing route to test collapse models by gener- 
ating spatial superposition states of the cen- 
tre of mass motion of very massive nanopar- 
ticles and possibly even microspheres. This 
is a combination of optomechanics with cen- 
tre of mass motion superposition states as 
in matter-wave interferometry. Optically 
trapped particles represent an almost ideal 
realization of a harmonic oscillator as al- 



ready mentioned by Ashkin (Ashkin, 1970) 



and more recently re-discovered for cavity 



optomechanics (Chang et al., 2010; Romero 



Isart et al., 2010). In comparison to mechan- 



ical cantilevers as discussed in the previous 
section there is no mechanical link acting as 
a dissipation channel for the mechanical os- 



cillation in such systems if implemented in 
a vacuum chamber to avoid collisions with 
background gas particles. Therefore the me- 
chanical quality factor is very large. Such 
trapped particles can be seen as an optome- 
chanical system and techniques such as for 
cooling the oscillation - which is now the cen- 
tre of mass motion of the particle in the opti- 
cal trap - need to be implemented. Two cri- 
teria to test CSL and other collapse models, 
are fulfilled: a high mass of the particle and 
a large size of the superposition which can be 
comparable to rc- This makes such exper- 
iments a strong competitor to the OTIMA 
cluster interferometer, see Fig.2. 

In addition these systems enable free cen- 
tre of mass motion of the initially trapped 
particles after switching off the trapping field 
and after generation of the spatial super- 
position. This allows for combination with 
matter-wave interferometric techniques and 
schemes. Recently, some ideas have been put 
forth how to perform tests of quantum su- 
perposition with so-called beads (balls of di- 
ameter 10 nm to 10 fim made of glass or 
polystyrene) ( |Romero-Isart et~aT. , 2011b[c ). 
The basic sequence for such experiments is 
to first optically dipole trap a single particle, 
use optical techniques to cool the centre of 
mass motion of the bead in the optical trap- 
ping potential. The next step is to generate a 
superposition state of the particle position by 
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a double-well optical potential by single pho- 
ton addressing of the first excited vibrational 
state, as theoretically described within cav- 
ity quantum electrodynamics (QED). Cool- 
ing must be sufficient to increase the size of 
the particle wave-packet to overlap with both 
wells, so that there is an equal probability to 
find the particle left or right - the coherent 
superposition by a measurement of a squared 
position observable. After switching off the 
trapping potential in free fall the spatial den- 
sity distribution of the particle in multiple 
subsequent experiments can be mapped and 
evaluated for a quantum signature by for 
instance state tomography through Wigner 



viruses by optical fields has been pioneered 
by Ashkin and others since the 1970s and is 
now very broadly applied in many fields of 



science (Ashkin, 1970 Ashkin and Dziedzic 



1987). Techniques which are typically sum- 



marized by the term optical tweezing include 
the broad fields investigating optical angu- 



lar momentum (Allen et al., 2003), optimiz 



ing the trapping, levitating and guiding of 
single dielectric particles by optical gradi- 
ent and scattering forces in various geome- 



tries (Ashkin, 2006 Chu, 1998) including the 



guiding through hollow core photonic crys- 



tal fibres (Benabid et al, 2002) and opti 



cal binding (Dholakia and Zemanek, 2010). 



function reconstruction ( Romero-Isart et al. 



2011c) or much simpler by interference pat- 



tern detection at a fixed detector position. 
This is basically a double-slit experiment ap- 
plied to very massive objects (polystyrene 
bead of about 30 nm diameter has a mass 
of 10 6 amu). To get a significant detection 
statistics the single particle experiment has 
to be repeated many times. The connection 
to test collapse models is worked out in detail 



Ashkin and co-workers demonstrated already 
the trapping of polystyrene and glass micro- 
spheres, of viruses and bacteria and even of 
complete cells in solutions and high- vacuum. 
They developed a vacuum loading system 
and they demonstrated the stable levitation 
of particles at a vacuum of 10~ 6 mbar for half 
an hour by a feedback stabilization technique. 
A very detailed summary of this field can 



be found in Ashkin's book (Ashkin, 2006). 



in (Romero-Isart, 2011) and the experiment 



is called MERID which is the shortcut for 
mechanical resonator interference in a dou- 
ble slit. 

The manipulation of microscopic particles 
as silica and polystyrene sphere but also bi- 
ological cells and even living organisms as 



The particle size is typically limited to be 
not smaller than 1 /mi to form a stable trap, 
while optical near-field techniques have been 
very recently used to trap single nanoparti- 



cles with the help of plasmonic (Juan et al. 



2009) or photonic crystal structures (Rah 



mani and Chaumet, 2006) in solution. An 
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application of such advanced trapping tech- 
niques in vacuum has to engineer the chal- 
lenge of particle-surface van der Waals (vdW) 
and Casimir-Polder (CP) interactions or in 
turn could be used to investigate those inter- 
actions. However it has been demonstrated 
very recently that even 30 nm particles can 
be optically trapped in tightly focused free 
beams under vacuum conditions when gra- 
dient forces dominate scattering forces and 



vacuum conditions to dissipate the kinetic en- 
ergy of the trapped particle. All experiments 
- such as the competing cluster and molecule 
interferometry experiments - have to be per- 
formed at ultra-high vacuum (UHV) condi- 
tions (p < 10 _10 mbar) to avoid collisional de- 



coherence of the superposition state (Horn 



berger et al, 2003). A further challenge 



with parametric stabilization (Gieseler et al. 



is that the interferometer has to be stable 
over the duration of many single particle ex- 
periments. One idea is for experiments in 



2012). 



space (Kaltenbaek et al, 2012). Centre of 



Experimental challenges: Cooling is again 
the key for this experiment. Here one needs 
to cool the centre of mass motion of a bead 
in an optical field ideally to the ground state: 
Ashkin pioneered the feedback stabilization 



mass motion trapping would also be possible 



with ions (Leibfried et al., 2003). The elec 



(Ashkin, 2006). The Doppler cooling us- 



ing whispering gallery modes of the parti- 



cle has been proposed (Barker, 2010). Re- 



cently, the cooling of the centre of mass mo- 
tion of a single 1 fim glass bead to 1 mK 
has been achieved by a fast feedback stabi- 



lization technique (Li et al, 2011) as well as 



the optical parametric stabilization of a sin- 
gle silica nanoparticle (30nm) under vacuum 



trie trapping fields would replace the optical 
trap, but optical fields would still be needed 
for the cooling. While this is true free prop- 
agation of the charged particles, such as af- 
ter switching off the trap in the protocols for 
bead superposition experiments as explained 
above - would not be possible. The Coulomb 
interaction will certainly dominate the mo- 
tion of the particle - it will not be a free 
motion. On the other hand a recent pro- 
posal with magnetic levitated superconduct- 
ing particles claims feasibility for large super- 



conditions (10 mbar) at 400mK (Gieseler 



positions (Romero-Isart et al., 2011a). 



et al, 2012). These are the first promising 



steps to realize the proposed experiments to 
test superposition of such large and heavy 
particles. Importantly, feedback stabilization 
techniques will enable to trap beads under 



To avoid the difficulty of ground state 
cooling, one possibility is to use the Talbot- 
Lau interferometer scheme. Here, as we know 
from molecule interferometry, the require- 
ments on cooling are lowered as a quantum 
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interference effect can be observed at low spa- 
tial coherence of the matter wave. Centre 
of mass motion temperatures of lmK (for a 
given TLI geometry) would be sufficient to 
observe interference. This will work with a 
single particle source, but also many parti- 
cles in parallel traps would be possible which 
could significantly reduce the operation time 
of the interferometer and therefore lower the 
stabilization criteria on the interferometer. 
We are looking forward to see more exciting 
developments in this rapidly progressing field 
of research. 

One significant advantage of cavity op- 
tomechanics with trapped particles is the in 
principle very large separation of left and 
right for the superposition state which can 
be tuned by the optical field. Furthermore 
the optical field can be switched off and the 
particle can propagate in free space - showing 
the signature of superposition: an interfer- 
ence pattern in the spatial distribution. The 
size of the beads in the proposed experiments 
is on the order of lOnm to lOOnm, exactly 
the same size and mass range where cluster 
matter-wave experiments such as OTIMA are 
heading towards. 

F. Environmental Decoherence 

While the aforementioned and discussed 
collapse models can be seen as an exotic de- 



coherence mechanism we here discuss deco- 
herence effects of the environment interact- 
ing with the particle in superposition. Colli- 
sions with background particles and thermal 
radiation of the superimposed particle itself 
are counted as the major processes to local- 
ize the superimposed particle. Both decoher- 
ing effects affect all the different experimen- 
tal schemes to perform mesoscopic quantum 
superposition experiments and set limits on 
particle (and experimental setup) tempera- 
ture as well as background pressure inside 
the vacuum chamber depending on the size of 
the particle. According to decoherence the- 
ory the superposition state is destroyed and 
the particle is entangled with the environ- 
ment whenever any interaction of the super- 
imposed particle with the environment has 
the sufficient resolution to localize - to mea- 
sure the position of - the particle, which-way 
information is read out. 

We note, that an intrinsic difficulty with 
the test of collapse models is that it is clear 
how to falsify a proposed model with respect 
to predicted parameters. If on the other hand 
no interference pattern is shown by the exper- 
iments all systematic effects related to envi- 
ronmental decoherence have to be excluded 
as reason for the quantum to classical transi- 
tion. Here tuning of one of the test parame- 
ters such as mass of the particle or size of the 
spatial superposition will help to study the 



environmental decoherence effects. 

Mathematically, decoherence is described 
(as for the case of collapse models) by the ef- 
fect on the off diagonal elements of density 
matrix of the system including the particle 
and the environment which are reduced by 
the decoherence effect as given by the mas- 
ter equation very similar to Eqn. |221 The 
effect is evaluated by the decoherence rate 



tion. The emission of thermal photons is seen 
as the most important effect as the internal 
temperature of the particle is typically higher 
than the temperature of the environment. As 
an example for C70 fullerenes there is still 
full quantum contrast for emission of thermal 
photons by the fullerene at about 1500K, as 



experimentally observed (Hackermuller et al. 



2004 Hornberger et al, 2005). The interfer- 



function T as given in Eqn. |221| More de- 
tails on the concept and formalism of stan- 
dard decoherence theory can be found in the 
references given in Section I. More details on 
estimations of decoherence effects and asso- 
ciated decay rates for superposition experi- 



ence visibility is rapidly reduced for temper- 
atures at above 2000K where the wavelength 
of the emitted photons is comparable to the 
size of the superposition which was about 
lfim in this experiment. For a more detailed 
discussion of this long and short wavelength 



ments can be found in ( Romero- Isart, 2011) 





2010 




Romero- 



and (Nimmrichter et al, 2011b). Both pro- Isart 2011) 



cesses, collision and black body photon de- 
coherence, have been experimentally investi- 
gated and compared to theory with fullerene 



Romero-Isart estimates in (Romero- Isart 



interference (Hackermuller et al, 2004; Horn 



berger et al, 2003). We will here summarize 
the most recent estimate from this literature 
to give boundaries to the experiments. 

1. Thermal decoherence 

The emission, absorption, scattering of 
thermal - black body - radiation by the parti- 
cle in superposition can localize the particle 
if the wavelength of that light is compara- 
ble or smaller than the size of the superposi- 



2011) an emission localization time, which is 
inverse to the superposition decay rate, of 
100ms at a temperature of 100K for a 50nm 
particle, but claims that this time is indepen- 
dent of the particle size. If this claim is cor- 
rect it contradicts the observation of fullerene 
interferometry at 1500K. An extrapolation of 
this relation to mesoscopic particles (10 6 amu 
- 10 8 amu) gives temperatures between 800K 



and 200K (Nimmrichter et al, 2011b). In 



any case the predicted temperatures will have 
to be reached for the particle and the envi- 
ronment. This may require the cooling of 
the internal degrees of freedom of the par- 
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tide which is an experimental challenge, but 
buffer gas techniques are in principle appli- 
cable to any particle and cool all degrees of 



freedom (Gerlich and Borodi 2009 Maxwell 



et al, 2005). 



2. Collision decoherence 

Here the collision of the superimposed par- 
ticle with any other particle present will read 
out which-way information. Such collision 
decoherence processes have been studied in 
depth for fullerene experiments and an elabo- 



rate theory has been developed (Hornberger 



2006 Hornberger et al., 2004 2003). Ap- 



plied to the mesoscopic range of 10 6 amu to 
10 8 amu particles in OTIMA this gives min- 
imum required pressures between 10 _8 mbar 



and 10 mbar (Nimmrichter et al., 2011b) 



Romero-Isart estimates for a lOOnm sized 
particle and collisions with (N 2 )-molecules 
at 10~ n mbar, a decoherence time of about 
100ms in MERID. A more detailed parameter 



pends on the mean free path of the particle 
under the given vacuum conditions and there- 
fore the size of the particle. 

Certainly, more work on the theoretical 
side is needed to investigate those decoher- 
ence effects further. This will be an impor- 
tant guidance for experiments. For now it 
seems that collision decoherence can be con- 
trolled for mesoscopic particles while main- 
taining extreme UHV conditions in the ex- 
periments. On the other hand thermal ra- 
diation can become a more serious issue for 
larger particles of 10 8 amu. 

G. Concluding Remarks on Laboratory 
Experiments 

Interference of beads levitated in optical 
fields and interference of large metal clus- 
ters are both promising experimental routes 
to test collapse models. Clearly there is a 
certain possibility that other experimental 
routes or variations and combinations of the 



set is given in (Romero-Isart, 2011) and is in two main proposals: OTIMA (Nimmrichter 



agreement with the above values for OTIMA. 
This pressure is possible to achieve in ultra- 
high vacuum experiments. The parameter set 
also means that single experimental sequence 
from preparation of the coherence through 
superposition and detection has to be done 
in 100ms, which seems feasible in OTIMA as 
well as MERID. This estimate strongly de- 



etal, 2011a) and MERID (Romero-Isart 



2011) are successful in observing a meso- 



scopic single particle superposition state. 
OTIMA and MERID are the most advanced 
experimental attempts reported in the liter- 
ature at this time. 

In our opinion an experimental test of col- 
lapse models such as CSL with the Adler 
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value for A and a mass bound of 10 6 amu 
is within reach in the next 5 to 10 years. 
This will only be possible with intense re- 
search and development of new technologies 
for the handling of mesoscopic - lOnm to 
lOOnm sized - neutral and charged particles. 
Conditions to control environmental decoher- 
ence seem feasible to be reached in the ex- 
periments. We hope to see a scientific com- 
petition to probe this quantum to classical 
transition in the coming years. It will be in- 
teresting to see if quantum mechanics again 
survives. 

H. Cosmological Bounds 

As we saw in Sec. II. F, stochastic collapse 
leads to a secular increase in the energy of a 
system. For a group of particles of mass M 



thermal equilibrium because the cooling due 
to the adiabatic expansion of the Universe 
and the recombination of the plasma is bal- 
anced by the energy input into the IGM that 
comes from astrophysical processes such as 
supernova explosions and quasars. An upper 
bound on the stochastic parameter A can be 
obtained by assuming that all the heating of 
the IGM is from the stochastic heating of pro- 
tons and this gives that A should be smaller 
than about 10~ 8 . More detailed discussions 
of cosmological and astrophysical bounds can 



be found in (Adler, 2007 Feldmann and Tu 



mulka, 2012) 



A subject that is recently beginning to 



draw attention (De Unanue and Sudarsky 



2008 Landau et al. 


. 2011 




Leon and Su- 


darskyl 2011 


Perez et al., |2006; Sudarsky, 



the rate of energy increase is given by ( Adler 
dE 3 h 2 M 



2007) 



dt 



4 r% m% 



(224) 



See also (Bassi and Ghirardi 2003 Pearle 



and Squires, 1994). If there is no dissipa- 
tion in the stochastic collapse model, such an 
energy deposit will heat the system and the 
absence of the observed heating can be used 
to put upper bounds on A. 

An important case is the ionized inter- 
galactic medium (IGM), which has a temper- 
ature of about 2 x 10 4 between the redshifts 
of z = 2 and z = 4. The IGM is kept in 



2007, 2011) is the possible role of wave- 
function collapse in the very early Universe. 
A possible mechanism for the generation of 
primordial density fluctuations which even- 
tually grow to form large scale structures 
is provided by the hypothesized inflationary 
epoch in the very early history of the Uni- 
verse, just after the Big Bang. Inflation may 
have been driven by a scalar field and the 
zero point fluctuations of the quantized scalar 
field serve as a possible source for generating 



the requisite density inhomeogeneities (Lyth 



and Liddle, 2009). But how do these quan- 



tum fluctuations become classical, as the Uni- 
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verse evolves? Decoherence accompanied by 
the many worlds interpretation has been pro- 



posed as one possible solution (Kiefer and 



Polarski, 2009). Another possibility is that 



classicality is introduced by the models of 
stochastic collapse reviewed here, and it will 
be important and interesting to understand 
what sort of bounds are placed on the CSL 
parameters by the quantum-to-classical tran- 
sition of density fluctuations in the very early 
Universe. 

I. Bounds from other physical pro- 
cesses 

The standard GRW and CSL values for 
the model parameters were reviewed in Sec. 
II. I [including the enhanced value for A pro- 



stochastic collapse, giving A < 10 -3 sec -1 , 
(ii) excitation of bound atomic and nuclear 
systems [cosmic hydrogen should not decay 
during the life-time of the Universe] : A < 1 
sec -1 ; (hi) proton does not decay : A < 10; 
(iv) rate of spontaneous 11 keV photon emis- 
sion from Germanium : A < 10 -11 sec -1 , (v) 
effect on the rate of radiation from free elec- 
trons : A < 10 -5 sec -1 . 

Another interesting result is that of 



(Jones et al. 2004), which uses the Sudbury 



posed by (Adler, 2007) based on latent im 



Neutrino Observatory data to place a limits 
on the ratio of collapse rates of neutron and 
proton. The result this analysis is that the 
ratio of neutron to proton collapse rates is 
equal to (neutron mass/proton mass) ±.008. 
That is, mass proportionality to 1% accuracy. 
We also mention a proposal of an experiment 
to test the anomalous random walk due to 



age formation in a photograph, and flBassil collapse (Collett and Pearle, 2003), which 



et al, 2010a) based on image formation in 



the eye]. Earlier in this section we discussed 
bounds coming from diffraction experiments 
and from cosmology. A few other upper 
bounds have been placed too, taking into ac- 
count how some other processes would be 



affected (Adler, 2007). In so far as r c is 



concerned, tentative but plausible arguments 
have been given that it should be in the range 
10 -5 to 10 -4 cm. 

Amongst the processes studied thus far 
are : (i) decay of supercurrents induced by 



however has not been performed so far. 

Spontaneous photon emission: Accord- 
ing to standard quantum mechanics, a free 
charged particle travels along a straight line 
and does not emit radiation. According to 
collapse models, the same particle — though 
being "free" — always interacts with the noise 
field. It undergoes a random motion and, be- 
ing charged, it emits radiation. In a similar 
way, also a stable atom emits radiation. Not 
only because it undergoes a Brownian motion 
in space, but also because its electrons have 
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a non-negligible probability of being excited 
and subsequently de-excited with the emis- 
sion of photons. Therefore, collapse models 
predict the spontaneous emission of radiation 
from matter. 

The emission rate has been computed 
to first order perturbation theory using the 
mass-proportional CSL model, both for a free 



particle (Fu, 1997) (see (Collett et al, 1995 



Pearle et al, 1999) for a previous analysis) 



and for an hydrogen atom (Adler and Ra- 



mazanoglu 2007). In the first case, the pho- 



ton emission rate per unit photon's momen- 
tum is: 



dk 



e 2 \h 



dk 



1-1 + 



(225) 

where e is the electric charge, eo the vacuum 
permittivity, mo the nucleon's mass and k the 
emitted photon's momentum. In the second 
case, the formula changes as follows: 

dk 
(226) 

where ao is Bohr's radius. For small k this ex- 
pression is suppressed with respect to the rate 
of a free particle (the electron and the proton 
radiation rates add incoherently), while for 
large k is approaches twice the free-particle's 
rate. 




Comparison with experimental data (Fu 



the collapse parameter A of the CSL model, 
only 6 orders of magnitude away from the 
GRW value. Therefore, it excludes an en- 
hancement of this value of 8 orders of magni- 
tude proposed by Adler. However, as proven 



in (Adler and Ramazanoglu 2007) the emis- 



sion rate strongly depends on the type of 
noise. In particular, for a colored noise the 
emission rate is equal to that of the white 
noise, times the Fourier transform 7(0;*.) of 
the correlation function of the noise: 



dYi 



dk 



dT t 



Colored Noise 



dk 



White Noise 



(227) 



where u k is the frequency of the emitted pho- 
ton. For example, a cut off in the frequency 
spectrum of the noise field at 10 18 Hz would 
highly suppress the emission rate and restore 
compatibility between Adler's value. This is 
a rather high cut off, much higher than typi- 
cal cosmological ones (~ 10 n Hz). Therefore, 
it is reasonable to assume that the 'physical' 
emission rate — assuming that collapse mod- 
els provide a correct description of physical 
phenomena — is lower than predicted by the 
standard mass proportional CSL model. 

Table I summarizes the bounds on the 
CSL parameter A coming from various lab- 
oratory experiments and cosmological data. 



1997) (see (Collett et al, 1995; Pearle et al, As a final note, although spontaneous pho- 



1999) places a very strong upper bound on ton emission currently provides the strongest 
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upper bound on the collapse parameter A, 
macromolecule diffraction experiments seem 
to represent the most significant type of tests 
of collapse models, not only because they 
directly test the superposition principle of 
quantum mechanics, but also because they 
are less sensitive to the type of collapse model 
(dissipative or non-dissipative, with a white 
or colored noise field). 

J. Tests of gravity induced collapse 

Experiments on molecule interferometry 
and optomechanics are per se also test of 
gravity based collapse models: if a viola- 
tion of quantum superposition were to be ob- 
served, the next task would of course be to 
analyze which of the collapse models is indi- 
cated - CSL, gravity, or perhaps something 
entirely different. Another test, which has 
received considerable attention in the litera- 
ture on K-model cited above, is to look for the 
anomalous Brownian motion induced by the 
stochastic reductions. Such motion, which of 
course could also be induced by spontaneous 
collapse, seems too tiny to be detectable by 
present technology, but further careful inves- 
tigation into current technological limitations 
is perhaps called for (?). 

The optomechanical cantilever experiment 



ticular attention with regard to gravity in- 
duced collapse. Related discussions on this 



experiment can be found in (Adler et al. 



2005 Bassi et al., 2005a Bernad et al, 2006). 



An experiment to establish violation of 



Bell inequalities has been carried out (Salart 



et al., 2008), assuming that the time of col- 



lapse is as determined by gravity induced col- 



lapse (Diosi 1987) 



proposed in (Marshall et al, 2003) and dis- 



cussed above in Sec. IV. D has received par- 



It is not clear at this stage that there is 
a unique experimental signature of gravity 
models which will distinguish it from gravity- 
independent models of spontaneous collapse. 

V. SUMMARY AND OUTLOOK 

In the early years following the devel- 
opment of quantum theory in the 1920s 
the Copenhagen interpretation took shape. 
Dynamics is described by deterministic 
Schrodinger evolution, followed by a proba- 
bilistic evolution when the quantum system 
interacts with a classical measuring appa- 
ratus and quantum superposition is broken. 
An artificial divide was introduced between 
a quantum system and a classical measur- 
ing apparatus, in order for one to be able 
to interpret results of experiments on atomic 
systems. While widely accepted, even in its 
early years the Copenhagen interpretation 
had worthy detractors including Einstein and 
Schrodinger, to whom it was immediately ap- 
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parent that quantum theory by itself never 
says that it does not apply to large macro- 
scopic objects, and a direct consequence is 
paradoxes such as Schrodinger's cat. 

Two broad classes of attitudes developed 
towards the theory. One, given the extraor- 
dinary success of the theory, was to not ques- 
tion it at all: since no experiment to date 
contradicts the theory, one should accept the 
Copenhagen interpretation and the associ- 
ated probability interpretation as a recipe for 
making predictions from theory and compar- 
ing them with experiment. 

The other was to take serious note of the 
following difficulties : (i) classical macro- 
scopic systems are also quantum systems, 
and the quantum-classical divide introduced 
by the Copenhagen interpretation is vague; 
(i) the observed absence of macroscopic po- 
sition superpositions is in conflict with a 
straightforward interpretation of the quan- 
tum superposition principle; (iii) Schrodinger 
evolution being deterministic, it is 'paradox- 
ical' that probabilities should show up when 
one tries to describe the outcome of a mea- 
surement. 

As appreciation of these difficulties grew, 
the Copenhagen interpretation took a back- 
seat, and today it is perhaps fair to say that 
the interpretation is no longer considered vi- 
able, and should be permanently put to rest, 
having well served its purpose in the early 



phase of quantum theory. 

What has emerged on the scene instead, 
is three classes of explanations which ad- 
dress the difficulties mentioned in the previ- 
ous paragraph: 

[i] Do not modify quantum theory, 
but change its interpretation: This is the 
many-worlds interpretation. Quantum lin- 
ear superposition is never broken, despite ap- 
pearances. The different outcomes of a mea- 
surement are realized in 'different' universes, 
which do not interfere with each other be- 
cause of decoherence. It seems to us that 
in this interpretation it is not easy to un- 
derstand the origin of probabilities and the 
Born probability rule. Moreover, it is not 
clear when the multifurcation occurs. 

[ii] Do not modify quantum theory, 
but change its mathematical formula- 
tion: This is Bohmian mechanics. There are 
additional degrees of freedom — particles' po- 
sitions in space — whose introduction implies 
that outcomes of measurements can in prin- 
ciple be predicted beforehand, and probabil- 
ities can be avoided. 

[iii] Modify quantum theory: Replace 
quantum theory by a different theory, which 
agrees with quantum theory in the micro- 
scopic limit, agrees with classical mechanics 
in the macroscopic limit, quantitatively and 
dynamically explains the absence of macro- 



scopic superpositions and the emergence of 
probabilities, and whose experimental predic- 
tions differ from those of quantum theory as 
one approaches the mesoscopic and macro- 
scopic regime. 

In so far as the empirical situation is con- 
cerned, all three explanations are accept- 
able today. The many-worlds interpretation 
is in fact perhaps the favoured establish- 
ment viewpoint, because it involves minimal 
change in standard quantum theory: every- 
thing can continue to be as such, and that 
which is not observed is attributed to paral- 
lel branches of the Universe which cannot be 
observed. 

We have here proposed that the third 
avenue mentioned above be pursued: mod- 
ify quantum theory. What happens dur- 
ing a quantum measurement is a stochas- 
tic process. Even though the initial condi- 
tions and evolution for a microscopic system, 
successfully described by Schrodinger evolu- 
tion, are completely deterministic, the out- 
come of a measurement is completely ran- 
dom! A straightforward resolution would 
be to face the evidence head-on and de- 
clare that in the dynamics, determinis- 
tic Schrodinger evolution competes with 
stochastic evolution/reduction. For micro- 
systems Schrodinger evolution completely 
dominates over stochastic reduction. For 
macro-systems stochastic reduction domi- 



nates Schrodinger evolution, giving evolu- 
tion the effective appearance of Newtonian 
mechanics. Somewhere between the micro- 
and the macro- the Schrodinger evolution 
becomes comparable in strength to stochas- 
tic reduction. In this regime, which experi- 
ments are now beginning to probe, new phys- 
ical phenomena are predicted, which can be 
explained neither by quantum theory, nor 
by classical mechanics. These predictions, 
which are vulnerable to falsification, are also 
the strengths of a modified quantum theory. 
They are benchmarks against which the do- 
main of validity and accuracy of the standard 
theory can be verified in the laboratory. 

To this effect, the quantitative phe- 
nomenological models of Spontaneous Col- 
lapse, such as GRW, CSL, QMUPL and 
others, have been rigorously defined within 
the well-defined mathematical framework of 
stochastic dynamics. The models success- 
fully incorporate a Schrodinger type evolu- 
tion, and a stochastic evolution - the demand- 
ing requirements of non-linearity, causality, 
non-unit arity and norm preservation are suc- 
cessfully fulfilled. Two new universal param- 
eters are introduced. One is a strength pa- 
rameter which scales with mass, and ensures 
that stochastic reduction is negligible for mi- 
crosystems, but significant for macrosystems. 
The other is a localization length scale which 
defines the linear extent of the region to 



which stochastic reduction localizes an ex- 
panding wave-function. While known physi- 
cal and astrophysical processes put upper and 
lower bounds on these parameters, there is 
still a large permitted part of the parameter 
space and it will now be up to future labora- 
tory experiments to confirm or rule out these 
parameter values. 

Keeping in view the phenomenological na- 
ture of these models, which have been devised 
especially to resolve the quantum measure- 
ment problem, it is highly desirable to search 
for underlying physical principles and theo- 
ries for these models. Theories which emerge 
for reasons of their own, and which are not 
designed for the explicit purpose of explain- 
ing measurement. Trace Dynamics does well 
in this regard: its goal is to derive quan- 
tum theory from a deeper level, instead of 
arriving at quantum theory by 'quantizing' 
its own limiting case [classical dynamics]. It 
is an elegant structure in which Schrodinger 
evolution is the equilibrium thermodynamics 
of a 'gas' of classical matrices, and the ever- 
present Brownian motion fluctuations of the 
gas provide the stochastic process which com- 
petes with the equilibrium Schrodinger evo- 
lution. Under appropriate circumstances, the 
Brownian motion becomes important enough 
to be noticeable, and is responsible for the 
breakdown of quantum superposition. There 
perhaps could not be a more compelling 



representation of 'determinism + random- 
ness' than 'statistical equilibrium + statisti- 
cal fluctuations'. What is missing still are 
two important pieces of the puzzle: why 
do the Brownian motion fluctuations become 
more important for larger systems, and what 
is the origin of norm-preservation? 

Keeping Trace Dynamics aside for a mo- 
ment, one turns to investigate if gravity could 
couple with quantum effects and lead to an 
intrinsic uncertainty in space-time structure 
in such a way as to enable stochastic reduc- 
tion in macro-systems. At first glance, this 
seems not possible at all: quantum grav- 
itational effects can only be important at 
the Planck scale, and Planck length is too 
small to be of interest in laboratory physics, 
whereas Planck mass is too large to play 
a role in the quantum-classical transition. 
However, as more than one analysis shows, 
a subtle combination of linear extent of the 
object [measured in Planck units] and its 
mass [again measured in Planck units] allows 
gravity to bring about stochastic reduction. 
Gravity predicts the quantum-classical tran- 
sition very much in the domain in which it is 
expected on other grounds. 

Gravity provides a much needed physi- 
cal mechanism which could underlie sponta- 
neous collapse models. However, a proper 
mathematical treatment for building a grav- 
ity based theory of reduction is not yet avail- 



able. It is quite possible that a generaliza- 
tion of Trace Dynamics that includes gravity 
could unify spontaneous collapse and grav- 
ity models. Doing so could also explain why 
the Brownian fluctuations in Trace Dynamics 
and spontaneous collapse become larger for 
larger systems. For, we have indeed explic- 
itly seen in gravity models that the stochastic 
effect increases with mass. 

The need for inclusion of gravity in Trace 
Dynamics also stems from reasons having to 
do with space-time structure. With hind- 
sight, it is apparent that only when posi- 
tion localization is complete for nearly all ob- 
jects in the Universe, it becomes meaning- 
ful to talk of a background classical space- 
time geometry. If position localization is not 
achieved, and quantum coherence is signifi- 
cant, indeed that would prevent a meaningful 
definition of classical spacetime. Under such 
circumstances, and if one does not want to 
use classical physics starting point for 
quantization, one would have to include in 
trace dynamics, a matrix structure not only 
for the matter degrees of freedom, but also 
for space-time and gravity. Doing so holds 
the promise that one will be naturally led to a 
concrete mathematical formalism for describ- 
ing gravity induced collapse. Investigation 
and development of these ideas is currently 
in progress. 

A big stumbling block is the construction 



of relativistic models of spontaneous collapse. 
It is difficult to say at this stage whether this 
block will eventually be overcome, or is an in- 
dicator of some incompatibility between dy- 
namical models of wave-function collapse and 
special relativity. The collapse of the wave- 
function is an instantaneous process, and is 
said to violate the 'spirit' of relativity [like in 
an EPR experiment] . Radical though it may 
seem, we should eventually not be averse to a 
possible modification of special relativity to 
make it consistent with spontaneous collapse 
theories. Nevertheless, at the moment there 
is no matter-of-principle reason why collapse 
models should be incompatible with a fully 
relativistic scenario. Perhaps a generalized 
trace dynamics in which space and time are 
non-classical might have something useful to 
contribute here. 

The development of modified quantum 
theory has received great impetus from the 
arrival of pioneering experiments on molecule 
interferometry and optomechanics which can 
test these modifications. Prime amongst 
these is perhaps the 1999 discovery of in- 
terference and the verification of superposi- 
tion in the fullerene diffraction experiment. 
This paved the way for the developments that 
took place in the next two decades. Interfer- 
ence has now been observed in molecules with 
7,000 amu and tremendous effort is afoot 
to push this frontier to a million amu and 



beyond. Great ingenuity is being invested 
in devising new experimental techniques and 
technology which help advance this frontier. 
These experiments undoubtedly hold a place 
beside experiments which ushered in quan- 
tum theory a century ago: the spectrum of 
black-body radiation, atomic spectra, pho- 
toelectric effect, and matter interferometry 
with electrons. A broad class of theories 
predict that new physics will be seen in the 
range 10 6 amu to 10 9 amu. Perhaps in two 
decades from now, this range will have been 
tested. If quantum theory is found to hold 
good through this regime, then chances are 
good that linear quantum theory is univer- 
sally valid on all mass scales: we must then 
be content with many-worlds / Bohmian me- 
chanics, lest a more convincing interpreta- 
tion of the standard theory should emerge by 
then. If confirmation of the predicted modifi- 
cations is found, this will be nothing short of 
a revolution; a new theory of dynamics will 
have been born, to which quantum theory 
and classical mechanics will be approxima- 
tions. 
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Upper bounds on the collapse parameter A of the CSL model 
(with noise correlation length rc ~ 10 -5 cm) 



Laboratory 
Experiments 


Upper Bound 
on A 

(unit = s) 


Distance (orders of 
magnitude) from the 
CSL value 

Acsl ~ lO-^s- 1 


Distance (orders of 
magnitude) from 
Adler's value 

AAdler ~ 10 9 S 1 


Matter-wave 
interferometry experiments 


icr 5 


12 


4 


Decay of 
supercurrents 
(SQUIDS) 


lir 3 


14 


6 


Spontaneous 
X-ray emission 
from Ge 


icr 11 


6 


Excluded 


Proton 

decay 


10 


18 


10 


Cosmological 
Data 


Upper Bound 
on A 

(unit = s) 


Distance (orders of 
magnitude) from the 
CSL value 

Acsl ~ lCT 1 ^- 1 


Distance (orders of 
magnitude) from 
Adler's value 

AAdler ~ 10 9 S 1 


Dissociation 
of cosmic 
Hydrogen 


1 


17 


9 


Heating of 
intergalactic medium 
(IGM) 


10- 9 


8 





Heating of 
Interstellar dust 
grains 


10- 2 


15 


7 



TABLE I The table gives upper bounds on A from laboratory experiments and cosmological data, 
compared both with the CSL value A CSL ~ 10~~ s -1 and Adler's value A Ad ior ~ 10 -9 s _1 (see Section 
II. I). The X-ray emission bound excludes the value A A dicr for white noise, but this constraint is 
relaxed if the noise spectrum is cut off below 10 18 s _1 . Therefore, the bound coming from X-ray 
emission is very sensitive to the type of noise, i.e. to the type of collapse model. Large molecule 
diffraction would confront A CSL for molecules heavier than ~ 10 9 Daltons, and would confront A A di cr 
for molecular weights greater than ~ 10 5 Daltons. (The molecular diffraction bound on A decreases 
as the inverse square of the molecular weight, provided the molecular radius is less than rc; see 
Section II.H.) 
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